Google 



This is a digital copy of a book that was preserved for generations on Hbrary shelves before it was carefully scanned by Google as part of a project 

to make the world's books discoverable online. 

It has survived long enough for the copyright to expire and the book to enter the public domain. A public domain book is one that was never subject 

to copyright or whose legal copyright term has expired. Whether a book is in the public domain may vary country to country. Public domain books 

are our gateways to the past, representing a wealth of history, culture and knowledge that's often difficult to discover. 

Marks, notations and other maiginalia present in the original volume will appear in this file - a reminder of this book's long journey from the 

publisher to a library and finally to you. 

Usage guidelines 

Google is proud to partner with libraries to digitize public domain materials and make them widely accessible. Public domain books belong to the 
public and we are merely their custodians. Nevertheless, this work is expensive, so in order to keep providing this resource, we liave taken steps to 
prevent abuse by commercial parties, including placing technical restrictions on automated querying. 
We also ask that you: 

+ Make non-commercial use of the files We designed Google Book Search for use by individuals, and we request that you use these files for 
personal, non-commercial purposes. 

+ Refrain fivm automated querying Do not send automated queries of any sort to Google's system: If you are conducting research on machine 
translation, optical character recognition or other areas where access to a large amount of text is helpful, please contact us. We encourage the 
use of public domain materials for these purposes and may be able to help. 

+ Maintain attributionTht GoogXt "watermark" you see on each file is essential for informing people about this project and helping them find 
additional materials through Google Book Search. Please do not remove it. 

+ Keep it legal Whatever your use, remember that you are responsible for ensuring that what you are doing is legal. Do not assume that just 
because we believe a book is in the public domain for users in the United States, that the work is also in the public domain for users in other 
countries. Whether a book is still in copyright varies from country to country, and we can't offer guidance on whether any specific use of 
any specific book is allowed. Please do not assume that a book's appearance in Google Book Search means it can be used in any manner 
anywhere in the world. Copyright infringement liabili^ can be quite severe. 

About Google Book Search 

Google's mission is to organize the world's information and to make it universally accessible and useful. Google Book Search helps readers 
discover the world's books while helping authors and publishers reach new audiences. You can search through the full text of this book on the web 

at |http : //books . google . com/| 



WXoJUlo \ofe>=\. \'^ '^ 



HARVARD COLLEGE 



mm 



SCIENCE CENTER 
LIBRARY 



By Exchange 




^^"^ 



f 



t 



GENERAL PRINOIPLES 



OF THB 



METHOD OF LEAST SQUARES, 



WITH APPLICATIONS^ 



BT 



DANA P. BARTLETT, S.B., 

Pbofbsbob of Mathematx€», Masbachubiittb Inva'itutjb 

of tschnoloot. 



TUIBJ) EDITIOI^. 



BOSTON 
Tbb Authob 
1916. 



,% \5 ,S 




Xf/-jnr:ii<^<^ 



wt, 1910. 
Bi DANA P. BAHTLETT. 



PREFACE. 



The preparation of this volame was undertaken with the 
view of presenting in as simple and concise a manner as 
possible the fundamental principles of the Method of Least 
Squares. While it is believed that everything essential to 
the solution of all ordinary problems has been included, no 
attempt has been made to develop at length those special 
methods and forms that are so useful and almost necessary in 
case large numbers of observations of certain kinds, such, 
for instance, as those met with in geodetic and astronomical 
measurements, are to be adjusted. 

Frequent references throughout the text, and more particu- 
larly the list of works given on page v of the Appendix, will, 
however, enable the student to extend his studies in what- 
ever special direction his profession may require ; it being 
expected that this book will in such cases be looked upon 
merely as an introductory treatise. All of the works men- 
tioned have been freely consulted in the preparation of these 
pages, and the author desires in particular to acknowledge 
his indebtedness for many of the examples. 

Daka p. Babtuett. 
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THE METHOD OF LEAST SQUARES. 



CHAPTER I. 

OENSRAI. PRINOIPLB8. 

1. In scientific investigations of all kinds it is frequently 
necessary to determine the values of certain quantities by 
means of actual measurements either with or without the aid 
of instruments. The observations may be made directly upon 
the values of the unknown quantities or upon certain functions 
of the unknowns. In the latter case the values of the required 
quantities must be obtained by computation from the observed 
values of the functions. In order to obtain more accurate 
values of the unknowns than would be given by a single 
measurement, or set of measurements, the observations are 
usually repeated either in the same way and under the same 
conditions or in a variety of different ways and under vary- 
ing conditions. 

Under these circumstances it will invariably be found that 
the different measurements give discordant results, the amount 
of the discrepancies varying with the character of the observa- 
tions; and the question that now presents itself is how to 
determine from these discordant observations the true values 
of the required quantities. From the nature of the case, 
however, we can not expect to obtain our values with absolute 
accuracy; all that we can hope for is to obtain those values 
which are rendered most probable after all the observations 
are taken into account, and, further, to determine the degree 
of confidence that can be placed in those values. 
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2. . The attainment of the above results constitutes the 
primary object of the Method of Least Squares. The method 
is also employed in comparing the relative worth of different 
measurements of the same quantity, and in determining the 
equation of a curve which shall suitably represent the relation 
between two variables in cases where the exact law connecting 
them is not known. 

Also, before making any observations, we may employ the 
method to determine how precise the component measurements 
of a series must be in order to yield a required degree of pre- 
cision in the final result; or, conversely, to determine what 
the precision of the final result will be, knowing the precision 
attainable in the component measurements. This latter appli- 
cation of the method will be treated at length in the course 
on " The Precision of Measurements." 

3. Errors. The cause of the discrepancies between the 
results of our different observations is that every observation 
that is a measure is subject to error. These errors are of two 
kinds, — Constant or Systematic Errors and Accidental Errors. 

4. Constant Errors are errors which in all measures of the 
same quantity, made with the same care and under the same 
conditions, have the same magnitude, or whose presence and 
magnitude are due to some fixed cause. These constant 
errors may be of several classes, which are designated as 
follows: — 

First. Theoretical Errors, such as those due to the refrac- 
tion or aberration of light, the effect of a definite change in 
temperature or moisture on our standards of measurement, 
etc. As soon as their causes are known the magnitude of 
these errors may be calculated and their effect eliminated 
from the observations. 

Second. Instrumental Errors, such as errors of division of 
graduated scales, defects in micrometer screws, eccentricity of 
circles, etc. These errors will be discovered by an examina- 
tion of the instruments and their effects eliminated from the 
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observations, either by a particular method of using the 
instruments or by subsequent computation. 

Third. Personal Errors. These are due to personal pecu- 
liarities of an observer, who always answers a signal too soon 
or too late, always estimates a quantity smaller than it is, etc* 
The character and magnitude of these errors may be deter- 
mined by a study of the observer, his "Personal Equation" 
may be obtained, and his observations thus corrected for this 
source of error. 

5. Mistakes. Although of a somewhat different character, 
these should be considered in connection with constant errors. 
A mistake is made when a figure 3 is read for a %ure 8, or 
when in reading a graduated circle which is numbered in both 
directions the angle is read 43° instead of the complementary 
angle 47°, etc. These mistakes are usually of such a charac- 
ter that they may be detected by an inspection of the observa- 
tions and a proper correction made. 

6. Accidental Errors are errors due to irregular causes, 
whose effect upon the observations is not determined by any 
circumstances peculiar to that particular set of measurements, 
and which cannot therefore be computed and allowed for 
beforehand. Such errors are those due to sudden changes in 
refraction owing to sudden and unobserved changes in tem- 
perature; unequal expansion of different parts of an instru- 
ment with change in temperature; shaking of an instrument 
in the wind, etc. But most important of all are those errors 
which arise from imperfections in the sight, hearing, and other 
senses of the observer, which render it impossible for him to 
adjust and use his instruments with absolute accuracy. 

After a full investigation of the constant errors, the observer 
should diminish the accidental errors as much as possible, both 
in number and magnitude, by taking every precaution and care 
in the measurements themselves. The problem now remains 
to combine the observations so that the remaining accidental 
errors shall have the least probable effect upon the results, 
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and it is to bring about this combination of observations that 
we employ the Method of Least Squares. 

When no more observations are made than are sufficient to 
determine one value for each of the unknown quantities, we 
must accept these values as the most probable ones. But if 
additional observations are made leading to discordant results, 
we can not take any one of them as the correct value, and in 
fact, as already stated, we shall probably not be able to obtain 
the true values of the unknowns. All that we can do is to find 
values of the unknowns which shall remove the discrepancies 
between the different observations and which shall be those 
values that are rendered most probable by the existence of 
the observations themselves. 

On first thoughts it may seem that these accidental errors, 
being due to so many different and unknown causes, will be 
beyond the scope of mathematical investigation. Neverthe- 
less, the theory of probability requires that these errors shall 
follow in magnitude and frequency a law that is capable of 
'vc;^ exact mathematical expression, and experience confirms the 

^ J correctness of this law. 






6 For more extended remarks on these subjects see — 
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v^ 1 5 Bolman, ** Discussion of the Precision of Measurements," pp. 1-14. 
Merriman, ** Text-Book of Least Squares," pp. 1-6. 
Chauvenet, '* Spherical and Practical Astronomy," pp. 469-473. 
Wright, ''Treatise on the Adjustment of Observations," pp. 11-18. 
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^ LAWS OF ERRORS OF OBSERVATION. 



7. The derivation of the general laws of the occurrence of 
errors of observation, and of the processes for determining 
the most probable values of the unknown quantities, will be 
based upon the following 

Axiom. If a series of n direct observations, Miy M^ . . . M^<^ 
is made upon the value of a quantity M^ all the observations 
being made with the same care and under the same circum- 
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Stances, the most probable value M^ of that quantity is the 
arithmetical mean of the observations. Or 

® n n 

8. The Beal Error {x) of an observation is the difference 
between the observed value of the measured quantity and the 
real value. 

9. The Besidual {v) of an observation is the difference 
between the observed value of the measured quantity and 
the value rendered most probable by the existence of the 
observations. 

10. Example, Eight observations are made upon the 
resistance of a coil of wire, the true resistance being 512. 
Find from these observations the most probable resistanoei 
and also the real errors and residuals. 



Observations. 


Beal Errors. 


Besldaals. 


M. 


X 


V 


512.4 


+ •4 


+ .30 


612.2 


+ .2 


-f .10 


611.9 


-.1 


-.20 


612.3 


+ .3 


+ .20 


611.8 


-.2 


-.30 


512.3 


+ .3 


+ .20 


611.9 


-.1 


-.20 


612.0 


.0 


-.10 


Mean — 512.10 




S« — .00 



Frcnn the observations, then, we should say that the most 
probable resistance of the coil is 512.10. It will also be 
noticed that the sum of the residuals is zero. That this is a 
general result following from the assumption of the arith- 
metical mean as the most probable value may be proved as 
follows: If the observations are J^, Jf^, . . . M^t the 
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arithmetical mean M^y and the residuals Vi, v^, . . . v^y then 
we have 

= 2ilf— 2^ since J!^= 

n 

.-. Sv=0 (2) 

11. Weighted Observations. The weight of an observa- 
tion expresses its relative worth compared with other obser- 
vations. Thus, if six observations are made upon the value 
of a quantity, five of which give the same result, while the 
sixth differs, in combining these two different results' to 
obtain the most probable value of the unknown, the first 
value ought to have five times the influence upon the final 
result that the second has, since it has taken five times as 
much labor and time to obtain it. Hence in general we may 
say — 

12. The Weight {p) of an observation may be considered 
as representing the number of times the observation has been 
repeated and the same result obtained. 

The weights assigned to observations may be due to a 
variety of causes, as difference in skill of observers, difference 
in the instruments used or the circumstances under which 
the observations are made, etc. But whatever the cause, the 
effect on the final values of weighting an observation will be 
the same as indicated in the 'preceding paragraph. 

13. Example, Suppose n observations, JIfi, J!^, . . . jS^y 
of weights Ply />s» • • • Pm ^^^ made upon the value of a 
quantity M. To find the most probable value Mq of the 
quantity. 

From the above interpretation of the meaning of weight, 
we may consider that the whole number of observations is 
Pi'{'Pi'\- • • •jPm or %py and that the result Mi has been 
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obtained in pi observations, M^ in p^ observations, etc.. 
Therefore, by (1) 

Ji^ is called the General Mean. 

If the residuals are v^ ^;s, . . . Vn, we have 

2:i>v=0 (4) 

Which shows that in the case of direct observations of differ- 
ent weights the sum of the weighted residuals is zero. 

14. If the observations are not made directly upon the 
values of the required quantities, the method of adjusting 
the results so as to obtain the best possible values of the 
unknowns will depend upon the laws which govern the dis- 
tribution of the errors of these observations. It is found in 
practice that the accidental errors of observations follow cer- 
tain well defined laws, and what these are may best be seen 
by taking an actual example. 

15. Example, One thousand shots are fired at a target 
which is divided into a number of horizontal sections by lines 
one foot apart, the centre line of the target being in the 
middle of one of these spaces. The shots were distributed 
as follows : — 



At Spaet. Shot$. 


In Spaet. 


atoti. 


InSpaee, 


Shot 


-f-5Jto+4j- 1 


+Hto+ i 


190 


-2ito-8i 


79 


+4*" +3* 4 


+ * " - i 


212 


-H « -4i 


16 


_).3J " +2i 10 


- i " -H 


204 


-4j- " -5i 


2 


■^^ « -\-H 89 


-H " -2* 


198 







In this case the errors are evidently the distances of the 
shots from the centre of the target. Further, as far as can 
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be jailed from these one thousand shots, if another ahot U 
fired the probability that this shot irill fall between the 
lines 

+ 5Jand+4J is .001 

4-4J " +3i " .004 

4-3J « -I-2J " .010 

+ 2* " +H " -089 

+ li « + i«.190 

+ i « - i " .212 
The sam of the above probabilities is unity, and, therefore, 
as far as the preceding shots show the lOOlet shot viU cer- 
tainly hit the target. 

16. Now using as abscissas tbe distances of the horizontal 
lines from the centre of the target, and as ordinates the num- 
ber of shots falling in the oorresponding spaces, we may 
oonstmct the following figure : — 



- J and 


-11 i. 


.204 


-U " 


-2i' 


.193 


-^ » 


-S}. 


.079 


-3J •• 


-4i' 


.016 


-4i " 


-H« 


.002 



And if the entire area of tMs figure is taken as unity then 
ha area of each rectangle will denote the probability of a 
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9 



shot, if fired, falling within the corresponding space of the 
target. 

The graphical representation of the accidental errors of 
observation will always give a figure similar to the above. 
Hence denoting errors by abscissas, and their frequency by 
ordinates, the law of error of any series of observations may 
be represented by a curve whose general form is determined 
by Figure 1. This curve is called the " Curve of Error," and 
is shown in Figure 2. 

r 



Figure 2. 




In order that this curve may represent exactly the distri- 
bution of the errors in any given series of observations it 
ought to meet the axis of JT at some definite distance to the 
right and left of the origin and coincide with the axis from 
there on, for in all actual observations there is a limit beyond 
which no errors occur. But as the exact point of meeting 
could not be determined for any given case, and as it would 
not b^ possible to obtain the equation of such a curve, we 
make it asymptotic to the axis of J^ taking care that the 
error thus introduced shall in any set of observations be so 
small as to be negligible. 
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17. ' An inspection of Figures 1 and 2 will now exhibit 
some of the general lawc of errors of observation and the 
corresponding properties of the curve of error. 



Laws of Error derived from an 
inspecHon of Figure 1. 



First, Small errors are more 
frequent than large ones. 

Second. Positive and negative 
errors of the same absolute mag- 
nitude are equally likely to occur. 

Third, The probability of the 
occurrence of very large errors is 
very small. 

Fourth, The frequency of any 
error depends upon the magnitude 
of that error. 



Bepresentation of these laws by the 
Curve of Error, 



The maximum point of the 
curve is on the axis of F. 

The curve is symmetrical with 
respect to the axis of F. 

The curve is as3rmptotic to the 
axis of X. 

The equation of the curve will 
be of the form 



y=4>{x) 



(5) 



18. If, now, the total area between the curve and the axis 
of JT be denoted by unity the probability that the error of 
any given observation will fall between the magnitudes x and 
x-\-dx will be represented by the area included between the 
curve, the axis of J^ and the ordinates of the curve at the 
errors x and x -|- docj or by 



y dx=<f>(^x) dx 



(6) 



And this probability will be known as soon as we find the 
form of the function 4>(x), 

19. The above expressions in (5) and (6) are the ones 
that we should use if we regard the curve of error as repre- 
senting the law of occurrence of errors of observation. If, 
however, we look upon the curve as expressing the law to 
which we must make the residuals conform, in order that the 
values of the unknown quantities obtained from them may be 
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the most probable valaes, we should. replace x hj v and use 
the expressions 

y=*(^) (7) 

and y€lv=<l>(v)dv (8) 

for (5) and (6), respectively. 

THE EQUATION OF THE CURVE OF ERROR. 

20. Let n observations, all of the same weight, with 
results Mi^ M^ . . • it^y ^ made upon any function or 
functions of a number of unknown quantities 2i, z^y • . . Zq\ 
and let the residuals of J!^, JS^, . . . M^ be v^^ V29 • • • ^n) 
and the probability of the occurrence of these residuals be 
^(Vi) dv^ ^(^2) dvy . . . ^(Vn) dv^ respectively. Then the 
probability of the simultaneous occurrence of all these 
residuals will be 

P=4>{v{) 4>{v^) . . . 4>(vn) (dv)^ (9) 

.% log P= log <l>(Vi) -|- . . . log <l>(Vn) -J- n log dv (10) 

Each different method that might be adopted for computing 
the values of the unknowns z^ s;^, . . . Zg would lead to a dif- 
ferent set of residuals Vi, t's, • . . f n » l>ut obviously that set 
of values of z^ z^^ , , , Zq should be considered the best which 
corresponds to the particular set of residuals Vi^ Vg, • • . v^y 
the probability of whose occurrence is greater than that of 
any other set. 

Therefore the most probable values of z^ z^y , . , Zq are 
those that make I^ in (9), or logP in (10), a maximum. 

The values of 2^, ^29 • • • Zq corresponding with this latter 
condition are those that satisfy equations (11). It maybe 
noticed that these equations also express the preliminary 
conditions leading to a minimum value of logP, but the 
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nature of the problem is such that a maximum value of JP 
evidently exists while a minimum does not, and it is there- 
fore unnecessary to investigate further the mathematical 
conditions for a maximum. Hence we have 



aiogP 1 d<t>(Vi) 



1 BH^n) ^^ 



dz^ *(vi) a^a ^^ .., . 



• • • 



aiogP 1 a^(vi) 



dZq 



1 ^*(vi) . _J_iiK)=o 



i^ = *.w£ (12) 

and for convenience we may put 

substituting in (11) we have 

*(f i) aij + *(«.) 5^ + • • • 'A(»«) aif = 






(14) 



dvi dVi dVf^ 

H^i) ai^ + *(^0 ar, + • • • *("•) ai, = <> 

These equations contain all the unknowns 2^, ^2y > • ^qy ^^^ 
there are as many equations as unknowns, hence as soon as 
we find the form of the function il/(v) we can solve these 
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equations for the most probable values of 2i, z^^ , , , Zq. Since 
we have considered the general case, and the above results 
are to hold true whatever the number of unknown quantities 
and the form of the functions observed, we may deduce the 
form of ^(v) by solving a special example. 

Example. Let n observations of equal weight be made 
upon the value of a single unknown e^, with results M^ M^^ 
. . . Mf^y and let the residuals be Vi, Vj, . . . v^. Then the 
most probable value of Zi is given by 

differentiating with respect to ^i, 

1=——=——= — ^ (\ 

Bzi dz^ * * * Bzx ^ ^ 

substituting (a) in (14), changing all the signs, we have 

*(^i) + H'^2) + • • • K'^n) = (b) 

But in this case, as was shown in (2), 

Vi + Vj + . . . v„ = (c) 

In order that (b) and (c) may both be true the functional 

symbol ^ must indicate multiplication by a constant. That 

is, in general 

\l/(v) = CO (15) 

Substituting this in (13) and (12), 

d<l>(v) , ^ dv 

^ ^ 1 ^*(v) ^v 

therefore , . x ;r — = cv ^r— 

4>{v) az az 

Integrating, log <l>(v) = icv^ + ^ 
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Since y==^(v) is the equation of the curve of error, (7), 
we may therefore write it 

But on examination of the curve, y is seen to be a decreas- 
ing function of v, and hence the exponent of e is essentially 
negative. Accordingly we will write our equation in the 
form 

y=he-h'^ (16) 

the values of k and h depending upon the character of the 
observations, but in all observations of the same kind and 
weight having the same values. 

This equation represents the law in accordance with which 
the residuals must be distributed in order that the best results 
may be obtained from our observations. But, as before 
mentioned, if we wish our curve to represent the most prob- 
able distribution of the real errors of observation we should 
write the equation in the form 

y=he-h^^ (17) 

Hereafter we shall use without further remark either form 
of the equation according to the aspect in which we are 
considering our curve. 

An inspection of the above equation will show that it 
satisfies all the conditions noted in discussing the form of the 
curve of error in paragraph 17. 

21. It is important to notice that in all discussions in the 
Method of Least Squares the number of observations is sup- 
posed to be large and always greater than the number of 
unknown quantities. As will be illustrated later on, para- 
graph 91, whenever this is the case there is a remarkable 
agreement between the results obtained in practice and those 
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indicated by the theory. And even when the observations 
are few in number the method still affords the best means at 
our command for their adjustment, the results obtained 
merely having a smaller weight than they would have had if 
derived from a greater number of observations. 

THE METHOD OF LEAST SQUARES. 

22. We are now in a position to see whence comes the 
name " Least Squares." 

In paragraph 20 it was pointed out that whenever we 
make a series of observations, each observation of the set 
having the same weight, the most probable system of values 
of the unknown quantities will be that which corresponds 
with the set of residuals the probability of whose occurrence 
is a maximum. That is, the best set of values of the 
unknowns will be that which gives a maximum value to 

But from equation (16) this reduces to 

Since the exponent of e in this expression is negative, 
evidently P will be a maximum when 

Vi^+t^a*+ • • • Vn^=^v^ is a mind/mum. (19) 

Hence the adjustment of observations by the Method of 
Least Squares is based upon the principle that the most prob- 
able system of values of the unknowns is that which renders 
the sum of the squares of the residuals a minimum. Hence 
the name. 

The conditions for a maximum value of P were expressed 
in equations (14), and since it has been shown that the 
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function ^ means multiplication by a constant, those equa- 
tions reduce to the following, called 

NORMAL EQUATIONS. 






(20) 



• • • 

• • • 



Bvi , „ dv, , .. dv. 



*'^+*«5^ + "-^-§^=« 



An inspection of these equations will show that they also 
express the conditions that will make the sum of the squares 
of the residuals a minimum. 

In the adjustment of observations the above are the funda- 
mental equations. In order to obtain the most probable 
values of the unknowns in any set of observations, all that is 
necessary is to form the Normal Equations for that set and 
solve them simultaneously. The examples already solved for 
direct observations are merely special cases of the above 
general solution. 



CHAPTER 11. 

THS Al>Jn8TMENT OF OBSERVATIONS. 

INDIRECT OBSERVATIONS. 

23. In the determination of the values of quantities by 
means of observations the functions of the unknowns that it 
is necessary to observe may be of any form, but if they are 
not linear the normal equations derived from them are likely 
to be complicated and difficult, if not impossible, to solve. 
Hence if the observations are not upon linear functions of 
the unknowns, the first step will be to reduce them to 
equivalent linear expressions by transformations depending 
on the character of the observed functions ; see paragraphs 
43 and 44. It will be necessary to consider, therefore, the 
method of adjusting observations on linear functions alone, 
and the procedure in cases of this kind may be illustrated by 
the following simple example. 

24. Example, S^^ /S^, /S^, are three solids whose masses 
are required. Not having standard weights enough to obtain 
all these masses directly, by varying the distribution of the 
solids in the pans of the balance the following observations 
are made : — 

aS'i = aS'j -f- 1'7 grams. 
8^ = 2.4 « 

A^a + /^8 = /^i + 1.0 « 
/^2 = /^3 4- 3.0 « 

If, now, the most probable values of /^j, jS^ /Sj, are repre- 
sented by 2^1, Z2y ^89 ^^^ the corresponding residuals of the 
observations by Viy Vj, . . . v^^ the values of these residuals in 
terms of «i, «2, z^ may be found from the above observations 
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by transposing all terms to the first members of the eqna^ 
tions, and we obtain at once the following, called 

OBSERVATION EQUATIONS. 



5zi — 622 — 2.5 = 

2 X (b) — 42i + 6za — 4.6 = 



(A) 



J5j — 2J2 — 1.7 = Vi 

gg — 2.4 = Vi 

— 2l + 2^2 -f. 28 — 1.0 = Vg 
^2 — 2^8 — 3.0 = V4 

Applying equations (20) as formulas we have the 

NORMAL EQUATIONS. 
(^ - 22 - 1.7) + (- 2Ji + 22 + «8 - 1.0)(- 1) = 

(2^^2,-1.7)(_l) + («2Ji + 22+28-1.0) 

+ (22-«8-3.0) = (B) 

(28-2.4)+(-2i+Z2+»8-1.0) 

-|-(«,-2J8~3.0)(-l) = 

Simplified, these become 

221 — 222— 28 — 0.7=0 (a) 

— 22i + 322 — 2.3 = (b) 

— 2i -[-328—0.4=0 (c) 

3 X (a) 621 — 622 — 328 — 2.1 = 

(c) — 2i -1-328 — 0.4=0 



«1 


— 7.1 — 






2i— 7.1 


(d) 


substitute (d) in (b) • 


22 5.5 


(e) 


« (d) in (c) 


28 — 2.5 


(f) 



An inspection of the work in this example will show that 
for the adjustment of observations of equal weight on linear 
functions of the unknowns we may derive the following : — 
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25. BiQe. For each observation write an " Observation 
Equation'*'^ ; then for each unknown form a ^^ Normal 
Equation^^"* by multiplying the first member of each observor 
tion equation by the coefficient of that unknown in that 
eqtiation^ adding the results and placing the sum equal to 
zero. Solve these equations simultaneously for the values 
of the unknowns. 

In solving for the most probable values of the unknowns 
the second members of the observation equations are very 
commonly written zero instead of Vi, v^ , . . v„. For this is 
the form in which the equations naturally appear, and if the 
observations were exact the residuals would actually all be 
zero. The method of solution is the same in either case. 

26. Observations of Unequal Weight. If the observations 
are not all of equal weight the same method will apply, 
except that in the formation of the normal equations each 
observation equation will be used the number of times denoted 
by its weight. Thus in the last example if the observations 
have the weights 4, 9, 1, 4, the normal equations will have 
the same form as in (B), page 18, but each part of each 
equation will be multiplied by the weight of the observation 
equation from which it is derived. This will give the 

NORMAL EQUATIONS. 

4(Zi-«,-1.7) + (-% + aj3 + 23-.1.0)(-l) = 

4(2i - «, - 1.7)(- 1) + (- 2i + 2, + 23 - 1.0) 

+ 4(2a-28 — 3.0) = 
9(2, ^ 2.4) + (- 21 + 2,4-23 - 1.0) 

+ 4(2a-2,-3.0)(-l) = 
or reduced 

52i — 62a — 2, — 5;8 = (a) 

— 521 + 922— 328— 6.2 = (b) 

— 2i — 32a +142,— 10.6=0 (o) 

the solution of which gives 

2i = 7.07 2a = 5.42 2, = 2.42 
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Further it will at once be seen that if jOi, J»2> • • • Pw ^^^ the 
weights of the corresponding observations, equations (20) 
take the general form: — 

WEIGHTED NORMAL EQUATIONS. 






(21) 



• • • 

• • • 



V 



Hence for the formation of the normal equations in weighted 
observations on linear functions of the unknowns, we have 
the following: — 

27. Bide. For each observation write an ^^ Observation 
Equation'*'* ; then for each unknown form a ^^Normal 
Equation^'* by multiplying the first member of each observa- 
tion equation by the coefficient of that unknown in that 
equation and by the weight of that equation, adding the 
results and placing the sum, equal to zero, Solve these 
equations simultaneously. 

28. The same result will be obtained if we begin by 
multiplying each observation equation by the square root of 
its weight and then proceed according to the first rule 
(paragraph 25). 

This result illustrates the important principle that multi- 
plying a set of equations by the square roots of their weights 
reduces them, all to equivalent equations of weight unity. 

29. Belation between the Weight of an Observation and 
the Value of h. If in paragraph 20 the n observations 
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have weights p^^ jOg, . . . Pw and the quantity h values 
Aj, Aj, . . . A,j, then equation (18) becomes 

= kjc^ . . . A:„6-(^i'^i' + Wf^2« + . . . An^f^n^) (ef,,)« (22) 

The most probable set of values of the unknowns is that 
which makes P a maximum, and P is a maximum when 

hi^ Vi + /ih* t?a^ + • • • hn Vn^ is a minimum. (25) 

The conditions for a minimum value of this expression are 
the following, which are then for this case the 

NORMAL EQUATIONS. 






(24) 



• • • • • • 

• • • • • • 



But equations (21) are also the normal equations for 
this case. Hence (21) and (24) must be identical, and 

Pi : P2 : . • • jPn = fei^ : th^ : * - ^W (25) 

That is, the square of h is proportional to the weight of 
the observation. Accordingly, since h increases in value as 
the quality of the observations is improved, it is called " The 
Measure of Precision." 
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Further, it follows from (23) and (25) that the most 
probable system of values of the unknowns will be that in 
which 

jPit^i*+i>2t^2*+ • • 'PnVn^ i8 a minimum. (26) 

And this is the most general form of statement of the 
principle of " Least Squares." The same principle is repre- 
sented in equations (21). 

30. Computation of Correotions. If large numbers occur 
in the observations it is better to compute the most probable 
corrections to apply to the observed values rather than the 
most probable values of the unknowns themselves. In this 
way we can often avoid a large amount of numerical work. 

31. JE^xample. Pi, Pgi ^sy -^4* ^sy are five points whose 
altitudes above the mean level of the sea are to be determined 
from the following observations of difference of level. 

Pi = 573.08 P4 - Pa = 170.28 

Pa - Pi = 2.60 P4 — Pg = 425.00 

Pa = 575.27 P6 = 319.91 

P, — Pa = 167.33 Pg = 319.75 
P4-P,= 3.80 

An inspection of these observations shows that we may put 

Pi = 573 + 2i P4=745 + 2;4 

Pa = 575 4-22 1^6 = 320 -f z^ (A) 

P3=742 + 23 

where z^ ^29 ^89 ^49 ^59 ^^^ small corrections whose most 
probable values are to be determined. We now have for 

OBSERVATION EQUATIONS. 

573 + 2i — 573.08 = or 2i - .08 = 

6764-2a — 573 — 2i— 2.60=0 or 2a — ^i — . 60 = 

575 + 2a - 575.27 =0 or 2, — -27 = 



or 


28 


-2a -.33 = 


or 


24 


-28— .80 — 


or 


24 


-2a -.28 — 


or 


24. 


— 26 =0 


or 




26+. 09=0 


or 




26 + . 26 = 
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742 + 2, - 575 - 23 - 167.33 = 
7^^^z^^742-z^-' 3.80=0 
745 + 24 — 575 - 2a - 170.28 = 
745 + 24 — 320 - 26 - 425.00 = 

320 + 26-319.91 = 

320 + 26-319.75 = 

From these we now form the 

NORMAL EQUATIONS. 

22i- 2a + .52=0 

— 2i + 422— 28—^4 — .26=0 

— 2a + 228- 24 + .47=0 

— 2a— 28 + 324— 26—1.08=0 

— 24+326+ .34 = 

and solving, 

2i= — .19; 2a=:.14; 28 = . 05; 24=. 43; 25 = . 03 

Substituting these in equations (A) we have for the most 
probable altitudes, 

Pi = 572.81 P8 = 742.05 ^6 = 320.03 
Pa = 575.14 P4 = 745.43 

If the original observation equations had been retained, 
the independent terms in the normal equations would have 
been 

570.48 240.26 163.53 599.08 214.66 

32. Significant Figures. The adjustment by the Method 
of Least Squares of observations which occur in practice, 
although not difficult, is apt to be long and laborious. Hence 
to reduce this labor as much as possible it is of great import- 
ance that careful attention should be given in the solutions 
to the proper use of significant figures. When in doubly 
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however, as to the proper number of figures to retain it is 
better to keep too many rather than too few, as the superfluous 
figures can be rejected at the end of the computation ; while 
if too few are retained the results obtained from the compu- 
tations will be worthless. 

For a general discussion of the subject of significant figures 
see Holman's << Precision of Measurements," pages 76 to 84, 
but for the present the following rules will suffice for most 
cases. 

Bnle 1. Jn casting off places of figures increase by 1 the 
last figure retained^ when the following figure is 6 or over, 

Bnle 2. In the precision measure retain two significant 
figures. 

Bnle 3. In any quantity retain enough significant figures 
to include the place in which the second significant figure 
of its precision measure occurs. 

Bnle 4. When several quantities are to he added or 
subtracted^ apply Mule 8 to the least precise and keep only 
the corresponding figures in the other quantities, 

Bnle 5. When several quantities are to be multiplied or 
divided into each other^ find the percentage precision of the 
least precise. If this is 

1 per cent or more^ use four significant figures, 
.1 " " " « " five " " 

.01 « " « « « six " " 

in all the work. If the final result obtained in this way 
conflicts with Hule 5, apply the latter, 

Bnle 6. When logarithms are usedy retain as many 
places in the mantisscB as there are significant figures 
retained in the data under Hule 5, 

The application of these rules is not always possible in the 
course of the work, since the precision measures may not be 
known until the end of the computation. But as a general 
rule it is sufficient in direct observations to retain one more 
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place of figures than is given by the individual observations, 
and in indirect observations to retain two additional places. 

CONDITIONED OBSERVATIONS. 

33. Conditioned Observations are those in which the 
unknown quantities must be determined not only so as to 
satisfy as closely as possible the observation equations, but 
also so as to satisfy exactly certain other conditions. These 
conditions must be less in number than the unknown quantities, 
otherwise the unknowns could be determined from the con- 
ditions alone. 

The adjustment of observations of this class may be reduced 
to the method already used for unconditioned observations in 
the following way. 

The observations are represented by " Observation Equa- 
tions," and the conditions by certain other equations, called 
" Condition Equations." 

Between these two sets of equations we will eliminate as 
many unknowns as there are conditions. From the resulting 
equations, which will be the same in number as the observa- 
tions, we will form in the usual manner the ^^ Normal 
Equations" for the remaining imknowns. Having solved 
these normal equations and substituted the results in the con- 
dition equations, we shall obtain the values of the unknowns 
first eliminated. 

All conditions of the problem are now fulfilled, for the 
condition equations are satisfied exactly and, moreover, 
according to the principle of Least Squares our results 
are those rendered most probable by the existence of the 
observations. 

As in the example last considered, it is often more 
advantageous to compute corrections to the observed values 
of the unknown quantities rather than the values of the 
quantities themselves. 



"-V 



( 
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34. Example, Find the most probable values of the angles 
of a quadrilateral from the observations. 



A — 101^ 13' 22" weight 3 
J9= 93 49 17 « 2 
C= 87 5 39 " 2 
2>= 77 52 40 " 1 



(A) 



JL+J94-(7+2> = 360^ 0' 58" 



The condition to be satisfied is in this problem 

u4 + ^ + (74-2>=360^ (B) 

Let Ziy z^ 089 ^4 ^6 ^^6 most probable corrections to add 
to the observed values. This gives for 



OBSERVATION EQUATIONS 

2i=0 weight 3 
a,= « 2 
2g=0 « 2 

24=0 « 1 



(C) 



and for the 



CONDITION EQUATION 

2i + «.+«8 + «4 + 58 = (D) 

Eliminating 24 between (D) and (C), the equations from 
which the normal equations are to be derived become 

«! = weight 3 

«,= " 2 J, 

2,= "2 ^ ^ 

J5j4.2,-|_g,-|_58=0 « 1 
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Applying the rule in paragraph 27, these give the 

NORMAL EQUATIONS 

42l+ 22+ 28+58=0 
21 + 322+ 28+58=0 (F) 

2i+ 22 + 328+58=0 

Solving, and substituting the results in equation (D), we 
find 

2i=- 8.29 28 =-12.43 .Qv 

2j=- 12.43 24= -24.86 ^ ^ 

Applying these corrections to the observations (A), the 
most probable values of the angles are 

A = 101^ 13' 13".71 

J9= 93 49 4.57 .^v 

(7= 87 6 26.67 ^ ^ 

D= 77 52 15.15 

Vote. Jn eliminating unknowns between the observation 
and condition equations care must be taken that the obser- 
vation equations are not combined with each other or multi' 
plied by any quantity. For if this is done the weights of 
the observation equations will be altered, (See §28.) 

35. In the above example it is evident that the corrections 
to be applied to the different observations are inversely as 
their weights. And, in general, when there is but one 
equation of condition, the observations expressing direct 
determinations of the unknowns, the corrections will be pro- 
portional to the coefficients of the unknowns in the equation 
of condition divided by the weights of the corresponding 
observations. A proof of this is given in paragraph 117. 
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The most common case is that in which these coefficients 
are all unity, as in the example just solved, and we may then 
derive the 

BiQe. Find the difference between the theoretical and 
observed results and divide this correction among the 
observations in the inverse ratio of their weights. 

In the last example the sum of the observed angles exceeds 
860** by 68". Therefore the correction to be applied to JL is 

1 
" ^^ >< \ r"^"i = - 58 X i = - 8.29 

EMPIRICAL FORMULAS AND CONSTANTS. 

36. In the work so far considered the observations are 
supposed to be made either directly upon the values of the 
unknown quantities or upon some function of the unknowns 
whose form, and the constants entering into it, are definitely 
known. But another sort of problem frequently occurs, 
in which observations are made upon the values of a certain 
variable and the corresponding values of some function of it, 
the exact form of the function not being known. The object 
in this case is the determination of the most probable form of 
the function and the values of the constants involved ; that 
is, the derivation of the algebraic expression best representing 
the law connecting the variable and function. 

This expression may be looked upon as the equation of a 
curve, abscissas denoting values of the variable and ordinates 
values of the function, and for all values of the variable 
within the range of the observations we may determine from 
it the most probable values of the function corresponding. 
But except in special cases, where the number of observations 
is large, where the law connecting variable and function is 
well defined, and where the equation obtained is an accurate 
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representation of this law, it cannot be assumed to apply 
beyond the range of the observations. And in no case would 
it be safe to make use of the curve very far beyond the limits 
of the observations. 

37. The Method of Least Squares will not assist in deter- 
mining the form of the function. This must be settled upon 
beforehand, either from theoretical considerations or by 
constructing a plot, using values of the variable as abscissas 
and of the function as ordinates, when the smooth curve 
drawn through the points thus obtained will indicate the 
form of equation to be used. 

It is to be observed that this is a method of trial and will 
not necessarily give the most probable form of the function ; 
and in fact we may not be able to obtain the form that would 
be absolutely best. Further, several forms of equation may 
be known which would represent well the plotted points. In 
such a case that should be considered the best in which the 
sum of the squares of the residuals is found to be the least. 

38. As soon as the form of the function is decided upon it 
should be reduced to the linear form, and the determination 
of the values of the constants involved is then a simple 
application of the preceding methods. 

As the " Observation Equations " in any given problem will 
all be of the same kind, it is usually advisable to write out 
the general form of the "Normal Equations" and arrange 
the computations in tabular form, while the retention of 
the proper number of significant figures is of particular 
importance in this work. 

39. A case that frequently occurs is that in which the 
quantity y is a constantly increasing function of the vari- 
able aj, or where the plotted curve is approximately parabolic 
in form. Here the equation 

y = A-^Bx+ Cx^+Dx'-^- . . . (27) 

may be taken to represent the relation between the variable 
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and the function. The larger the number of terms taken in 
the second member, the more accurately may the equation 
obtained be made to represent the results of the observations ; 
but the labor involved increases rapidly with increase in the 
number of terms, and if the plot shows a very nearly straight 
line the first two terms alone may suffice. 

40. Mxample, In measuring the velocity of the current 
of a river the following results were obtained :— 



Depths. 

X 

1 


Yelocities. 
V 

4.86 


2 


6.14 


3 
4 


5.16 

4.86 


5 


4.24 


6 


3.36 


7 


2.16 


8 


0.67 



The velocity at the surface is 4.260. Find the equation of 
a curve which will express the relation between x and V. 
Plotting the observations we find the curve 

V 
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This is approximately parabolic in form and passes through 
the fixed point (0, 4.25). Therefore the relation between x 
and V may be expressed by the equation 



F= 4.26 + ^x4- (7ar» 



(A) 



and substituting in this the corresponding values of x and V 
as given by the observations, we shall have eight observation 
equations from which the most probable values of B and C 
are to be computed. 

All of the observation equations being of the form (A) we 
have the 

NORMAL EQUATIONS 



0^ Saj*+ ^ 2aJ« + 4.25 SiB" - S T^= 



(a) 
(b) 



For computing the coefficients in these equations it is most 
convenient to arrange the following table. 



X 


V 


Vx 


x" 


Vx^ 


a;' 


X* 


1 


4.86 


4.86 


1 


4.86 


1 


1 


2 


5^.14 


10.28 


4 


20.56 


8 


16 


3 


5.15 


15.45 


9 


46.35 


27 


81 


4 


4.85 


19.40 


16 


77.60 


64 


256 


5 


4.24 


21.20 


25 


106.00 


125 


625 


6 


3.36 


20.16 


36 


120.96 


216 


1296 


7 


2.16 


15.12 


49 


105.84 


343 


2401 


8 


0.67 


5.36 


64 


42.88 


512 


4096 


36 




111.83 


204 


525.05 


1296 


8772 


Sa; 




SFa; 


Sar» 


tVx" 


So* 


SaJ* 
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Substitnting these results in (a) and (b) we have 

8772 (7-f 1296 j5 + 341.96 = (c) 

1296(7+ 204J5+ 41.17 = (d) 

and solving, 

(7 = - .1493 B = .7465 (e) 

Therefore the required equation is 

V = 4.25 + .7465aj - .1493a5« (f ) 

Whenever the quantities in the observations are so large 
that the use of logarithms is desirable, these can best be put 
in the same columns directly over the natural numbers. 

41. It is to be remarked that in solutions like the above 
we assume, from our method of forming the observation and 
normal equations, that the observations on the values of the 
function are alone subject to error, the observations on the 
variable being supposed to be exact or to have errors so small 
as to be negligible. 

42. Feriodio Phenomena. If as the variable increases the 
function passes through recurring values, that is, if y is a 
periodic function of a;, some form of trigonometric equation 
would be the prober one to select. For instance, a good form 
to use is 

y = A-\-Bsin ^x + C cos ^x (J8) 

where Ay JB^ C are the constants whose most probable values 
are to be found, and m is the number of units of x comprised 
in the entire cycle of values of y. The quantity m is to be 
determined from an inspection of the observations, and if it 
appears that several values of m might be used, all should be 
tried, and that which leads to the smallest sum for the squares 
of the residuals is to be considered the best. If the several 
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cycles are not similar and regular, additional terms involving 
multiples of x will have to be added to equation (28). 

43. Special Treatment of Exponential Equations. If the 
equation selected is not in the linear form as regards the 
unknown constants the general method of procedure is 
given in paragraph 44, but in some special cases a more 
simple reduction is possible. A quite common case is the 
following: — 

Suppose the relation between x and y is expressed by the 
equation 

y = ifcx"* (29) 

the problem being to obtain the best values of k and m. 
Taking the logarithms of both members of (29), 

log y = log k-\-m log X 

Denoting log k by k\ this equation becomes 

m log X'-\-k' — log y = 

which is in the linear form as regards the unknowns m and 
k'^ and the normal equations may now be formed in the usual 
manner. 

The most convenient way of determining whether equation 
(29) is a suitable one to select or not is to plot the corre- 
sponding values of x and y on logarithmic cross-section 
paper. If (29) is a proper equation to use the plotted points 

will lie upon a straight line. 

REDUCTION OF OBSERVATION EQUATIONS TO THE 

LINEAR FORM. 

44. When the observation equations are not linear as 
regards the unknowns, the only practicable method of pro- 
cedure, as already mentioned in paragraph 23, is to reduce 
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them to that form. This reduction m&j be e£Eected as 
follows : — 

Let the observation equations be 

(A) 



fn{Zi Zj, . . . Zq) = Jfjj 

in which Zj, Zg, . . . Z^ represent the unknown quantities and 
Mx^ M^^ , , . Mn the observations, the functions being of 
known form. 

Let Zi', Za', . . . Zq be approximate values of Zj, Zg, . . . Z^, 
found by trial or the solution of a sufficient number of the 
observation equations, and let the most probable values of 

^l9 ^2) * * * ^q ^^ 

Z^ + Z,, Z^-{-Z^,...Zq'^Zq (B) 

2i, gj, . . . Zq being small corrections whose values are to be 
determined by the Method of Least Squares. The first 
observation equation in (A) may then be written, 

/l(Zi' + 2i, Z^ + Z^, . . . Z^^Zq) = M, 

Expanding the first member by Taylor's Theorem, denoting 
/i(Zi', Z3', . . . Zq) by Atj, and neglecting terms containing 
powers of 2i, z^ . » . Zq higher than the first, this becomes 



'9' 



*^+ aZi' ^^ + dZ^' 2a+ . . . ^^ / Zq — Ml 



If now we represent the coefficients of Zj, ^a? • • • ^Jg, by 
Oij ^ly • . • ^1, also A^i — Ml by m^, and treat the other 
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equations in (A) in the same manner, our observation equa- 
tions will take the form 

Ol^i -|- h^Z2 + • • • 2\^q + ^'^l = ^1 
«2«1 + ^2«2 + • • • 2%^q +^^2 = ^2 



• • • 



(C) 



«n2Ji + K^2 + • • • 2n^q + ^n = V» 

The second members reduce to the residuals since the 
quantities 2i, 22> • • • ^q represent merely the most probable 
values of the corrections. Equations (C) can now be solved 
in the usual manner and the most probable values of 
Z^, Za, . . . Zq found by substituting the results in (B). 

45. Example. From the following observations find the 
most probable values of x and y. 

sin X -f- cos 2y == 1.5 

cos 05 4" 3 sin y = 1.7 (A) 

aj^ -f 5y = 2.1 

By trial it is found that 42° and 18° are approximate 
values of x and y. Then in accordance with paragraph 44 
we put 

a5 = 42°4-2i y=18°+22 (B) 

and expanding the different functions, we find 

Ai = sin 42° + cos 36°, -g^ = cos 42°, ^^ = - 2 sin 36° 
= 1.48 = .74 = - 1.18 

dk dk 

*a = COS 42° + 3 sin 18°, g^' = -sin42°, ^" = 3 cos 18° 
= 1.67 = - .67 = 2.85 



= 2.11 = 1.47 = 5 



Ai - Jfi = - .02, Aa - Jfa = - .03, A, — Jf, = .01 
= mi, = ma, = m. 
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Therefore we have for 

OBSERVATION EQUATIONS 

.74si - I.I822 - .02 = Vi 
- .6721 + 2.8522 — .03 = V, (C) 

1.4721 -i- 52a + .01 = tJ, 

From these are obtained the 



NORMAL EQUATIONS 



(») 



3.I621+ 4.5722 + .020 = 
4.5721 4- 34.5l2a - .012 = 

Solving, 

2i= -.00845 2a = .00147 (E) 



These results are in circular measure. Reducing to degree 
measure we have 

2i= -.29'.0 2a=5'.l (F) 

Substituting these in (B), 

X = 41^ 31'.0 y = 18° 5M (G) 

46. If the observations are conditioned, precisely the same 
method will be followed in reducing all the equations to the 
linear form. The rest of the solution will then be as usual. 

If the values found for 2i, 22, . . . Zq should turn out to be 
so large that the terms involving their second and higher 
powers can not be neglected as assumed, the process must be 
repeated using the values of Zi, Z2, . . . Zg first obtained as 
approximations. 

In a few cases the equations may be reduced to the linear 
form by some special artifice of a simple character, as in 
paragraph 43. In such cases the method of expansion by 
Taylor's Theorem should not be used. 



CHAPTER m. 

THE PRXSOISION OF OBSERVATIONS. 

47. The work so far considered has treated solely of the 
methods by means of which the most probable values of the 
unknown quantities may be determined from a series of 
observations. But in general something more than this is 
desired. We wish to know, if possible, how much reliance 
can be placed upon the results obtained, and how they com- 
pare in precision with other determinations of the same 
quantities. Preliminary to the discussion of the precision of 
our results it will be necessary to consider more fully the 
Curve of Error. 

48. The Constant jb. The law of distribution of errors of 
observation has been shown to be represented by a curve 
whose equation is 

y = ^6"-^*^" 

In this, if 05 = 0, y = h (30) 

Therefore the constant k represents the intercept of the 
Curve of Error on the axis of Y. It is not, however, an 
independent quantity but is determined by the value of A, 
as will now be shown. 

49. To Find h in Terms of h. Since, as shown in para- 
graph 18, the total area between the curve and the axis 
of -3r is denoted by unity, we have 

k f^e-^^^ dx=l or k f V^"^' ^« = | (a) 
This may be written 

Jo 2A ^ ' 
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Let t = hx, ... dt = hdx. Also, when aj=oo, «=qo 
and when « = 0, « = 0. ' 

•■• r''^''^'=r'-^'^^<^ (c) 

Multiplying this equation by 

•^O Jo 

we nave 



— 1 r* 

~2Jo 1 



da; 



+ a^ 



= i[tan-..];=J 



•• X"''"*'^' = T (81) 

Substituting (31) in (c) and (b) 

^ — A A 

2 ~" 2A "'* * = ^ (82) 

Therefore equations (16) and (17) become 

v'-nr 
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50. It has been shown that the quantity A is a <^ Measure 
of Precision " of the observations, and hence a' determination 
of its value in each case would enable us to compare the 
relative reliability of different measurements. In practice, 
however, it is not found convenient to compute the value of 
h directly, and so this quantity is used only in developing 
the theory of the subject, while in the comparison of observa- 
tions certain other quantities now to be derived are used as 
precision measures. These latter quantities are called the 
"Mean of the Errors," the "Mean Error," and the "Probable 
Error," respectively, and may be computed directly from the 
observations. Further, as it will be found that they all bear 
a definite relation to A, the value of this quantity can be 
determined from them if desired. 

In the following discussions no distinction is made between 
positive and negative errors of the same numerical magnitude, 
and unless otherwise stated the observations are all of the 
same weight. 

THE MEAN OF THE ERRORS OR AVERAGE DEVIATION. 

51. The Average Deviation (a.d.) of an observation is 
the arithmetical mean of the errors all taken with the positive 
sign. 

Since from (6) the probability that the error of a single 
observation will fall between x and x-\-dx hs 

^{x) dxj 

if n observations are made, the number of errors falling 
between these limits is 

n il>(x) dx. 
Hence the sum of all the errors of the observations is 

nix ^(a;) dx or 2ni x ^(sb) dx. 
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Dividing this last expression through by n we have 

.00 2A /»• 

05^(05) dx = -— \ 

Vir*^o 






= - -T= r"^'"^'*' (-2A»a5) efe 



or a.d. = — — (84) 



THE MEAN ERROR. 

52. The Mean Error (fi) of an observation is the square 
root of the arithmetical mean of the squares of the errors. 

The total number of errors being n, the number falling 
between x and x-\- ckc is, as just shown, 

n <l>(x) dx^ 

and the sum of the squares of these errors is 

n x^ <l>(x) dx. 

Therefore the sum of the squares of all the errors is 

n i Q^ ^(^) dx 
h /•• 

But as shown in paragraph 49, 



.-. |i« = A p ^-;i«x» ^ ^ (a) 






(fe) 



Differentiating (b) with respect to A, 



•/ — OB A 
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and replacing the integral in (a) by its value as determined 
from equation (c), we have 



THE PROBABLE ERROR. 

53. The Probable Error (r) of an observation is an error 
such that one-half the errors of the series are greater than it 
and the other half less than it. Or it is an error of such a 
magnitude that the probability of making an error greater 
than it in any given observation is just equal to the probability 
of making one less than it, both probabilities being one-half. 

The probability that the error of an observation will fall 
between x and x-\-dx being <l}(x) dxj the probability 
that the error will fall between the limits r and — r is 



JP = r4>(x)dx = ^ Ce-^^^dx (36) 

9 

If r is the probable error, JP is one-half, or 



2A 



Ce-^'^'dx = t (37) 



rr ^0 



and from this definite integral r is to be founds 

Let t = hx^ .•. dt = h^. Also when aj = r, we have 
t = Ar, and when aj = 0, ^ = 0. Substituting these results 
in (37), we have 
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Denote hr by p. Then by interpolation in a table of 
values of this integral, the value of hr in (38) is found to 
be 

p = hr = .47694 (39) 

^ = *: = -*- <*•> 

54. When t is small the values of | e"^ dt may be 

•/ 

found by expanding e"** into a series and integrating the 
successive terms. Thus, by Maclaurin's Theorem, 

^ ^ [2^ ~ |£ ~ 
.3 ^ 6[2^ 7(3^^ 



RELATIONS BETWEEN |L, r, a.ll., h^ AND p. 

p 1 

55. From (40) ^ = 7, and from (34) a.d. = 



••• r z= p a.d, \v 

or r = .8153 a.d. (41) 

1 
Also from (35) p. = 



r = pp)/2 
or r = .6745 |i (42) 



THE PBBCI8I0N OF OBSEBVATIONS, 43 

The relation between the valnes of /yi, r, a,d.^ and h may 
be conveniently expressed as follows : — 

1 r 
|i^ = •- = -= a.d. 0F (43) 

or arranging in tabular form, 





H- 


r 


a.d. 


/* 


1.0000 


1.4826 


1.2533 


r 


0.6745 


1.0000 


0.8453 


a.d. 


0.7979 


1.1829 


1.0000 



From this it will be seen that 

fi > a.d. > r (44) 

56. Farther, since in (25) it was shown that p oc h\ it 
follows at once from (43) that 

That is, the weights of different determinations of a 
quantity vary inversely as the squares of their Mean Errors^ 
their Probable Errors^ or their Average Deviations. 

It is to be observed, however, that the determination of the 
relative weights of quantities from a comparison of their 
precision measures according to (45) applies only when the 
quantities are of the same kind and subject to the same con- 
stant errors, if any of the latter exist. (See § 98.) The 
applications of (45) are numerous and important. 

57. Mcample A. Suppose n direct observations, all of 
the same weight, be made upon a quantity, and that the 
probable error of a single observation is r. Then since the 
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weight of the arithmetical mean is n, its probable error r^ 
will be given by 

— = -- or ro= — (46) 

Or in general, suppose S is any precision measure of an 
observation of weight jt>, and suppose Po is the weight of a 
second similar quantity or observation, then the correspond- 
ing precision measure So of the latter will be 



= «Ai 



(47) 



The case of most conmoion occurrence is that in which 
/> = 1, and then we have ^ -, 

80 = -^ , ;. -- -- (48) 

JSxample B. A line is measured five times and the aver- 
age deviation of the mean (A. D.) found to be .016 feet. 
How many additional measurements are necessary in order 
that the A. D. of the mean may be reduced to .004 feet ? 

Let X be the total number of observations required. Then 

a; : 5 = .000256 : .000016 

a; = 80 

Consequently the number of additional measurements 
required is 75. 

Eo&ample C, In two determinations of the quantity X 
there were obtained 

Xi = 427.320 ± 0.040, X^ = 427.30 ± 0.16 

Find their relative weights, and the most probable value 
of X and its probable error. 
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N^ote. The above is the method commonly employed to 
denote that the probable errors of the observations are 0.040 
and 0.16. 

From (46) ^ = -— = — 

From (3), the most probable value of X will be given by 

= 427 + '' X ''^ + ■'' 
= 427.319 

and the weight of Zq being 17, by (48), 

To = = .039 

Therefore we should write the result 

Xo = 427.319 ± .039 



REPRESENTATION OF |l, a.a.^ AND r ON THE CURVE OF 

ERROR. 

58. To find the points of inflection of the Curve of Error 
we have 

y = ke"^^^ 

ax 

For a point of inflection, ^— j = 0, [or 

Ih^k e-^"*' (2 A^a-a _ 1) = 

1 
X = = M by (36). 
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That is, the Mean Error is represented by the abscissa of 
the point of inflection of the Curve of Error. See OM in 
figure 2, page 9. 

Next, for the abscissa of the centre of gravity of the area 
to the right of OTJ we have' 

C y X dx 
Xo= OJ) =z ^' 



o 



dx 



2 A >^0D 

~V7 Jo ' 



e-^^x dx 



for, (§ 18), f% ^^=\ 



Integrating, 



Xq = — ~ = a.d. by (84). 

hSfir 



Finally, if an ordinate PP' be drawn so as to bisect the 
area to the right of the origin between the curve and the 
axis of -XJ the Probable Error will be represented by the dis- 
tance of this ordinate from the axis of Y, For the proba- 
bility of the occurrence of an error less than the amount OP 
is then equal to the probability of the occurrence of an error 
greater than OP. This being the case, OP is the probable 
error by definition. 



CHAPTER IV. 

OOBffPUTATION OF THB PRSOISION BCBASURBB. 

DIRECT OBSERVATIONS. 

59. Observations of Equal Weight. Given n direct obser- 
vations all of the same weight on a single quantity M^ to 
find the Mean and Probable Errors and Average Deviation 
of a single observation and of the Arithmetical Mean. 

Let the observations be J/^, J^, . . . J!df^. 

" " arithmetical mean be Mq. 

" " real errors be a^, a^, . . . x^^, 

« " residuals be Vj, v„ . . . v^. 

Denote the mean and probable errors and average devia- 
tion of a single observation by /i^ r, and a,d.y respectively, 
and the corresponding quantities for the arithmetical mean 
hy /i<h ^oj and A.J), Then by definition, 

Vi= Ml -- Moy t?, = Jfi — Mq^ ... Vn = Mn — Mo 



and 



=v^ 



If JU^ represented the true value of M^ the residuals 
would be the same as the real errors, and we should have 



.=V¥ 



and if n is large this formula is practically exact. But when 
n is small a more accurate expression is necessary. To 
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obtain this let Mo + Xq be the true value of M. Therefore 

«i = ^1 — (Mo -\-Xo) = Vi — Oo 
aj, = Ji; — (Mo + JBo) = Va — Oo 

• •• ••• ••• 

• •• ••• ••• 

aSn = Jf; — (Mii + ^) = ^n - «o 

Squaring, adding, and dividing by n 

^' = /*' = ^ (^^' - 2 asaSv + najo") 

= ^ -{- xo' for by (2) Si; = 

The value of a^ is not known and can not be found 
exactly, but it is approximately equal to the mean error of 

Moi that is, by (48), to fio = -^. 

Substituting this value in the above we have 

^ n ~^ n 

(n - 1),A» = Sw* (49) 



" n— 1 



(«•) 



by (48) m> = SJ ^Z"' .s (") 



n(n — I) 



/ 

by (42) r = .6746 ^ ^ri - («> 



and n = .6746 \ ^,^\^ (68) 



«(n.— I) 
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60. In order to avoid the use of the squares of the 
residuals we may proceed as follows : From (49) 

S»« = (n - 1) ^ 
On the average, the values of the residuals will then be 



//I— 1 



• • • 



In - 1 



Adding and dividing by n, neglecting the signs of the 
residuals, 

%v / w — 1 SaJ / n — 1 

n \ n /I \ n 

a.d. = — (54) 

>Jn {n — 1) 

by (48) A.D. = — (55) 

n^n — 1 

by (41) r = ^ (5«) 

Vt^(t^— 1) 

J, t8453 2!f^ .^M. 

and fo = -z^zziz (W) 

n Vn — 1 

The mean errors may also be computed from the above by 
using the table of equivalents in paragraph 55, but this is 
not customary, formulas (50) and (51) being used for this 
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purpose. Results derived from (50) are to be regarded as 
more accurate than those obtained from (54), the latter 
being a second approximation. 

61. ^Example. From the following measurements on the 
length of a base line find the most probable length and the 
values of the various precision measures. 



M 


V 


«* 


455.35 


.02 


.0004 


.85 


.02 


4 


.20 


-.18 


169 


.05 


-.28 


784 


.75 


.42 


1764 


.40 


.07 


49 


.10 


-.23 


529 


.30 


-.03 


9 


.50 


.17 


289 


.30 


.03 


9 


3.80 


+ .70 


.3610 


455.330 


-.70 




Jf, 


S« — 


S»* 



By (50) and (51) 



= v 



.3610 
9 



r = .6745 /I 
By (54) and (55) 



= .20 



= .13 



/A 



= .063 



VlO 
ro = .6745 li^ = .042 



a.d, = 



1.40 
VT&O" 



= .15 



A.D. = 



1.40 
"30" 



= .047 



r = .8453 a.d. = .13 



Ta = .8453 A.D. = .040 
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and we should write for the most probable length of the 
base line 

Jfo = 466.330 ± .042 

62. Observations of Unequal Weight. Using the same 
notation as above, with slight modifications, we will 

Let Mq represent the General Mean. 
" Ply i^2> • • • Pn ^6 the weights. 
<< a.c?.i, a.d.^y • . • a.d»n be the average deviations. 
^ fii, iHi • • • I'-n ^^ ^6 mean errors. 
" **!> ''2? • • • ''n he the probable errors. 
** a.d.^ fji^ r, and v refer to observations of weight unity. 

Then by (48) 

If the ^^ Observation Equations " are formed for this case 
they will be 

J^ — JMi = t>i, M2 — Mq = Vfy . . . Mn — Jli = t?». 

And, as was shown in paragraph 28, if these equations are 
each multiplied by the square root of the weight of the corre- 
sponding observation, they will all be reduced to equivalent 
equations of weight unity. On performing this operation 
it will be seen that the residuals of the new equations 
become 

t>i >rpiy v^ yTpt, ... v^ ^n 

And evidently to these reduced observations the formulas 
of paragraph 59 apply. 



1 
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Therefore 





j: 


Sjpv* 




n-l 




"Zpv* 




p^{n - 1) 




Siw* 



r = MU |i (58) 



rj = .«744 |ij (5f ) 



SjXn - I) '• = •«« "^ (••> 



Alcio, by a method similar to that used in paragraph 60, it 
may be shown that 



St^V^ 



»*c^* 




>Jn(n- 1) 


n.m.n...t^ 


— 


SvVi> 


\M/%vV*f^ 


Vl>jk w(n - 1) 


A.T>. 




svVjp 



r = MSia.d. (61) 



fj^ = ,8I5S a.d.ic (62) 



^^pn{n — I) 



n = .8153 ^.1>. (63) 



Formulas (61), (62) and (63) are not of much value in 
practice unless all the weights are perfect squares, for other- 
wise no labor is saved in the computations, since the square 
root of each weight will have to be determined, and the results 
obtained can not be considered as reliable as those given by 
formulas (58), (59) and (60). 
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63. Example, Given a series of observations on M, the 
difference in longitude between two stations. To find M^ 



M 


P 


pM 


V 


v^p 


«» 


pv* 


V 4' 30" 


4 


120 


1.6 


3.2 


2.6 


10.4 


41 


1 


41 


9.4 


9.4 


88.4 


88.4 


43 


1 


43 


11.4 


11.4 


180.0 


130.0 


37 


9 


888 


5.4 


16.2 


29.1 


261.9 


48 


4 


192 


16.4 


32.8 


269.0 


1076.0 


34 


16 


544 


2.4 


9.6 


5.7 


91.2 


25 


9 


225 


6.6 


19,8 


43.6 


392.4 


46 


1 


46 


14.4 


14.4 


207.4 


207.4 


28 


25 


700 


8.6 


18.0 


13.0 


325.0 


24 


4 


96 


7.6 


15.2 


57.7 


230.8 


1» 4' 81".6 


74 


2840 




150.0 




2813.5 


-Mi 


Sp 


Sjtwn 




%v<i/p 




Sj»t>* 



By (8) 



S pM 2340" 

JUi = =f— = 1» 4' -f 



Si> "74 

= 1»4'31".6 



= v/ 



2813.5 



9 



= 17.7 






\/u 



= 2.1 



r = .6745 a = 11.9 



150 
a,d, = = 15.8 

^90 



r 
3 



= - = 4.0 



. _ a.d. ^ _ 
-4.i>. = -zz:= 1.8 

V74 



ro = .6745 ^ = 1.4 = .8453 A.D. = 1.5 



Mo = r4'31". 6 ± 1".4 
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FUNCTIONS OF OBSERVED QUANTITIES. 

64. Theorem. Given any number of quantities and their 
Mean and Probable Errors and Average Deviations, to find 
the Mean and Probable Errors and Average Deviation of any 
function of the quantities. 

Let the quantities be Jt/^, J/^, . . . Mq. 
^' mean errors be fA^ fi^y > • . fiq. 
" function be Jf=/(J/i,JI^ . .. -M^)- 

" mean error of Jf be JE'. 
" probable error of M be H. 
" average deviation of Jf be j9. 

The derivation of the general formula will be simplified if 
we consider first a few special forms of functions. 

65. Case Z Suppose JK" = Jtti ± Jfj. 

The number of observations from which Mi and M^ and 

hence /ai and /a^j have been determined is not necessarily 

known, but we may assume that for each quantity it is any 

large number n, and that the real errors of the observations 

are 

for Jfi, xi\xi",xi"\ 



/// 



• • 



Then the real errors of My computed from the separate 
observations on Mi and M^ will be 



HI 



Xi'±X^\ Xi"±X^", Xi"'±X^' 



M^ = 



n 

SsBi*±2Sa5iaJ2 + Sa52* 



n 
or H^ = ^i^ + m^ (j4) 

since in the most probable case the term Sa^fl^ will dis- 
appear, as there most likely will be as many positive as 
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negative products of the same absolute magnitude of the 
form i^i X2> 

By successive applications of the above the same principle 
may be extended to cover the algebraic sum of any number 
of quantities. So that if 

then i?« = m« + m« + . . . |L^a = S|i« («5) 

Since probable errors and average deviations differ from 
mean errors merely by a constant factor, we shall likewise 
have 

U« = rx' -{- r,^ -{- . . . Tq^ = Sr» (06) 

1>» = a.d.x^ + a.d.^^ + . . . = S a.d.^ (67) 

66. Example. Given the telegraphic longitude results, 

h, HI. seCm sec* 

(A) Cambridge west of Greenwich, 4 44 80.99 ± 0.23 

(B) Omaha « " Cambridge, 1 39 15.04 ± 0.06 

(C) Springfield east " Omaha, 25 8.69 ± 0.11 

Find X, the longitude of Springfield, and its probable error. 

Z = -4 + ^ — (7 = 5* 58"* 37'.34 ± 0'.26 
for by (66) B == ^^(.23* + .06« + .IV) = .26 

67. Case IL Suppose JH = aiMi 

Using the same notation as in paragraph 64, the real 
errors of M will be 

«i«i'> <hP'h!'t <hPh'"j • • • 

2&* = = Oi Ml 






or JE = aiVn. (68) 
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68. Case III. Suppose M =^ aiMx ± a^M% > . . UqMq 

By combining (65) and (68) 

also JJ^ = SaV^ (») 

and n^ = Sa^a.d.' 

69. ^Example, The length of a bar at 20^ Centigrade is 
found to be 75".0041 ± 0".0037, and its expansion per 
centigrade degree 0".0036 ± 0".0018. What is the length 
of the bar at 56*" Fahrenheit ? 

20^ C. = - X 20 + 32 = 68** P. 

5 ^ 

5 5 

The expansion for 1** F. will be - X .0036 ± - X .0018 

or .0020 ± .0010 



•*• 



L = 75.0041 - 12 X .0020 

= 74.9801 
B = V.0037«+ (12 X .001)« 

= .013 
Z = 74''.980 ± 0".013 

70. Case IV. Suppose Jf = f (JIfi, JMTj, . . . J£^). 

Let J^ = Oi -f- ^> -3^ = 02-)- ^a, . . . Mq = a^ -)- m^ 
where Oi, o^, • • . a^ are arbitrarily assumed quantities very 
nearly equal to M^^ it^, . • • Mg^ so that m^, m^ . , . m^ are 
so small that their second and higher powers may be neg- 
lected. Then the errors of -3^, -3^, . . . Mq may be con- 
sidered to be in m^, m^ . . . m^, and hence /Hj fi^ • • • fiq 
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may be regarded as the mean errors of the quantities 
m„ mg, . . . niq. We now have 

M = / (Oi -f ^i» «2 + ^a> • • • ^g + ^g) 

Expanding this expression by Taylor's Theorem, and 
denoting /((hj a2? • • • ^g) ^7 ^'9 we have 

M = M' 

dM' . dM' . dJf' 

+ ^^ -dST + ^» -d^ + • • • ^^ d^ (^) 

-j- negligible terms in the second and 
higher powers of mi, m^j . . . niq. 

Then the mean error of Jf will be the same as the mean 
error of terms (A), and by (69) this is given by 

or what is practically the same thing, by 

- - = -• [li.T +••••-• [:-iJ <'•> 
^=<-«'[li.]'+---.{lfJ' 

71. Example A. Two sides of a right triangle are meas- 
ured with results 

a — 49.53 ± 0.69 
h = 60.38 ± 0.93 
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Find the length of the hypotenuse c and its probable 
error. In this example 

M= ^Mi'+JiQ = V«* +b^ = 70.65 

dM a dM b 

ailfj ^ a" + *» ' 9-^; v' a* + 6* 

By (70) i2«= V^|-5. + r,«^^, 

|_70.65j ^ |_70.66j 






R = .78 
and c = 70.65 ± 0.78 

Example JS. If the probable error of x is r, what is 
the probable error of the common logarithm of a; ? In this 
case 

^ = logio ^1 = logio « 
dM d logio 85 logio e 





dJi^ da; a; 


••• 


n- ^^"% 



a; 



Example C, If the weight of a; is jo, what is the weight 
Po of sin a? ? 

Denoting the mean errors of x and sin a; by /i and E^ 
respectively, we have, by (45) and (70), 



_^ . r^ sinaj"]^ 

^ = ;.« -^— = ,i« COS« 



X 



and 
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^ _ jU^ 
p E^ 



P 

Pq = — -— = p sec* X 

cos* a? 

72. Equation (70), which expresses the law of propaga- 
tion of error in functions of observed quantities, is one of the 
most important in the whole theory of the Method of Least 
Squares. Upon it in particular is based the discussion of the 
"Precision of Measurements." This subject treats, in the 
first place, of the methods of finding the precision of a 
quantity obtained by computation from a series of measured 
quantities ; and, in the second place, it investigates the pre- 
cision with which the component measurements of a series 
must be made in order to obtain a required degree of pre- 
cision in the final result. The following simple example will 
illustrate the character of the solutions : — 

73. Example. In the determination of a current by a 
tangent galvanometer we have 

I = 10 — tan <^ 
Gf 

where I is the current in amperes, ST the horizontal com- 
ponent of the earth's magnetic force, G the galvanometer 
constant, and ^ the angle of deflection. Given the errors 
^9 Sg, $89 in i9^ (? and tan ^ to find the error A in Z 

By (70) 

-v[||J+v[|i]-+v[^T 



100 «, . . 100 -^«. , . 100- „. 
= — V tan» .^ + -^j- V tan« ^ + -^P- V (a) 
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Dividing this equation by I^ = tan* ^ we have 

m"=[ij+[i>[wj <-> 

That is, the square of the percentage error in Z is equal 
to the sum of the squares of the percentage errors in Sy 6j 
and tan ^. Hence if 

ff is determined within .4 per cent, 
tan <^ " " " .1 " " 



then — = V -16 -f .04 + .01 = .46 per cent, (c) 

Next, suppose the value of T is required to within .1 per 
cent. To find the necessary accuracy in the determinations 
of J9^ &, and tan ^ when the error in each of these quan- 
tities is to have the same influence upon the total error. 

From (b) we shall now have 

.001. = s [ij = 3 [A]- = , [^J <a, 



• • 


8, .001 „„„,„„ 
-± .000577 


• 


8, — .00068^ 




8, = .00058 G 




St — .00058 tan ^ 



(e) 



It is comparatively easy to obtain the necessary accuracy 
in the measurements of Q and tan <^, but difficult in the 
case of S 

For additional work of this kind see Holman's '< Precision 
of Measurements." 
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74. Combination of Functions of the Same Variables. 
It is to be noticed that equation (70) applies only when M 
is a function of independent quantities. If Jt/^, M^^ . . . Mq 
are merely different functions of the same quantities we 
must proceed as follows : — 

Let JIfx = ^ (2^1, ^29 • • • ^k) 

J^ = ^ (2l, Zg, . . . z^) 

M = f ( JIfi, Jf,) 

If any single observations of Zi, z^^ , , , zjg are subject to 
errors Xiy x^^ . . , X]tj the corresponding errors in Mi and 
M^ will be 

for Miy Xi = OiflJi -}- «8«2 + • • • ^*«* (a) 

" JMi, ^ = Oi'aJi -|- <*a'«2 + • • • «* «* (b) 

Where a^, a,, . . . aj^ are the differential coefficients of 
3fi^ and ai\ a,' . . . a^.' the differential coefficients of J[^ 
with respect to 2i, sja* • • • 25*. The corresponding error in 
M will then be 

X = AXi + JL'Xj (c) 

Where A and ^' are the differential coefficients of M 
with respect to Mi and M^. Substituting in (c) from (a) 
and (b)y 

X = (Aoi + ^Vi) aji + (Aa^ + JLV,) aja . . . 
= o.Xi -j- ^aja + . . . XaJt 

Then if the number of observations or values of X be 
denoted by n, 

SX* San* %xJ 

E^ — aVi^ + iS^Mi* + • • • ^V*^ (d) 
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since in the most probable case the product terms will 
cancel out. 

Expanding (d) we have 

JS^ = (Aa^ + A'oi')^ /iii« + ( JLoa + A'a^'y ^« + . . . 

= A\a,%^ + OaV . . .) + ^'\a^'W + <hW . . .) 

75. Example. As a very simple problem take 

Jfi = 22i, M^ = 32i, /Ai = 0.1, 

and Jf = J/i -|- J8^. 

Then -4 = 1, -4' = 1, Oi = 2, a/ = 3. 
By (71) ^« = 4 X .01 + 9 X .01 + 2 X 2 X 3 X .01 = .26 
or E = 0.5 

In this particular example the result may be found directly 
from (68) by substituting at first in M the values of Mi 
and M^. 

Thus Jf = 22i -f 325i = 52i 

-^ = 5/i.i = 0.5 

If Ml and M^ had been independent quantities, by (64) 
or (69) we should have had 

E = V(2 X 0.1)2+ (3 >< 01)2 
= 0.36 

INDIRECT OBSERVATIONS. 

76. The determination of the precision measures of the 
unknown quantities in case the observations are indirect 
involves a knowledge of the weights of the unknowns, and 
consequently the method of computing these weights must 
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first be demonstrated. It will be assumed at present that all 
the observations are of weight nnity. 



FIRST METHOD OF COMPUTING THE WEIGHTS. 



Let the observations be 









in which Jf^, Jf^, . . . Jf^ denote the actual observations, 
and ^1, ^29 • • • ^9 ^^ most probable values of the unknown 
quantities. Let 

A?i — J/i = mi, Ajj — J^ = mj, ... k^ — M^ = f^n 

Then the above equations give rise to the 



OBSERVATION EQUATIONS 

«i2i + *i2a + • • • 3'i2« + ^1 = 



«a«i + *a«a + • 



<hfi\ + Ma + • 



3'a«g + ^a = 






^n2« + ^n = Vf 



By the rule, paragraph 25, we now form the 

NORMAL EQUATIONS 

2i S a* + ^2 2 ^^ + • • • ^q^H + S *m 



(A) 








(B) 



... 



... 



2i S afi' 4" ^a S ^2' + . . . 2g S 2'" +2 fi'^ = 
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Multiply the first of (B) by Qiy the second of (B) by 
Qt9 ' ' ' and add the results. Then 



2i(CiSa* + Qtiab -\- . 


. . CgS aq) 


• • • • • • 


• • • 


• • • • • • 


••C,S?») 


+ CiS am + CaS ftm + . . 


. Qg%qm 


— 





(C) 



I^t Qi, Qs, ... Qg be determined so that the coefficient 
of Zi in (C) shall be unity, and the coefficients of 2^, 
2s, . . . Zg each equal to zero. That is, let 

CiSa' + Q^ab + . . . Qg^aq -1 = 
Qaab + C»S^' + . . . CgS^y = 

... ... ... I JL^ 1 

... ... ... 

Qi^aq + ftS^y + • • • Qq^q* = <> 

Then (C) becomes 

^ + CiSam + Qi^bm + . . . C«S?m = (E) 

Equations (D) may be derived from the normal equations 
(B), if in them we replace 2i, z^^ ... Zg by §i, ^2? • • • ©g? 
S am by — 1, and S ^wi, S cm, . . . S qm, by zero. Hence 
the solution of the normal equations with these changes will 
give the values of §j, Ca> • • • Qg* 

TO SHOW THAT Qi IS THE RECIPROCAL OF THE WEIGHT OF Zi. 
Expanding the coefficients of (E) we have 
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«i + Qi{<h^i + <^t^ + • • • «n»^) 
+ ft(^i»^i + ^8^3 + • • • ^n^n) 



= 



^nT^n) 



or collecting the coefficients of mi, mg, . . . m^ 



• • • 



• • • 



= 

For convenience in writing we will let 



(E') 



ai 






o» = QiOn + Cs^n + • 



• QqSl 



Cg^i 



Zi -|- aiWli -|- 0)^2 -|- . . . Onplf^ z=z 



(F) 



(G) 



Multiply the first of (F) by Oxj the second by a,, ... 
and add the results. Then 

Soa = CiSa* + Q^^ab + . . . CgSa^^ 

= 1 by first of (D) (H) 

Multiply the first of (F) by b^ the second by ^t> • • • 
and add the results. Then 
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%lHi = Qi%ab + CaS^ + . . . Qq^bq 

= by second of (D) 

Likewise Sca = (H') 



• • • 



%ga = 

Multiply the first of (F) by ai, the second by a,, . » . 
and add the results. This gives 

= Q, by (H) and (H') (I) 

Let fA be the mean error of an observation of weight unity. 
Let fis^ be the mean error of «i. 
Let Ps^ be the weight of 2?!. 

The mean errors of m^, m^ . . . m^ are the same as the 
mean errors of J/^, J^, . . . Jf^i and each is accordingly 
equal to /a. Therefore from (G), by (69) 

ihit^ = oLiW + f^W + • • • a»V* 

= Ci/^« by (I) . .. (J) 

2 •• 

But by (48) ^' = ^ .,, v.V 

Comparing this with (J) we see at once that 

«^ = K (^> 

Therefore, for the First Method of computing the weights 
we have the following : — 
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77. Enle I. In the normal equation for Zi write — 1 
for the absolute term ^am^ and in the other equations 
zero for each of the absolute terms S bm,^ S cm, . . . S qm. 
The value of z^ found from these equations^ is the recipro- 
cal of the weight of the value of z^ obtained by the solution 
of the normal equations. 

To find the weights of gj, «8> • • • ^g> proceed in a similar 
way^ forming a corresponding set of equations for each 
unknown. 

SECOND METHOD OF COMPUTING THE WEIGHTS. 

78. Write equations (B) of paragraph 76 in the form 

«i S «* + 22 S a^ + ... Zq^aq -\- S am = A 

Zx%ab '\- z^^b^ -{-,.. Zq%bq -^^ %bm = B 

... ... ••• f dLj I 

... ••• •.• 

Zx%aq -\- z^'^bq -{-... Zq'^q^ ^ "Zqm = Q 

Then in the solution by the preceding method, equation 
(E) becomes 

«i + CiS aw + gjS bm -^ . , . ^^S qm 

= Q^A J^ Q^B -\- ... QqQ (M) 

in which, as was proved in (K), gi is the reciprocal of the 
weight of Zi. Whatever method of elimination is employed 
in the solution of the normal equations, the coefficient of A 
in the value of Zx must necessarily be always the same. 
Hence we have 

79. Enle II. Write Ay B^ . . . Q instead of zero in the 
second members of the normal equations and carry out 
their solution in any convenient way. Then the most 
probable values of Si, z^^ . . . z are given by those terms 
in the results which are independent of A^ B. . . . Q, 
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The weight of z^ is the reeiprocai of the coefieierU of A 
in the voUue of z^. The weight of z, t« the reeiprocai of 
the coefficient of B in the vaiue of z^, etc^ etc, 

THIRD METHOD OF COMPUTING THE WEIGHTS. 

80. From the second, third, . • • equations of (L) find 
the values of z^^z^ , , , z^ in terms of Zi and snbstitate 
in the first of (L) without redaction. Then the first of 
(L) becomes 

Bzi i=z T + A + terms in JB, C, . . . Q 

Where T is the sum of all the numerical quantities result- 
ing from the substitutions. Dividing through by Hy 

«i = -^ + -^ + terms m JB,C, ... Q (N) 

T 

in which -— is the most probable value of z^ and, as was 

shown in deriving the second method, 

- = e. 
B ^' 

B = p^ (O) 

From this follows at once 

81. Bule III. Substitute in the normal equation for z^ 
the values of z^^z^^ . . . Zq in terms of z^ as found from 
the remaining equations. Then before freeing of fractions 
or introducing any reduction factor^ the coefficient of z^ 
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in this equation is the weight of the value of Zi obtained 
in the solution. 

To find the weights of »2, »8> • • • ^qy proceed in a similar 
way with the normal equations for each of these unknowns. 

For the solntion of an example by the three different 
methods see paragraph 84. 



THE MEAN ERROR OF AN OBSERVATION. 

82. The next step will be to derive fiy the mean error 
of an observation of weight nnity. In the following demon- 
stration the equations referred to by letters are those in 
paragraphs 76 to 81. 

Let the real values of ZiyZ^,..Zq be 

«i + «i» 23 + «tj • • • 2J^ + «ff 
and substituting in (A) we have 

ai(«i + «i) + ^i(2» + aj,) + . . . qi(Zq + Xq) + mi = Ai 
«i(«i -F Oi) + *«(«« + «») + ••• M^q + «g) + ^j = A, 

• •• ••• ... IJT I 

• •• ••• ••• 

«*(«! + «i) + ^n(«a + aJ,) + . . . qn(Zg -f a;^) + m„ = A» 

where A^, A,, . . • A,» are the real errors of Mi^ M^^ . . . M^^ 
or of; m,xj m29 . . . m^. 

Multiply the first of (P) by Oi, the second by o^, . . . 
and add the results. This gives 

2i S a* + z^^ab -\- ... Zq%aq '\' S am 
+ aji S a* + a^ S «^ + . . . aj^ S ag' = S a A 

But by the first of (B) the first line in this equation is 
equal to zero, and therefore 
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Xi'S a^ -\' X2% ab -{- . . . Xq% aq — % aA = 
Also ajj S a^ + iCa S ^* + . . . aJg S ^g- — S ^ A = 

• •• ••• ••• I V^ 1 



• • • 



Xi^aq -{- Xi^bq -\- . . . Xg^q^ — S g^A = 

These being of the same form as the normal equations 
(B), the value of jBi resulting from their solution will be 
of the same form as that of Zi in the solution of those equa- 
tions, with only the substitution of — A for m. From 
(G) we shall therefore have 

Xi — ctiAj — OfAs — ... OnAfi = (R) 

Multiply the first of (P) by Vi, the second of (P) by 
t^s, . . . and add the results. Then 

(Zi + Q^)'Siav + (2j + x^)%bv + . . . (2^ + Xg)%qv 

-\- Sww = SAt) 

and multiplying the first of (A) by a^, the second by 
Os, . . . and adding the results 

S at; = 2i S a* + 22 S a^ + . . . 2g S fits' + S ^^ 
= by the first of (B). 

Also ^bv = 

• . • 
-Siqv = 

Substituting these values in the above, we find that 

Sww = SAt) (S) 

Now multiply the first of (A) by Vi, the second by 
Va, . . . and add the results. This shows that 
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«! S «v + 2, S ^v + ... 2^ S 2'v + S mv = S V* 
and as above, 

S av == Xbv . , . = %qv = 
Combining this result with (S) we have 

Smv = Sv* = SAv (T) 

Next multiply the first of (A) by Ai, the second by 
A}, . . . , add the results and compare with (T). This gives 

gJi S aA -f- ^2 S ^A -f- • ' • 25^ S g'A + S mA 

= SAv = Sw* (U) 

And finally, multiplying the first of (P) by Ai, the second 
by As, ... and adding the results, we have 

Zil^aA -{- Z2%bA -\- ... Zg%q^ -\- S mA . 
+ «! S aA + aJa S ^A + • • • ^^S 3'A = S A^ 

Therefore from (U) 

Sw' + aJiSaA + »2S^A + . . . Xq^qA = SA* (V) 

» _ SA^ _ ^ a^SaA + . . . g^SgA 

We must now find the mean values of the terms Xi 2 aA, 
x^ 2 ^A, ... Xq'S qA. Expanding 2 aA, 

S C(A = ^Ai -|- OsAj -|- . . • ^A,^ 

from (R) Xi = aiAi + osA, + . . . o^An 

Multiplying, 

Xi 2 aA = aiat Ai^ -{- a^a^A^ -|- . . . OmomAh' 
-|- terms in A^A^ A^As, . . • 
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In the most probable case these product terms vanish, and 
snbstitnting for Ai*, Aj*, . . . A„^ the mean value /a^, we 
have 

Xi%a^ = /a' 2 oa 

= ,.* by(H) 

Similarly aijSdA := ^* 



• • • 



From (W) then, q being the number of unknowns, 



n 


-^ n 


\ 


s/-, 


n— q 






\l ^' 


V' 




^ Pz(n 


-Q) 


-r .1 


WiK 4/ 


S«» 



(7J) 






r, = .6745 |i, = ,6745 V / ^ (^^) 

83. By a method similar to that used in paragraph 60, we 
may denve 

St? 
^•d* = — (75) 

si n{n — q) 

a.d.z = — z==zzz==r (76) 

^ Pzfii^n - q) 

84. JExample. In illustration of the above processes 
take the example in paragraph 24, where we found for 



\ 
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OBSERVATION EQUATIONS 



(A) 



and for 



«i — ^2 — 1.7 = ^^1 

2, — 2.4 = w, 

«i + 2a + 28 — 1-0 = t?8 

22 — 28 — 3.0 = V^ 

NORMAL EQUATIONS 

2 2i — 2 2j — 28 — 0.7 = 

2 2i + 3 22 — 2.3 = 

2j -|- 3 gj ^ 0.4 = 



SOLUTION FOR THE WEIGHTS BY THE FIRST METHOD. 

For finding the weight of z^ the above normal equations 
would be written 






(a) 





(b) 





(c) 



22i — 2z,— z,- 


-1 — 


(a') 


— 22i + S a. 


— 


(b') 


- Si + 3a, 


— 


(«') 



Solve for 2i 

3 X (a') 6 2i — 6 22 — 3 28 — 3 = 

(C) - 2i + 3 28 =0 

II III, 

5 2i — 623 —3 = 

2 X (V) -^ 4 2, + 6 22 = 

2i -3 = 

1 

2i = 3 .-. p^^ = y ^^/^ 

For finding the weight of 22 the equations are 

22i-2 2a- 2, =0 (a") 

-221 + 3 2, -1 = (V) 

- 2i +3 2, =0 (c") 



/ 
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Solve for z^ 

3 X (a") 6 «i - 6 Sa - 3 2, = 

(c") - gi + 3 gs = 

5 gi — 6 gj =0 
6 

D 

5 3 

Substitute in (b") g, = - .-. p^ = - (d") 

«> 

For finding the weight of g, the equations are 

2 gj - 2 ga — g4 =0 (a'") 

- 2 gi + 3 g, =0 (V') 
_ zi +3g, -1 = (c'") 

Solve for g, 

4 
from (V') 2 g, = - gi 

3 

substitute in (a'") 2 gi — 3 gj =0 

2 X (C") -2g^ + 6g8-2 = 

3 gg - 2 = 
g, = I ... p^ = J (d-) 

SOLUTION BY THE SECOND METHOD. 

The normal equations will now be modified so as to appear 
in the following form : — 

2gi - 2g3 - g, - 0.7 = ^ (a) 

- 2 gi + 3 gj - 2.3 = ^ (b) 

- gi + 3 gg - 0.4 = C7 (c) 
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Solving, 

3 X (a) 6 2i — 6 2j — 3 », — 2.1 = 3 ^ 

(c) - gi + 3 gg - 0.4 = C 

52j_6ga — 2.5 = 3^ + C7 

- 4 gi + g ga - 4.6 = 2.B 

gi — 7.1 = 3^ + 2^ + (7 (8) 



gi = 7.1 and /?^ = ~ (d) 



Substituting (8) in (b) 



3 g, = 16.5 + 6^ + 5^ + 2(7 



3 

g, = 5.5 and ^^ = ~ (e) 

o 



Substituting (8) in (c) 






3g, = 7.5 + 3^+2^ + 2(7 



3 

g, = 2.5 and p^ = - (f) 



SOLUTION BY THE THIRD METHOD. 
The normal equations are now taken in their original f omu 

2 gi ~ 2 g, — 2, — 0.7 = (a) 

- 2 gi + 3 ga — 2.3 = (b) 

- gi + 3 g, — 0.4 = (c) 
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To obtain «i and its weight we proceed as follows : 

Zi A 

from (c) 2, = — + ~ 

2 2i , 2.3 
from (b) «,= — + — 

Substitute in (a) 

4 4.6 «i -4 A ,7 A 

' 3 ' 3 3 3 

Collecting terms, — 5~ ^^ ^ 

Zi = 7.1 and p^ = - (d) 



. « 



For Z2 
3 X (a) + (c) 5 2i - 6 22 - 2.5 = 

O 

12 
Substitute in (b) _ — z, - 1.0 + 3 a, - 2.3 = 



Collecting terms, -22 — 3.3 = 

5 

2, = 5.5 and p^ = - (e) 

For 28 
3 X (a) + 2 X (b) 2 2i - 3 2, - 6.7 = 

3 ,6.7 

"^ = 2"'+ 2" 

3 6 7 

Substitute in (c) — - 2, - -^ + 3 23 - 0.4 = 
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^ « . 3 7.5 
Collecting terms, - z^ = 

3 

2, == 2.5 and p^ = - (f) 

It is evident that the three methods give identically the 
same results and that the work is about the same in each case. 
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Substituting the values found for s^, z^ and z^ in the 
observations equations (A), we have 

2 



3 



7.1 - 5.5 - 1.7 = Wi = - .1 .01 = Vi 

2.5 - 2.4 = t?2 = + .1 .01 = v^' 

- 7.1 + 5.5 + 2.5 - 1.0 = t?3 = - .1 .01 = w«* 

5.5 — 2.5 — 3.0 = v^ = .0 .00 = v^^ 



In this example /^ = 4, ^ = 3. 



.03 = %v' 



By (72) 



= V^T = ''' 



By (74) r = .6745/* = .12 

By (73) A^ = y I^ = .30 r^ = .20 

^ = y/:^ = .22 r^ = .15 

;*^ = l/:^ = .14 r^ = .10 

f 

3 
By (75) «.<?. = — = .15 

r = .8453 a.d. = .13 
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85. Observations of XTnequal Weights. If the observations 
are not all of the same weight the formulas and operations 
are merely modified in the usual manner and equations (72) 
and (75) take the more general form 



y n — 



(H) 



a.d. = , ^^^- (78) 

V n(n — q) 



CONDITIONED OBSERVATIONS. 

86. Suppose there are given n observations, n' condi- 
tions and q unknown quantities. Then by paragraph 33, 
the method of solution is to eliminate n' unknowns between 
the "Observation" and "Condition" equations, leaving 
q — n' independent unknowns in the " Normal " equations. 
Consequently formula (77) now applies to this case and it 
would be written 



= v/.r^#^ (») 



n — q -{- n' 



also |i« = V ^ .^ ^^:: . ^., (80) 



Pz(n - g + n') 



The weights of the q — n' unknown quantities can be 
found by any one of the three methods already given and 
then the mean error of each unknown may be computed by 
using formula (SO). If the mean errors of the n' quantities 
that were first eliminated is wanted, their weights must be 
determined by eliminating a different set of n' quantities 
from the original observation equations and solving the neces- 
sary sets of equations for these weights. The first method of 
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solntion for the weights would perhaps be best here as the 
actual values of the unknowns have already been found. 

87. Example. Given the 

OBSERVATION EQUATIONS 



(A) 



Zi ^ z^ — 8.0 = weight 1 

z^ _ 2^ ^ 1 5 _ Q "4 

2a — 2.2 = « 3 

Zj^ + 2g — 3.4 = « 2 

and the 

CONDITION EQUATION 

23 _ 2, — 0.5 == (B) 

To find the most probable values of 2^, z^ and z^ and 
also their mean errors. 

Eliminating z^ between (A) and (B) there remains 

2i -j- «a — 3.0 = weight 1 

2i _ 2, -|- 1.0 = "4 

2a — 2.2 = " 3 

2j _|- 2jj _ 2.9 = "2 

From these we have the 



NORMAL EQUATIONS 

7 25i - 2a - 4.8 = 
— 2i + 10 22 — 19.4 = 



(C) 



Solving, 



z^ = 0.98 p^^ = 6.9 

2a = 2.04 /)^ = 9.9 (D) 

from (B) 23 = 2.64 
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Now eliminating z^ between (A) and (B), we find 

«i + 2^8 — 3.5 = weight 1 

2i - 23 4- 1.4 = « 4 

gg _ 2.7 = "3 

25j -|- 23 -- 3.4 = « 2 

From these we derive the new set of 



(E) 



(P) 



NORMAL EQUATIONS 

7 2i — 28 — 4.3 = 
— «i -|- 10 28 — 24.4 = 

and solving for z^ 

28 = 2.54 p^ = 9.9 (G) 

Substituting the values of 2^, 2^ and 2, in equations (A), 
we have the residuals 



V 


v' 




pv^ 


.02 


.0004 




.0004 


.06 


.0036 




.0144 


.16 


.0256 




.0768 


.12 


.0144 




.0288 

.1204 r- %pv^ 


Since in this example 


n . 


4, 


n' = 1, ^ = 3, 



By (79) , = S/:^ = .25 






/*«! = 



/*»! = 



f*a. = 



/* 



Va:9 



= .09 



= .08 



= .08 
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The first significant figure in the mean errors is so large 
that it is not worth while to retain the second place as usual. 

If desired, the probable errors and average deviations can 
now be computed by the usual formulas. 

88. In case the observations are made directly upon the 
values of several quantities subject to certain conditions, we 
have 71 = g', and equation (79) reduces to 



= v/^' («) 



from which the mean error of any observation may at once 
be computed from its weight. 

89. Example. Taking the example in paragraph 34 on 
the measurement of the angles of a quadrilateral, we had for 



OBSERVATION EQUATIONS 



2l = 





weight 3 


Zj = 





« 2 


SSs = 





« 2 


24 _ 





« 1 



(A) 



for the 



CONDITION EQUATION 

2l + 22 + 2^3 + 2, + 58 = (B) 



and for the 



NORMAL EQUATIONS 




4 2i + 22 + 28 + 58 


— 


2i + 3 «a + 28 + 68 


— 


2i + «a + 3 «8 + 58 


— 



(C) 



Solving these equations for the values of the unknown 
quantities and also for their weights, we have 



2i = 


— 8.3 


Pz^ 3.5 


2a — 


- 12.4 


Pz, = 2.5 


28 = 


12.4 


Pz, — 2.5 
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and forming a new set of normals containing 24, we find, on 
solving for that quantity, 

254 == — 24.9 p^ = 1.75 

These valnes of 2^9 ^29 ^9 ^4 ^^ ^.Iso in this case the 
residuals of the observations, and therefore to compute the 
precision measures we have 



V 


. v^ 


pv^ 


8.3 


68.9 . 


206.7 


12.4 


153.8 


307.6 


12.4 


153.8 


307.6 


24.9 


620.0 


.620.0 



«' = 1 



1441.9 = Sjw* 



By (81) ^ = \/^^ = 38 

M.. = 7=z = 20 r^ = 13 

V3.5 

f*«t = M«. = -7=1 =24 r^ = r^ =: 16 

V2.5 

/*«. = -7= = 29 r^ = 20 

Vl.75 

We shall accordingly write for the most probable values of 
the angles of the quadrilateral 



A 


= 101° 


13' 


14" 


± 


13" 


JB 


— 98 


49 


5 


± 


16 


C 


— 87 


5 


27 


± 


16 


J) 


— 77 


52 


15 


± 


20 



In the original solution in paragraph 34 it is evident that 
the results were carried out to a greater number of places of 
significant figures than the character of the observations 
warranted. 



CHAPTER V. 
busobllanbous theorems. 

THE DISTRIBUTIOrf OF ERRORS. 

90. Having developed the processes for the adjustment of 
observations according to the Method of Least Squares, it 
will now be interesting to show how closely the distribution 
of errors found in actual practice corresponds to the theo- 
retical distribution upon which our methods of solution are 
based. 

By formula (36), the probability that the error of a single 
observation will be numerically less than a is 

p = M fV^Vaj (82) 

Let ^ = A a;, .*. <?^ = A dx. Also when a; = 0, 

« == ; and when aj = a, < = Aa = p ^. Sub- 

r 
stituting in (82), 

Values of P for values of the argument — are given in 

Table I. Also for any series of observations this quantity 
P will represent the fraction of the entire number which 
should have errors less than the amount a. Hence if P is 
multiplied by the whole number of observations the result 
will be the number of errors which should be less than the 
limit a. 
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91. Example, Forty measurements on the diameter of 
Saturn's ring were made by Bessel, with the following 
results : — 



M V 


M 


V 


M 


V 


M 


V 


38".91 -.40 


39".35 


+ .04 


39".41 


+ .10 


39".02 


-.29 


39 .32 +.01 


39 .25 


— .06 


39 .40 


+ .09 


39 .01 


-.30 


38 .93 - .38 


39 .14 


— .17 


39 .36 


+ .05 


38 .86 


-.45 


39 .31 .00 


39 .47 


+ .16 


39 .20 


-.11 


39 .51 


+ .20 


39 .17 — .14 


39 .29 


-.02 


39 .42 


+ .11 


39 .21 


-.10 


39 .04 — .27 


39 .32 


+ .01 


39 .30 


-.01 


39 .17 


-.14 


39 .57 +.26 


39 .40 


+ .09 


39 .41 


+ .10 


39 .60 


+ .29 


39 .46 +.15 


39 .33 


+ .02 


39 .43 


+ .12 


39 .54 


+ .23 


39 .30 - .01 


39 .28 


-.03 


39 .43 


+ .12 


39 .45 


+ .14 


39 .03 - .28 


39 .62 


+ .31 


39 .36 


+ .05 


39 .72 


+ .41 



From these the most probable value of the diameter is 
found to be 

B = 39".308 ± 0".022 

the probable error of a single observation being r = 0".136. 
Compare the theoretical and actual distribution of errors 



etween 


0".00 


and 


0".05 


«( 


.05 


u 


.10 


C( 


.10 


a 


.20 


«c 


.20 


(C 


.30 


u 


.30 


<( 


.40 


over 


.40 







In the following table the first column gives the successive 
values of the limiting error a, the second column the values 



a 



of -, and the third column the corresponding values of P. 
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The fourth column contains the differences between the suc- 
cessive values of P, and by multiplying each of these 
differences by 40, the number of observations, we have the 
quantities in column five, which are the numbers of errors 
that according to the theory should fall within the corre- 
sponding limits. Column six shows the actual number of 
residuals occurring between these limits. 



a 


a 
r 


P 


d 


n 


n' 


0.00 
0.05 
0.10 
0.20 
0.30 
0.40 

QO 


0.000 
0.368 
0.735 
1.471 
2.206 
2.942 

00 


0.000 
0.196 
0.380 
0.679 
0.863 
0.953 
1.000 


0.196 
0.184 
0.299 
0.184 
0.090 
0.047 


8 
7 
12 
7 
4 
2 


9 

6 

12 

8 
3 

2 



This is a close agreement considering that the number of 
observations is not very large. Also the number of errors 
greater than the probable error should be equal to the num- 
ber less than it. On counting the residuals we find twenty-one 
less than 0'M36 and nineteen greater. 



THE REJECTION OF OBSERVATIONS. 

92. After a series of measurements have been made, it is 
frequently found that one or two of the observations differ 
widely from the others, and hence it becomes a matter of 
great importance to establish, if possible, some criterion by 
which we may determine whether such discordant observa- 
tions should be rejected or not. We are not concerned here 
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with the question of the detection of a mistake or constant 
error, which a consideration of the circumstances of the 
observations or of the instruments might reveal, but it is 
assumed that there is nothing whatever to guide us except 
the mere fact of the unusual size of the residuals of the 
observations under discussion. To reject an observation 
merely because it differs considerably from the others is 
entirely unjustifiable, while to retain it without any investi- 
gation is a neglect of the evidence furnished by the observa- 
tions themselves. 

The adoption of any rigid criterion based upon the magni- 
tude of the residuals is perhaps more satisfactory from a 
mathematical standpoint than from that of a practical observer, 
and some of the latter are of the opinion that no observation 
should be rejected entirely, even the most widely discordant 
ones being given a certain weight. In this latter case, 
however, the Theory of Probability will furnish a guide as to 
the proper weights to assign to the different observations. 

Of the various criteria that have been proposed, that 
developed by Pierce (see Chauvenety page 568) is perhaps 
the most complete. The derivation and application of this 
criterion is, however, somewhat long and complicated, and 
for all ordinary cases the following simple methods will give 
practically as good results. 

93. Criterion for the Eejection of a Single Doubtful 
Observation. It was shown in (83) that in a series of n 
observations the number of errors numerically less than a 
should be wjP, and therefore the number of errors greater 
than a should be 

n - nP = n(l - P) (84) 

If the value of the expression in (84) is less than one-half, 
the occurrence of an error of magnitude a will have a 
greater probability against it than for it, and hence the obser- 
vation corresponding may be rejected. 
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Accordingly the limit of rejection, ce, of a single doabtfol 
observation is obtained from the equation 



n (1 - P) = - 



or 



P = 



in - I 
tn 



(M) 



94. JExample. Fifteen observations on the value of an 
angle are made. Ought any of the observations to be 
rejected ? 



M 


V 


V* 


«' 


»'* 


v" 


«"« 


2" 23'.90 


..^ 


.30 


.090 


-.41 


.168 


-.33 


.109 


23.76 





.44 


.194 


-.65 


.303 


-.47 


.221 


25.21 


+ 


1.01 


1.020 


+ .90 


.810 






24.68 


+ 


.48 


.230 


+ .37 


.137 


+ .45 


.203 


23.96 


— 


.24 


.058 


-.35 


.123 


-.27 


.073 


24.26 


+ 


.06 


.004 


-.05 


.003 


+ .03 


.001 


24.83 


+ 


.63 


.397 


+ .52 


.270 


+ .60 


.360 


24.07 




.13 


.017 


— .24 


.058 


-.16 


.026 


23.98 


— 


.22 


.048 


-.33 


.109 


— .25 


.063 


24.14 





.06 


.004 


-.17 


.029 


-.09 


.008 


24.40 


+ 


.20 


.040 


+ .09 


.008 


+ .17 


.029 


24.38 


+• 


.18 


.032 


+ .07 


.005 


+ .15 


.023 


24.59 


+ 


.39 


.152 


-.28 


.078 


+ .36 


.130 


24.10 




.10 


.010 


+ .21 


.044 


— .13 


.017 


22.80 


. 


1.40 


1.960 










2° 24 .20 






4.256 




2.145 




1.263 


Jfo 




s«* 




S«'» 




S«"» 
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Using all the observations we find 



Jfo = 2^ 24'.20 r = .6746 i/i:?^ = .37 

^ 14 

By (86), P = |1 = .967 

By Table I, ^ = 3.17 .-. a = 1.17 

As the residual — 1.40 is larger than a, we reject the 
last observation. 

From the remaining observations we now compute a new 
mean value and a new set of residuals* And we find 



M'o = 2** 24'.31 r' = .6746 J^:^ = . 



27 
13 



By (86), P = g = .964 

By Table I, - == 3.11 .-. a = .84 



The third observation may accordingly be rejected. 
From the thirteen observations that remain we find 



M"o = 2° 24'.23 r" = .6745 l/i:?^ = .22 

^ 12 



P = ^ = .962 
26 



- = 3.08 .-. a = .68 
r 



Therefore no more observations are to be rejected. 
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95. The Huge Error. In cases where the number of 
observations is not unusually large, a simple and safe criterion 
for the rejection of a doubtful observation is found in the use 
of the " Huge Error." 

This is an error of such a magnitude that 999 out of every 
1000 errors are less than it and only 1 as large as or greater 
than it. 

Therefore the probability that the error of any given 
observation will be less than the " Huge Error " is .999, and 
from Table I, 

when P = .999, -^ = 4.9 

r 

a = Buge Mrror = 4.9 r 

= 3.3 |Ji (86) 

= 4»1 a.d. 

Then in any limited series of observations, if an error 
greater than the huge error is found, we should reject the 
observation corresponding. 

See also, Holman^ page 30 ; Wr%gh% page 131. 



CONSTANT ERRORS. 

96. Throughout our discussion of the methods of adjusting 
observations so as to obtain from them the most probable 
values of the unknown quantities, all constant errors are 
supposed to have been eliminated before the Method of 
Least Squares is applied in deducing the results. 

If this is not done, and each observation is subject to the 
same constant error, the final result will be affected by an 
equal amount, and in short, the Method of Least Squares is 
not capable of removing or reducing the effect of errors of 
this kind. All that is accomplished by the use of the method 
is to reduce to a minimum the effect of the Accidental Errors. 
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Hence it will be seen that although by increasing the 
number of observations of a given kind we may increase the 
precision, that is, reduce the probable error, of our final 
result as much as we choose, yet we do not in this way 
necessarily increase the accuracy of the determination. 

But if the unknowns can be determined in several ways, or 
under a variety of different circumstances, with various 
instruments, or by different observers, then it is most proba- 
ble that the constant errors of the different sets of measure- 
ments will be grouped about the true values of the unknowns 
according to the exponential law of error. Accordingly a 
combination of such observations will enable us to increase 
not only the precision, but also the accuracy of the final 
result^ the constant errors of the different sets tending to 
cancel each other in the same way that the accidental errors 
of a single set do. 

It is for this reason that determinations of a quantity from 
observations made in a variety of ways are more valuable 
than those obtained merely from different sets of measure- 
ments of the same kind. 

97. The probability of the existence of a constant error 
may often be expressed in the following manner. 

Example. A standard 100 ohm coil is compared with 
a Wheatstone's bridge and the mean result found to be 
100.90 ± 0.20. To find the probability that there is an error 
in the bridge between -j- 0.30 and -{- 1.50 ohms. 

Suppose the result 100.90 ± 0.20 is treated as a single 
observation, and we find by an application of (83) the prob- 
ability that the error of this observation is numerically less 
than 0.60 ohms. ■ 

Here ^ — '^ = 3.00 .-. P = .967 

r ,20 

Hence, as far as is shown by the observations, the proba- 
bility that 100.90 ohms is within 0.60 ohms of the true 
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value is .957. But since it is known that the true resistance 
is 100 ohms, it follows that there is the same probability 
that there is a constant error in the bridge between -{- 0.30 
and -f- 1.50 ohms. 

98. Combination of Determinations having Different Con- 
stant Errors. In case two or more determinations of a 
quantity, together with their probable errors, are obtained, 
the method of combining them so as to secure the best final 
result was considered in paragraph 56, and in Example C, 
paragraph 57. But it was there assumed that all the results 
were subject to the same constant errors, while if this is not 
true the probable errors of the separate determinations bear 
no relation to their weights, and accordingly in such cases 
another process must be adopted. 

To determine whether the difEerent measurements may 
fairly be considered to have the same constant errors we may 
proceed ab follows : — 

Let the determinations of the quantity M be 

^1 ± n (a) 

J^2 ± rt (b) 

and let the difEerence between these results be 

c? = JMi — JM; (c) 

Then the probable error of d is by (66), 

i2 = V^i' + ^2^ (d) 

If d is of such a magnitude that an accidental error as 
great as it may reasonably be expected, we may assume that 
the constant errors of (a) and (b) are the same, and pro- 
ceed as in Example C, paragraph 57. 

But if the probability of making two determinations which 
differ by the amount d is very small, we had best consider 
Ml and M^ to have the same weighty provided there are no 



92 METHOD OF LEAST SQUABES. 

special reasons for regarding one better than the other. The 
final value of Jtf will then be the arithmetical mean of Mi 
and M^ and its probable error will be found by (53). 

99. Example A. An angle is measured by a theodolite 
and by a transit with results 

ByTheodoHte, 24° 13' 36".0 ± 3".l 
By Transit, 24** 13' 24" ± 14" 

What is the most probable value of the angle and its prob- 
able error ? 

Referring to the preceding paragraph, 

ri = 8.1, r, = 14, d := 12, 

and the probable error of d is 



B = )/ 3.1^ + 14« = 14 

Then from Table I the probability that the accidental 
error of a determination will be at least as large as 12 is 
found from 

a 12 

- = — = .86 .-. 1 - P = .67 
r 14 

That is, there is more than an even chance that two such 
determinations of the angle will differ by as much as 12. 

Hence it is fair to assume that the two determinations are 
not affected by constant errors of different magnitudes, and 
they would be combined as in Example C, paragraph 57. 

Example JB. Suppose the zenith distance, M^ oi9k star, 
observed at two different culminations, is found to be 

Jfi = 14° 53' 12".10 ± 0".30 
Jfa = 14° 53' 14".30 ± 0".50 

What is the best final value ? 

Here d = 2.2, R = ^ -09 + .25 = .58 

and for - = -^ = 3.8, 1 — P = .01 
^ .58 
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Therefore the chance that the difference in the two deter- 
minations, due to accidental errors, will be as large as 2.2 is 
only one in a hundred. It is to be concluded then that the 
constant errors of observations at the two culminations differ 
by about 2.2, and as there is nothing to show that one meas- 
urement is more accurate than the other we will give them 
both the same weight and take the mean. Then the best 
value for the zenith distance is 

Mo = 14° 53' 13".20 ± 0".74 

..o CO, . 12".10 + 14".30 

For Mo = W 63' 4- -^ 

^ 2 

= 14** 53' 13".20 



. / 2.42 
By (53) ro = .6745^2"^^ = .74 

For a more extended treatment of this subject see Johnsorij 
" The Theory of Errors and Method of Least Squares," 
chap. vii. 

THE WEIGHTING OF OBSERVATIONS. 

100. In case the relative worth of observations is not 
settled by methods already discussed, the proper weight to 
assign to each quantity in the final adjustment can only be 
determined from a full knowledge of all the circumstances of 
the measurements. Even then considerable experience in 
the particular work in hand is required before the best values 
for these weights can be assigned. The weight given to a 
quantity should never be considered final, but always subject 
to revision whenever new information with regard to the 
quantity is obtained. Thus an observation which at first is 
supposed to deserve a high degree of confidence is often 
found on later investigation to possess very little value, and 
vice versa. 

See also Wright^ page 118. 
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OTHER LAWS OF ERROR. 

101. Although in the great majority of cases the distribn- 
tion of errors follows the exponential law thus far considered, 
there are a few special cases in which somQ of the suppositions 
made in deriving that law do not hold, and hence for the 
adjustment of such observations the corresponding special 
laws of error must be determined. 

For instance, in applying the exponential law we assume a 
large number of observations, that each observation is sub- 
ject to the same law of error, that small errors are more 
likely to occur than large ones, and that positive and negative 
errors are equally probable. Now it is easy to conceive of 
cases where only positive errors can occur, or where the 
probability of the occurrence of a small error may not be 
greater than that of a larger one, etc. If we can determine 
the different sources of error in any case and the relative 
effect of each upon the quantity sought, we shall arrive at 
the law of error for that particular set of observations. The 
case of most conmion occurrence is the following. 

102. Suppose all errors between the limits a and — a 
are equally probable, and that there are no errors beyond 
these limits. Then if y r= ^(a) is the equation of the 
Curve of Error, and its area is represented as in para- 
graph 18, we have 

X </»(«) dx =: 1 (a) 

a 



or 



2 <l>(x) C^'dx = 1 (b) 



since by the supposition made <f>(x) must be a constant. 
Integrating and solving for ^(a), we have 



y = ^(x) = -1- (-87) 

2 a 
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To find the Mean Error we have by definition as in para- 
graph 52 






7? dx 



a 



• • 



ii = « 



v» 



The Probable Error is derived from the equation 

CfKx) dx =z 

I Cdx = 1 

a Jo 9. 



.•. 



^ a 



Finally, for the Average Deviation we have 



a.d. =z i X <l> (x) 



dx 



= -/' 



X dx 



•\ 



a.d. = 



a 
i 



And the Curve of Error has the form 



(88) 



(89) 



(••) 
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That the average deviation and probable error are in this 
case equal to one half of a may also be seen from the defini- 
tion of these quantities. 

Example, In taking a logarithm from a four place table, 
what is the probable error of the mantissa ? 

In this case the maximum error is .00005, and all 
errors between — .00005 and .00005 are equally probable. 
Therefore 

r = .000025 

103. The only other special case of common occurrence is 
that in which the error of a quantity is due to two sources, 
each of which can with the same probability assume all 
values between a and — a. Here it may be shown that 
the curve of error consists of two straight lines whose equa- 
tions are 

la — X , ia + X ,„. 

y = — and y = ■ (91) 

2 

Also |i« = -- a% r = (2 - \^) a (M) 

o 

For a more extended discussion of special laws of error, see 
Wrighty paragraphs 31 to 39. See also Example 151. 

104. Contradictory Observations. Suppose three observa- 
tions of a quantity give results 55, 56, 91. It is obvious that 
to take the arithmetical mean as the most probable value in 
this case would be contrary to the evidence furnished by the 
measurements, while in so small a number of observations it 
would not be allowable to reject the value 91 entirely. With 
such a series of observations no satisfactory solution can be 
obtained, but probably the best thing to do would be to take 
the value which has as many observations less than it as it 
has greater, or 56. 



CHAPTER VI. 

GAUSS'S METHOD OF SUBSTZTUTION. 

105. The most laborious part of the application of the 
Method of Least Squares to the adjustment of observations 
consists in the formation and solution of the <' Normal Equa- 
tions," and this labor increases enormously with increase in 
the number of observations, of unknowns and of conditions. 
It is not at all unusual to find that the adjustment of a single 
get of observations takes several weeks, even with all the aid 
that can be obtained from tables of logarithms, of squares, of 
products, and of reciprocals, and also from the use of calcu- 
lating machines. In general the computations can be per- 
formed more rapidly and with less fatigue by using a machine 
or a table of products than by using logarithms, but a combi- 
nation of methods is often desirable. 

106. In any case however it is of the utmost importance 
that the formation and solution of the normal equations 
should be effected in a systematic way, and that as far as 
possible checks on the numerical work be carried along in 
the computations. By a slight amount of additional work 
checks upon the results at successive stages of the solution 
may be obtained by means indicated in the following dem- 
onstrations, while the most satisfactory form for the solution 
of the normal equations is given by the << Method of Substitu- 
tion" proposed by Oausa. 

This method will now be explained; and it is to be observed 
that it is customary in this subject to enclose a quantity in 
brackets when the sum of a number of quantities of the same 
kind is to be denoted. Thus [a^] means the same as ^ab. 

For the sake of simplicity in demonstration it will be 
assumed that the observations are all reduced to weight unity. 
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107. Checks on the Pormation of the Vormal 
If, as in paragraph 76, we take for 



OBSERVATION EQUATIONS 



«f2i 4- 5,2, -f- 



^n^l + K^ + 



• • 3'a«g + »W, = 






• • • 



• • 9n^q -\- ^^ = % 



we shall have for 



(A) 



NORMAL EQUATIONS 

[aft] 2i -f- [ftft] 2, -|- . . . Ibq^ Zq -}- [ft^] 









(B) 



Let 



«i + *i + • • • 3'i + '^i = *i 

a, -|- 5, -|- • . , g'j -}- m, = «2 

••• ••• ••• 

••• ••• ••• 



(C) 



... [a] + [ft] + . . . [g] + [m] = [«] (D) 

Multiplying the first of (C) by mi, the second by m^ . . . 
and adding, there results 

lam'] + Ibm] + . . . [ffw.] + [mm] = \_sm] (M) 

Next multiplying each of equations (C^ by its a and 
adding, and then each by its h and adding, etc., we have 
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[aa] + [aft] + • • • C«ff] + [«^ = [««] 
[aft] + [ftft] + . . . [ftff ] + [6w] = [»«] 

••• ••• ••• ( vvl 

• •• ••• ••• 

[a* ] + [ftff] + • • • [ff«] + [«^] = [ff«] 

Equation (93) will be satisfied if the absolute terms in 

the normal equations are correct, and equations (94) when 

the coefficients of the unknown quantities are correct. These 

check the formation of the normal equations. 

108. The Keduoed Normal Equations and the Elimination 
Equations. 

The value of z^ in terms of the remaining unknowns, 

derived from the first of equations (B), is 

Substituting this in the remaining n — 1 equations, they 
become 



• • • 



And letting 






OT 



• • • 



Iqm, 1] = [jm] - ^^^ 
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the above equations take the following form, which, being 
the same as that of the original normal equations, they are 
called the 

FIRST REDUCED NORMAL EQUATIONS. 

[W, 1] «, 4- [^c, 1] «3 + • • . [^^, 1] 2^g + [^^> 1] = ^ 
[bCy 1] 22 + IcCy 1] gg + . . . Icq^ 1] 2g + [cm, 1] == 

• •• ••» ••• I ^^9 

• •• ••• ••• 

{_bq, 1] 2a + Icq, 1] 2, + . . . Iqq, 1] z^ + [^m, 1] = 

An inspection of equations (F) will render it easy to form 
a rule for writing out any one of them. 

Now by means of the first of equations (6), eliminating 
^2 from each of the others in the same way that 2i was 
eliminated from the normal equations, there results the 

SECOND REDUCED NORMAL EQUATIONS 
Icc^ 2] 28 + • • • [^^5 2] Zq + Icm^ 2] = 



(H) 



Icq, 2] 23 + . . . Iqq, 2] Zg + Iqm, 2] = 

In which 

[*c, 1] Ibc, 1] 



[cc,2] = [cc,l] - 



[W, 1] 



(I) 



Cm. 2] = [,„..!] _ [%1K^^>13 



Continuing this process we shall finally arrive at the single 
equation 

Iqqy q-1] z^ + [^m, q-1] = (J) 

from which the value of Zg is determined. 
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The value of Zq.i will then be obtained by substitnting the 
numerical value of z^ in the first of the preceding set of 
equations, and so on, until finally z^ is obtained from the 
first of the original normal equations. The equations from 
which the unknowns are actually determined are then the 
following, called the 

ELIMINATION EQUATIONS. 

\aa\ »i + [^J^] «2 + • • • [/iq] »q + [am] = 
[66, 1]»2 + ...[&«, l]«<z + [6w,l] = 

• •• ••• ••• f V V J 

• • • • • • 

Iqq, q-1] «« + [gm, q-1] = 

It may be seen from the rule in paragraph 81 that 
[<qq^ q— 1] is the weight of 2?g, and the weight of any 
unknown might be found at the same time as its value by 
making it the last in the order of elimination, but except in 
special cases the weights had best be obtained by the general 
process of paragraph 115. 

109. Check on the Solution of the Normal Equations. 
Multiplying the first of the Observation Equations (A) by 
mi, the second by ms, . . . and adding the results, we have 

[mv] = [am] z^ + [^^] 2^2 + • • • [S'^] ^q + [rnin] 

But in equation (T), paragraph 82, it was shown that 
[mv] = [vv]. Therefore 

[yv] = [am] z^ -}- [^] ^2 + • • • [j^] ^q + [w*^] 

Substituting in this the value of z^ from the first of (96), 
we get the result 

[vv] = [hm^ 1] 22 + [cm^ 1] 23 + . • • [S'^j 1] «g + [w»w, 1] 
in which 
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[oJ] [am] 



[ftm, 1] = [«m] - 



[«.] 



(K) 



[mm, 1] = [mm] _ MM 

being similar in form to eqoationB (F). 

Next eliminating 2, in a like manner, we get 

\vv] = \cmy 2] 2?g -|- — [3'^> 2] «^ -}- [mm, 2] 

and oontdnuing this process it finally appears that 

110. Arrangement of the Computations. In computing 
the coefficients that appear in the ** Auxiliary" or <^ Reduced 
Normal E<][uations" it is most convenient to arrange the work 
in tabular form. The arrangement of the solution will be 
illustrated for an example containing four unknowns but it 
will be evident that the process can be extended to cover any 
case. 

Let the Observation and Normal Equations be represented 
by equations (A) and (B), there being only four unknowns 
^19 ^S9 ^S9 ^49 Ai^d arrange a table as on the next page. In this 
scheme of solution the upper lines of the rows in the first 
compartment contain all the quantities that appear in the 
Normal Equations, together with [mm] and the quantities 
in the column headed a which are used in checking the 
results in accordance with equations (93) and (94). The 
other compartments contain the corresponding quantities 
for the Reduced Normal Equations, and the first line in 
each compartment gives the coefficients in the Elimination 
Equations. 
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Scheme A. — Solution of the Normal Equations. 





a 


» 


e 


<f 


» 


9 


a 


logM 


log[«»] 


M 

log[«] 


[«l] 
log[«l] 


[«m3 
log [am] 


iog[«3 


b 


logJ» 


[W] 
log^.C-*] 


logJ»[ae] 


[W] 

A,[ad} 

log^Corf] 


[*«] 

-i» [«»] 
logJ»[a«3 


J»[a«3 
logJ»[a«] 




log^ 


1 


logJ,C«] 


[«£] 

J.[orf] 

Iog^[«f] 


[«.3 

^^cC-3 
logJ.[am3 


^«[«3 
logJ.[«3 




logJtf 




logJ^[«rf] 


C*»3 
- A4[am^ 
, log -4rf [am'] 


[A3 ^ 

log^rf[a»3 




log^ 




[m«3 

A^[am\ 

log ^» [am] 


;[~3 

4.[«3 
iogj,[«3^ 






log[M.l] 


log[4e.l] 


Iog[«;i] 


[fcmi] 

log[ftm.l] 


[*«.i3 
iog[&.,i3 




logB. 




log f?. [fc, 1] 


^.CHi] 

log^[W,l] 


[«».13 

*[*«m] 

log B^ [*m, 1] 


[e#.13 

B. [«^ 1] 

iogB,[&«,i3 




Iog'4 




log/7«[Hl] 


[rfm.l] 

i;«[i>mi] 
logii^[fcml] 


CA^i3 

-R«[«M3 
rogB^[AM3 




logoff. 




' [mm,l] 

^[HO 
log B, [&m 13 


[tw. 13 

A. [AM] 

log ^[A.. 13 






[ce,2] 
log[...2] 


log [at. 2] 


[«»»2] 
log [cm. 2] 


C«^23 ' 
log[c*,23 




logCj 




[ritf.2] 

C*[«*23 

logC4[«i;23 


C4[cm.%^ 
logCrf [cm, 2] 


[A, 23 
' O4 [cm, 23 . 
logO'rf[c«,23 




logC. 




[mm, 2] 

C,[«m,2] 

log C;, [cm, 2] 


[m.,23 

CL.[c23 

log C, [cf. 2J 






[i«.83 
log[iM.83 


[rf«,83 
log[rfm.83 


[i^t.83 
Iog[rf*,8] 




logA, a l6f-«« 




[mm, 83 
logZ),[rf«,83 


[««,83 

A. [A. 8] 

logA.[A.83 


- 


«. - w = 


[mm, 4] 


[»w.4] 
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The logarithms of the quantities in the first row of each 
compartment are also written in, and from these by proper 
subtractions are obtained the logarithms in the margin, where 



A. - f°*^ 


[ac] lam} 
laay •" -^ laa-\ 

„ Ibc, 1] „ Ibm, 1] 

■"" — Ibb, 1]' • • • ^o — [M, 1] 

• • • • • • 1 




• • • • • • 

Idm, 3] 
** ^-^ Idd, 3] 



(L) 



Now in each compartment adding the logarithms at the 
margin to each of the logarithms in the first row of that com- 
partment we obtain the corresponding logarithms written in 
the other rows. The numbers represented by these logarithms 
are next written above them, and if each of these quantities 
is then subtracted from the one above it the result will be 
the corresponding quantity in the compartment below. Some 
of the squares in each compartment are left vacant as the 
quantities belonging to them have already appeared above. 

111. Application of Checks. Also, by (93) and (94), in 
the first compartment the quantities in the first lines of the 
last column should be equal to the sum of all the quantities 
in the first lines of the corresponding rows plus the quantities 
similarly situated above the first terms of the rows. Similar 
checks will apply in each compartment; for if from the 
second of equations (94) we subtract the product of the first 
equation multiplied by A^ we have 

Ibb, 1] + [ftc, 1] + . . . [^, 1] = [bs, 1] (97) 

In the same manner we may show that a corresponding 
check holds throughout, so that finally we shall have 

[mm^4] = [w«, 4] (M) 
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The last compartment of the table is added to give this 
final check and the value of [yv\ in accordance with equa- 
tion (96). 

If l^e multiplications and divisions are simple or if a table 
of squares or products or a computing machine is used the 
logarithms will of course be omitted from the scheme of 
solution. 

112. Example, In order to illustrate the systematic for- 
mation of the coefficients that appear in the Normal Equations 
as well as the solution of the latter by the above method we 
will take the 



OBSERVATION EQUATIONS 



«1 + 2a + 28 + «4 = 0.1 
21-4- 2a — 28 — 2* = 0.6 
gj ^ 222 — gg -|- ^4 = 0.1 



2i 


— 2a — 2^4 = 


0.3 




2a - 2, + 24 _ 


— 0.1 


2i 


- 2, _ 


0.4 



First form a table containing the coefficients in these equa- 
tions and also the sums 8, As a first check the sum of the 
quantities in colunm 8 should be equal to the sum of all the 
quantities in all the other columns. 



I. 


Coefficients in 


THE ObSEBVATION EQUATIONS. 


No. 


a 


h 


C 


d 


m 


8 


1 


- 1 


1 


1 


1 


- .1 


1.9 


2 


1 


1 


— 1 


- 1 


- .6 


— .6 


3 


1 


— 2 


1 


1 


— .1 


- 1.1 


4 


1 


— 1 





- 2 


- .3 


- 2.3 


5 





1 


^ 1 


1 


.1 


1.1 


6 


1 





- 1 





.4 


- .4 


Sum. 


3 





- 3 





- 1.4 


— 1.4 
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From these we now compute the coefficients in the Normal 
Equations (B), and also the necessary quantities for the 
check equations (93) and (94). 

II. Coefficients in the Normal Equations. 



Ifo. 


aa 
1 


ab 
-1 


ac 
-1 


ad 
-1 


ctm 


as 


bb 
1 


be 
1 


bd 
1 


bm 


1 


.1 


-1.9 


-.1 


2 


1 


1 


-1 


-1 


-.6 


- .6 


1 


-1 


-1 


-.6 


3 


1 


-2 


-1 


1 


-.1 


-1.1 


4 


2 


-2 


.2 


4 


1 


-1 





-2 


-.3 


-2.3 


1 





2 


.3 


5 














.0 


.0 


1 


-1 


1 


.1 


6 


1 
5 



-3 


-1 
-4 



-3 


-.4 


- .4 




8 




1 



1 


.0 
-.1 


Sum 


-1.3 


-6.3 



JVo. 


bs 


cc 

1 


cd 
1 


cm 
-.1 


cs 


dd 
1 


dm 
-.1 


ds 


mm 
.01 


ms 


1 


1.9 


1.9 


1.9 


-.19 


2 


- .6 


1 


1 


.6 


.6 


1 


.6 


.6 


.36 


.36 


3 


2.2 


1 


-1 


.1 


1.1 


1 


-.1 


-1.1 


.01 


.11 


4 


2.3 








.0 


.0 


4 


.6 


4.6 


.09 


.69 


5 


1.1 


1 


-1 


-.1 


-1.1 


1 


.1 


1.1 


.01 


.11 


6 


.0 


1 
5 






.4 

r 

.9 


.4 




8 


.0 
1.1 


.0 


.16 
.64 


.16 


Sum 


6.9 


2.9 


7.1 


1.24 



If the coefficients in the Observation Equations are large, 
so that it becomes convenient to use logarithms, other tables 
corresponding to I and II would be formed to contain these 
logarithms. 
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Substitutiiig these quantities now in the general tabular 
scheme of paragraph 110 we have the results on the following 
page. The work of the first compartment is performed 
without the use of logarithms, as the numbers are simple. 
The quantities in the last column should all be zero accord- 
ing to the check equations, what small differences there are 
being due to the rejection of figures beyond the third place 
in the decimals. The decimal points in logarithms to which 
correspond negative numbers have been replaced by the 
letter n. 

The demonstrations that have been made now enable us to 
see at once from an inspection of the results in this table 
that 

S4 = 0.238 

Pz, = 1.6 
[vt>] = .007 

Therefore substituting in equations (73) and (74) we 
have 

fi^ = Ml r^ = .032 

If the two values of [vt;] obtained in the solution had 
differed at all we should have taken the mean of the two. 

113. If the Elimination Equations (95) are divided by 
[ao], [^5,1], [cc, 2], [cf<?, 3], respectively, they become 

Zi + Ab»2 + Ac^z + Ad»4 + Afn = 9 

«8 + Caz, + C« = ^"^ 

And the solution for the unknowns can be effected most 
conveniently by arranging the computations in the manner 
illustrated on page 109. 
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ScHBMB A. Solution of 
a b c 



THB Normal Equations. 
d m 8 B 



a 


5 


-3 


-4 


-3 


-1.3 


-6.3 


00000 


b 
At, 


-.6 


8 
1.8 


1 

2.4 


1 

1.8 


- .1 

.78 


6.9 
3.78 


c 


-.8 




5 
3.2 



2.4 


.9 
1.04 


2.9 
6.04 


d 
Aa 


-.6 




8 
1.8 


1.1 

.78 


7.1 

3.78 


m 

An 


-.26 




.64 
.34 


1.24 
1.64 


b 




6.2 
0.7924 


-1.4 
0,1461 


- .8 
9,9031 


- .88 
9„9445 


3.12 
0.4924 









c 
Be 


9,3537 




1.8 
.316 
9.4998 


-2.4 

.181 
9.2568 


- .14 

.199 
9.2982 


-2.14 
- .705 
9^8479 


d 
Ba 


9,1107 




6.2 
.103 
9.0138 


.32 
.114 
9.0552 


3.32 
- .403 
9„6049 


m 
B„ 


9,1521 




.30 
.125 
9.0966 


- .40 

- .443 
9„6463 


c 




1.484 
0.1715 


-2.581 
0,4118 


- .339 
9^5302 


- 1.436 
0^1569 


1 
1 

1 


d 


0,2408 




6.097 

4.489 

0.6521 


.206 

.689 

9.7705 


3.373 

2.496 

0.3972 


m 
<7« 


9,3587 




•176 
.077 

8.8889 


.043 

.328 

9.5156 


d 


9,3769 — log -84 


1,608 
0.2063 


- .383 
9„5832 


1.227 

0.0888 


2 







.098 

.091 

8.9601 


- .285 

- .292 
9„4657 


tn 


0,238 = 24 [»w] 


= .007 


.007 


.007 
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Scheme B. — Solution op the Elimination Equations. 



--Z>m 


-Q24 




~ ^d24 

- ^« 28 

— ^6 2a 


24 
l0g2J4 


28 
log 28 


2a 
log 2a 


2i 


logCa 
log^d 
log^d 


log^. 
^ogAc 


log ^6 




log Aj, 2, 


log ^e «8 

log Ac 28 


log Caz^ 
log-Bd 24 
log Aa 24 




\ 



114. Filling oat this table for the example just solved, we 
have 



.238 


.228 


.142 


.260 




.414 


.031 


.143 






.145 


.514 
.191 


.238 


.642 


.318 


1.108 


9.3769 


9.8075 


9.5024 




0,2403 


9,3537 


9,7782 




9,1107 

9„7782 


9,9031 






9,2806 


9,1612 
9,7106 


9,6172 

8,4876 








ft 




9,1551 
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115. The Weights of the Unknowns. In order to deter- 
mine the precision measares of s^, z^ and Sj, it would 
next be necessary to compate the weights of the latter qaan- 
tities. The demonstration of the processes by which these 
weights may be found will not be taken up here, as the best 
method to adopt varies a good deal with the character of the 
example, but a statement of the results in the general form 
of solution will be given. 

By treating the Elimination Equations in a way similar to 
that used in deriving equation (E) of paragraph 76 from 
equations (B) of the same paragraph, we may show that 

z, ^ B^ + C^p, + JD^p, = (100) 

where the a's, ^'s, y's, are determined from the equations 
^d + Baa^ + ^do, + a, = ^d + ^d A + A = 

^c+^cai + a, = A+A = (101) 

^6 4- ai = O;, + ya = 

Then by an application of the principles of Rule 1, para- 
graph (77), it may be shown that 



1 




1 


P^^ 




[aa] 


1 






Ph 






1 






Ph 






1 







a,» 



a,* a,» 



~r rj.!. n ~r r^^ »1 "r 



[66, 1] ^ Ice, 2] ^ Idd, 8] 



1 . i8,« . fi 



~r r^^ »n "T 



a A 2 



[66, 1] ^ ice, 2] ^ idd, 3] 



1 _, Vo 

lee, f] "*" Idd, 8] 

1 

Pxt ~ Idd, 8] 



, (W*) 
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These equations can of coarse be extended to cover any 
number of unknown quantities, and tabular schemes for the 
computations of the a's, ^'s, y's, . . . and of the weights 
Px^f Pz^y Pz^i ... can readily be arranged. 

For a general demonstration of these results, and also for 
a discussion of special methods of solution, consult — 

clbAn«on, <*The Theory of Errors and Method of Least Squares,** 

chap. ix. 
Wright^ ** Treatise on the Adjustment of Obsenrations," chap. iy. 
Cliauvenetj ** Spherical and Practical Astronomy,** pp. 580-649. 



THE METHOD OF CORRELATIVES. 

116. The method of adjusting "Conditioned Observa- 
tions " explained in paragraph 33 is perfectly general, but 
where there are many conditions to be satisfied the solution 
is apt to be very laborious. For the case that occurs most fre- 
quently in practice, in which the observations are direct and 
equal in number to the number of unknown quantities, the 
process of solution devised by Gauss and called the " Method 
of Correlatives" is the most convenient. This method is 
derived as follows : — 

Let q observations, -3/1, -3/^, . . . Mq^ of the respective 
weights />i, />j, . • . jp^, be made directly upon the values 
of q unknown quantities, and let the most probable values 
of the unknowns be 

Zi = Jfi -|- Vi, 2a = -3/i -|- Vj, . . .Zq =z Mq-}- Vq. 

Where Vi, ^2, . . . Vg, are the most probable corrections to 
apply to the observed values as well as in this case the 
residuals of the observations. 

If the n' condition equations are not linear they may be 
reduced to that form by the method of paragraph 44, so that 
we may assume for our 



(A) 
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CONDITION EQUATIONS 

^1^1 + ^a^« "["••• ^q^q + W»i = 

biVi -|- ^»^a + • • • ^q'^q + »W, = 

• •• ••• ••• 

• •• ••• ••• 

llVi -f- ^l?s 4" • • • ^q^q + ^n' = 

In which the quantities m^, m,, . . . m^s wonld all be zero 
if the observations were exact. It is to be observed that the 
coefficients a, d, . . . / are not arranged in the same order in 
these equations as they are in the observation equations of 
paragraph 107. 

The values of i^i, Vs, ... Vq^ must be determined so as to 
satisfy the above equations and also by the principle of Least 
Squares, so as to make 

JPi^i* + Pi^i^ + • • • Pq'^q^ = * minimum. 
Corresponding with a minimum value of this last we have 

PlVidVi -f- P%V^dV2 + • • . PqVqdVq = (B) 

for all possible simultaneous values of dvi^ dv^ . . . dvq ; that 
is, for all values which satisfy the equations, 



(Zidvi -}- a^dv^ -}-... dqdvq = 
hidvx + ^a^v, + . . . hqdvq = 



(C) 



lldVi -f- h<^V2 + . . . iqdVq = 

obtained by differentiating equations (A). 

Therefore, denoting the first member of (B) by R and the 
first members of (C) by /S'l, /^j, ... jS^'^ it will be nec- 
essary that 

H _ kiSi - A?,/^, - . . . k^^^n' = (D) 

where k^ A;,, ... k^i , are undetermined coefficients. 
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This last eqaation will be satisfied if the coefficient of each 
difEerential in it is made equal to zero, that is, if 

p^Vi = kiCti -j- k^bi -|- • • • ^n'h 



(103) 



JPq^q *"~ "'i^q "7" ^2^q ~t~ • • • t^n'^q 



All that remains therefore is to find values of Vi, v,, ... Vq 
and kiy k^ ... kn,y which will satisfy simultaneously equa- 
tions (A) and (103), and that this may be done is easily 
seen from the fact that we have the same number of equations 
as unknowns. 

Substituting the values of Vi, v^ ... Vq from equations 
(103) inequations (A) we have the following : 



_ «, €141 , - — ab , , ^ ^^ . A 

fciS— + k^^— + ...*JnH— + mi = 

• •• ... •.. I aVv J 

• •• ••. ••• 

The solution of these equations gives at once the values of 
kiy A?2, ... kn'y which are called the "Correlatives'* of the 
Condition Equations. The values of Vi^ t?2, . . . Vq are then 
found by substituting the values of the k^s in equation (103). 

117. As equations (104) are of the same general form as 
a set of Normal Equations, Gausses Method of Substitution 
can be advantageously employed in the solution. 



114 METHOD OF LEAST SQUABES. 

When there is but a single equation of condition the second 
members of equations (103J reduce to their first terms, and 
equations (104) reduce to the single equation 

A?iS— -f mi = (105) 

and th(B values of t^i, v^ ... Vq in (103) become 

JPI Pi 

2 ~~zr 2 



p p 



(!«•) 



1» 

It is from these results that the rules in paragraph 35 are 
derived. 

118. Example. Suppose we have given the observations 



M, - 


2.02, 


weight 3 


M, - 


4.13, 


« 2 


M, = 


2.52, 


" 5 


M, - 


2.67, 


" 7 


M, - 


2.84, 


" 4 



(a) 



and let the most probable values of the unknowns be repre- 
sented by 

2i = itfi + Vi, 2j := Jfi + ^a» • • • ^6 = -^ + ^6 (^) 

Also suppose that the unknowns are subject to the conditions 

«, ^z^ = 1.6 ^^> 
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Then expressing these conditions in tenns of the corrections 
by means of (a) and (b), we have the 



CONDITION EQUATIONS 

^1 + ^2 + «8 + ^4 + ^6 + -18 = ^ ... 

Va - V4 - .04 = ^^) 

Referring to paragraph 116, we see that in this example 

n' = 2, mi = 0.18, m^ = — 0.04 

For the purpose of computing the coefficients in equations 
(104) we next arrange the following table. 



p 


a 


h 


aa 
P 


ab 
P 


bb 
P 


3 
2 
5 

7 
4 





1 

- 1 



1 

"s 

1 
2 
1 
5 
1 
7 
1 
4 




1 
2 



_ 1 

7 






1 
2 



1 
7 










599 
420 


5 
14 


9 
14 
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Substitating these results in equations (104), we have 

599 5 

*i H A?a + .18 = 

420 ^ ^ 14 * ^ 

5 9 
— A?i H A?, - .04 = 

14 * ^ 14 

Solving, A?i = - .1647 

A?, = .1537 

Then from equations (103) we get at once 



(e) 



(f) 



i>i = — .0549, i>4 = — .0455, 

V, = — .0055, vs = — M12. (g) 

v, = —.0329, 

And by substituting these results in (b) we can obtain the 
values of Si, 2,, 2,, s^, 25. 

The above is the solution of Example 70, page 127. 
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1. An urn contains five black balls, three red balls and 
two white balls. If three balls are drawn from the urn what 
different combinations may result, and what is the probability 
of each ? 

2. In a single throw with a pair of dice what is the 

5 
probability that neither ace nor doublets will appear ? 

y 

3. Foui* cards are drawn from a pack. What is the 
probability of getting four aces ? Of getting one of each 
suit? 

4. From a lottery of thirty tickets, marked 1, 2, ... 30, 
four tickets are drawn. What is the probability that num- 

2 

bers 1 and 15 will be among them ? 

145 

5. Find the odds against the appearance of 7 or 11 in a 
single throw with a pair of dice. 7 : 2 

6. I toss up n coins. What is my chance of getting just 
one head ? 

7. In a single throw what are the relative chances of 
throwing 9 with two dice and with three dice ? 24 : 25 

8. From 2 n counters marked with consecutive numbers 
two are drawn. What are the odds against having an even 
sum ? /^ : /i — 1 

9. In two trials with a single die what is the probability 
of throwing (a) an ace the first time only? (b) at least one 
ace? 

10. Find the probability of throwing doublets one or 

91 
more times in three trials with a pair of dice. - — 

216 
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11. Find the probability of throwing exactly three aces 

in five trials with a single die. 

3888 

12. A certain stake is to be won by the first person who 
throws 5 with a die of twelve faces. What is the chance of 
the sixth person ? 

13. A and B play chess. A wins on the average two 

games out of three. What is A's chance of winning just 

80 
four games out of the first six ? 

^ 243 

14. A and B shoot alternately at a mark. A hits once in 
n times and B once in n — 1 times. Find their chances of 
first hit, and the odds in favor of B if A misses on his first 
shot. Even. 7» : /^ — 2 

15. In how many trials will it be a wager of 4 to 3 that 
double five will be thrown with a pair of dice ? 30 

16. Find the probability of throwing one and only one 

ace in two trials with a single die. — 

18 

17. If I have three tickets in a lottery of four prizes and 

eight blanks, what is my chance of drawing a prize ? — 

55 

18. Find the probability of throwing at least four aces in 

six trials with a single die. 

23328 

19. On an average seven ships out of eight return to port. 
Find the chance that out of five ships expected at least three 

•11 ^ 16121 
will return. ^ 

16384 

20. In a lottery containing a large number of tickets, 
where the prizes are to the blanks as 1 : 9, find the chance 

of drawing at least two prizes in five trials. 

100000 
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21. In a purse are ten coins, all nickels except one which 
is a five-dollar gold piece ; in another are ten coins, all 
nickels. Nine coins are taken from the first purse and 
placed in the second, and then nine coins are taken from the 
latter and placed in the former. If you now had your choice 
which purse would you take ? 

22. A and B engage in a game in which A's skill is to 
B's as 2 : 3. What is A's chance of winning at least two 
games out of five ? 

23. If A's skill at a certain game is double that of B, 
what are the odds against A's winning four games before B 
wins two? 131:112 

24. A party of twenty-five take seats at a round table. 
What are the odds against any two specified persons sitting 
next to each other ? 

25. A has three shares in a lottery in which there are three 
prizes and six blanks. B has one share in another where 
there is but one prize and two blanks. What are their 
relative chances of getting a prize ? A : B = 16 : 7 

26. Expand through the terms involving A* and A;*, the 
expression 

-1— + (y + *)• • 

aj -f- A 

When 35 is 1 and y is -, does — h y* increase or 

2 X 

diminish when x and y begin to increase at the same rate ? 

27. Given f {Xy y) = «*(« -}- y)*, expand the expres- 
sion (05 + A)a (a + y + ky. 

28. Find the value of logio a -f~ ^^^ ^' when a = 1001 
and b = 0.1°. Give the result first to five places and then 
to seven places of significant figures, in each case without the 
aid of tables. 4.000433 

29. Transform to the new origin, (2, ~ 3, 1), the equa- 
tion, a5* + y*-)-»«— 4aj-f6y-2g-ll = 0, the 
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axes remaining parallel to the original ones. The equations 
of transformation are x = 2 + *'> y = — 3 + y'> 
2=1 + 2'. «« + y« -f 2« = 26 

30. Find the minimum value of a* -f" «y + y* — aaj — by. 

(ab - g' - b^) 

3 

31. Find the values of as, y and z that render a maximum 
or a minimum the function a^+y^-f-g^-f-a — 2 z — xy. 

32. Find the values of x and y that render a maximum 
or a minimum the expression sin x -|- sin y -|~ ^^^ (^ ~h y>)* 

33. Find the co-ordinates of a point the sum of the squares 
of whose distances from three given points, (a^, yi), (a^, yg)* 

(3585 ys)* 18 a mmimum. x = ^ ' — ' ^^ 

3 

34. Given the volume, a*, of a rectangular parallelopiped, 
find its shape when its surface is a minimum. 

35. Find the volume of the greatest rectangular parallelo- 
piped that can be inscribed in the ellipsoid 

^ y* 2* 1 ^^^^ 

-, + j2 + -i- Y^ 

36. In the Physical Laboratory apparatus for illustrating 
the Estimation of Tenths, the reading of the micrometer head 
for a certain setting is computed to be 2.3038, which may be 
taken as exact. Setting the apparatus by the eye the follow- 
ing readings are obtained : — 

2.314 2.324 2.310 2.519 2.326 
2.320 2.302 2.313 2.305 

What are the accidental and real errors and the residuals ? 
Are there any constant errors or mistakes ? If there are 
constant errors are they of the first, second, or third class ? 
How do you tell ? 
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37. Are the following observations such as to call for an 
application of the Method of Least Squares in their adjust- 
ment? 

aj+y— j5-j-^= ^ 3a5 — y— z — w=l 

05 — 2y4-22 + ^= — 1 X — y -{- ^z-\-^u = 9 

38. In the case of direct observations, what other quanti- 
ties besides the arithmetical mean might reasonably be 
assumed to give plausible values of the unknowns ? Why is 
the arithmetical mean preferred to these ? 

39. Find the most probable value of a quantity M from 
the observations 

216.58 



216.27 


216.16 


216.04 


216.19 


216.44 


.29 


.43 


215.99 


.39 


.51 


.33 


.09 


216.23 


.14 


216.94 



Also test the result by finding the sum of the residuals. 

M = 216.251 

40. In the determination of a certain wave length, Row- 
land made the following observations. Find the most prob- 
able wave length. 

4.524 4.515 4.513 4.507 4.501 4.485 4.517 4.493 4.505 
.500 .508 .511 .497 .502 .519 .504 .492 

4.5055 ± 0.0017 

41. Ten measurements of the density of a body gave 
results as follows. Find the most probable density. 



9.662 


9.664 


9.677 


9.663 9.645 


.678 


.659 


.662 


.680 .654 

9.6689 ± 0.0022 



42. In a triangulation of the U. S. Coast Survey an angle 
was measured twenty-four times with results 
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116« 43' 44".45 


48".90 


47".40 


47".85 


51''.75 


50 ^5 


49 .20 


47 .75 


50 .60 


49 .00 


50 .95 


48 .85 


51 .05 


48 .45 


52 J5 




51".S0 


49".05 


46".75 






51 .05 


50 .55 


49 .25 






51 .70 


49 .25 


53 .40 





nnd the best valae for the angle. 116'' 43' 49".64 ± 0".28 
48. The following observations were made with a sextant 
in order to determine the latitude of a place : — 



48" 4' 46" 


7" 


59" 


52" 


47" 


36" 


24 


28 


89 


52 


15 


40 



What is the most probable latitude ? 43'' 4' 27" ± b" 

44. Determine the quantity M from the observations 

M = 81.55 .41 .68 .25 .27 .77 .13 .86 .10 .03 
Wts. 25 25 16 9 9 9 4 1 1 1 

45. An angle M is measured with the following results : 

3f p Jf" p J£ p Jf" p If p 

45".00 5 42".50 5 27".50 3 36".25 2 45".00 2 

31 .25 4 37 .50 3 43 .33 3 42 .50 3 40 .83 3 

45 .00 3 38 .33 3 40 .63 4 39 .17 3 

Find the most probable value of the angle and test the 
result by computing Sjtw. 39".78 ± 0".94 

46. In one hundred measurements of angles made in the 
primary triangulation in Massachusetts by the Coast Survey 
there were found between 

6".0 and 5".0 1 error. 1".0 and 0".0 26 errors. 

5 .0 « 4 .0 2 errors. .0 « — 1 .0 26 « 

4 .0 « 3 .0 2 « - 1 .0 « - 2 .0 17 « 

3 .0 " 2 .0 3 « - 2 .0 « - 3 .0 8 « 

2 .0 « 1 .0 13 « — 3 .0 « — 4 .0 2 ** 

Plot these results as in Figure 1, page 8. 



EXAMPLES. 128 

47. In sixty-six determinations of the velooity of light 
made at Washington the percentage of errors of different 
magnitudes was found to be as follows : — 



Over — 7.0 


.8% 


Over 


7.0 


0.0% 


Equal to — 7.0 


.8 


Equal to 


7.0 


1.6 


-6.0 


1.6 




6.0 


3.5 


-5.0 


2.0 




6.0 


4.2 


-4.0 


3.9 




4.0 


5.8 


— 3.0 


7.8 




3.0 


8.6 


-2.0 


12.8 




2.0 


12.4 


— 1.0 


17.1 




1.0 


17.1 



Plot the results and draw a smooth curve. 

48. Given the observations 

V 

I 

aj _ y + 22 = 3 405-1- y -f 42 = 21 

3a; -f- 2y — 52 = 6 — aj -f- 3y -j- ^^ = 14 

find the most probable values of x^ y^ and z. 

X = 2.470 ± .088 

49. being the level of the sea and Pi, P^ Pg, three 
points whose altitudes are to be determined, the following 
observations are made : — 



Pi above — 10 ft. 


P, above P. = 9 ft. 


Pa « Pi — 7 


Pi « P. = 2 


P, « = 18 





Find the most probable altitudes. P, = 8.50 ± .29 

50. The altitudes of A above Oj JB above Ay and JB 
above are found by measurements to be respectively 
12.3, 14.1, and 27.0 feet. . What is the most probable value 
of each of these differences in level? A = 12.50 ± .17 

51. Measurements of the ordinates of points on a straight 
line corresponding to abscissas 4, 6, 8, 9, are made with 
results 5, 8, 10, 12. What is the most probable equation of 
the line in the form y = mx -f- 6 ? 5 = — 0.29 
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52. Find the altitades in Example 49 if the observations 
have weights 5, 3, 6, 2, 4, respectively. 

53. Solve the example in paragraph 31, giving the obser- 
vations the weights 25, 25, 4, 4, 4, 4, 4, 4, 1, respectively. 

Elevation of Pj = 320.25 

54. Find the most probable values of 2i, s^, and z^ from 
the observations 

Zi = 552.10 wt. 16 «! — 2j = .75 wt. 1 

— J5, -|- 2j, = .15 « 9 2Ji + 2J, — 2J, = 552.05 « 1 

2, = 551.23 "4 J5i — 2, = .70 « 1 

2, = 551.30 « 4 

2, = 551.2345 

55. In the triangolation of Lake Superior there were 
measured at station the angles 

F P = 62^ 59' 40".33 wt. 5 

F B = 64 11 34 .92 « 7 

i?' jB = 100 20 29 .12 « 4 

P B = 37 20 49 .55 « 7 

ROB = 36 8 55 .86 « 4 

Required the adjusted values of the angles. 

F OP = 40''.28 ± 0".34 

56. In the U. S. Lake Survey the following angles were 
measured at station North Base : — 

(1) Crebassa— Middle 55** 57' 58".68 wt. 3 

(2) Middle — Quaquaming 48 49 13 .64 '' 19 

(3) Crebassa — Quaquaming 104 47 12 .66 « 17 

(4) Quaquaming — South Base 54 38 15 .53 <' 13 

(5) Middle — South Base 103 27 28.99 " 6 

Find the adjusted values of the angles. 

(1) = 58".965; r = 0".28 
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57. Adjust the following observations of differences in 
level : — 

Altitude of A 401.3 wt. 16 C above B 72.5 wt. 9 

A above B 220.8 " 16 J. « jB 222.0 « 1 

J. « (7 150.2 « 4 Altitude of jB 180.7 « 1 

58. In *^ Conditioned Observations" can the number of 
observations required be less than the number of unknown 
quantities? Why must the number of conditions be less 
than the number of unknowns ? 

59. From the following measurements of the angles formed 
at the centre of a disk by four radial lines, find the most 
probable values of the angles. 

A = 104° 25' 13" C = %&" 33' 20" 

jB = 98 13 47 2> = 70 48 23 

A = 104° 25' 2".25 

Also solve giving the observations the weights 5, 2, 1, 4, 
respectively. 

60. Four observations on the angle ^ of a triangle gave 
a meap of 36° 25' 47", two observations on B gave a 
mean of 90° 36' 28", and three on C gave 52° 57' 57". 
Adjust the triangle. -4 = 36° 25' 44".2 ; r = 7".7 

61. Five angles at a station are measured, and also their 
sum. The observed sum differs from the sum of the five 
observed parts by the amount d. What are the adjusted 
values of the angles ? 

62. The three angles of a spherical triangle are measured 
with results 

A = 46° 17' 38".32 B == 73° 35' 16".15 

(7 = 60° 7' 5".16. 

Adjust the triangle, knowing that the spherical excess is 
2".475. A = 39".3; /* = 1".6 
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63. At the Btation Pine Mountain the following angles 
were observed between surrounding stations : — 

Jocelyne — Deepwater 66** 11' 52".600 wt. 8 

Deep water — Deakyne 66 24 15 .553 «« 3 

Deakyne — Burden 87 2 24 .703 « 8 

Burden — Jocelyne 141 21 21.757 " 1 

Find the most probable values of the angles. 

64. Solve Examples 55 and 56 by the method of '* Con- 
ditioned Observations." 

65. ^ is a station whose altitude is known to be 5240.1 
feet. £ and G are floats on a lake, and 2> is a signal point. 
From the following observations determine the most prob- 
able altitudes of jB, G and D, 

G below A 720.1 wt. 3 JB below A 719.7 wt 3 
i> « ^ 200.3 « 5 JB '' D 520.9 « 2 
G '' I) 520.4 « 2 

66. Given the following observations, subject to the con- 
dition 2i -{~ ^9 = ^89 fii^d the most probable values of 
^1, ^3, and 2(. 

22i — 2, + 2, = 3.0 22, — 2, = 1.0 

22i — 32, = — 4.5 2i + 22, = 5.1 

2i + 2, = 3.8 

67. The chemical composition of a specimen was found 
by several observers to be as follows : — 

JPb = .52 Other substances = .09 Au and Ag = .39 
Ag = .27 JPb and Ag = .78 Impurities = .10 

Au = .11 , P6 and impurities = .62 Au = .12 

From these observations find the most probable composition 
of the specimen. 
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68. From the following observationB what are the best 
values of the unknowns, supposing that y and z must be 
equal ? 



a; + y — 5.2 


Wt. 4 


y + 2 4.2 


wt. I 


X — 3.0 


« 9 


z — 2.0 


" 4 


35 — 2 1.1 


« 1 







23» 


41».041 


wt.80 


42 


14 .875 


u 7 


47 


27 .718 


" 8 


23 


46 .816 


« 7 


5 


12 .929 


« 5 



69. In determining the difference in longitude between 
various cities the results obtained were 

(1 ) Cambridge — Washington 

(2) Cambridge — Cleveland 

(3) Cambridge — Columbus 

(4) Washington — Columbus 
I (5) Cleveland — Columbus 

Adjust these observations. 

70. The capacity of a condenser is known to be 14.0 m. f. 
It is divided into five sections, a, 6, c^ d^ 6, and it is known 
that the difference between b and d is 1.5 m. f. Find the 
most probable capacities of the sections from the observa- 
tions 

a = 2.02 wt. 3 d = 2.67 wt 7 

b = 4.13 « 2 e = 2.84 « 4 

i ^ = 2.52 « 5 ^ ^ 1 9651 

71. If the unknowns in the following observations are 
subject to the condition x -\- 2y -\- Sz = 36, what are 
their adjusted values ? 

X = 4.3 wt. 1, y = 5.7 wt. 4, 2 = 7.3 wt 9 

X = 3.77 

72. A cannon is discharged horizontally from the top of 
a bluff. Observations on the time, and distance of fall of the 
ball gave the results 

t = 0.5 1.0 1.5 2.0 seconds 
« = 1.2 4.0 9.1 15.0 metres 
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What curve, passing through the point of departure of the 
ball, will represent the above observations ? 

73. An Argand burner shows the following efficiencies 
with varying rates of gas consumption : 

g = 2.0 2.3 2.8 3.3 4.0 4.5 5.0 feet 
B = 2.1 2.4 2.5 3.0 3.2 3.8 4.1 

Find the equation of the straight line which best rep- 
resents the relation between g and JS. The measurements 
on g are without appreciable error. 

74. Observations are made upon the expansion of Amyl 
alcohol with change in temperature as follows : — 

V = 1.04 1.12 1.19 1.24 1.27 cu. cm. 
t = 13.9 43.0 67.8 89.0 99.2 C. degrees 

If V =z 1 -\- JB t -\- C t^ expresses the law connect- 
ing the volume and temperature, find the most probable 
values of JB and C. 

75. In a Hooke's joint where the angle between the axes 
is 45°, X being the angular rotation of the driver, and y 
that of the follower, from the following measurements find 
the equation of a curve that will represent the relation 
between x and y — • as. 



X 


y — X 


X 


y — X 


X 


y — x 


0** 


o°.o 


80° 


- 5°.8 


140° 


8°.8 


20*> 


- 5°.7 


90° 


- 2°.0 


160° 


5°.3 


40° 


— 9°.9 


100° 


2°.3 


180° 


o^o 


60° 


- 10°.4 


120° 


8°.0 







y — JB = — 0.85 — 9.82 sin 2x -f 0.92 cos 2x 

76. A series of observations extending over a period of 
thirty years was made by Quetelet to determine the daily 
variation in temperature at Brussels. The mean results of 
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the measarements are given below. From them derive an 
equation to express the temperature at any time of the year. 



Jan. 4°. 66 


May 9°. 83 


Sept. 8°.16 


Feb. 5^42 


June 10^.09 


Oct. 6^.56 


Mar. 6M7 


July 9°.71 


Nov. 5M0 


Apr. 8^.59 


Aug. 9M4 


Dec. 4^41 


y = 


7.369 4- 0.9854 sin 


30aj - 2.7084 cos 30a? 




4- 0.0100 sin 


60aj 0.1950 cos 60a5 




0.0133 sin 


90aj+ 0.1783 cos 90a; 



In this answer the values of x begin at the 15th of Janu- 
ary, and represent the time in months. 

77. The law connecting the time of vibration of a pendu- 
lum with its length is assumed to be of the form, T ■=. m X**. 
From the following observations find the most probable values 
of m and n. 

T = 12.9 11.6 10.4 9.7 5.3 4.6 
Z = 164.4 132.9 107.6 93.5 28.4 20.6 

X is in centimetres, T in tenths seconds. n z= 0.5000 

m = 1.0044 

78. Determine the equation of a curve which will repre- 
sent the following observations : 



X 


y 


X 


y 


X 


y 


0.0 


0.00 


1.5 


1.09 


3.0 


8.65 


0.5 


0.04 


2.0 


2.56 


3.5 


13.72 


1.0 


0.31 


2.5 


4.99 


4.0 


20.47 



79. Determine the equation of a curve which will repre- 
sent the relation between x and y in the observations. 



X 


y 


X 


y 


X 


y 


X 


y 


0.0 


4.51 


0.3 


4.09 


0.6 


3.03 


1.2 


0.92 


0.1 


4.44 


0.4 


3.76 


0.8 


2.24 


1.5 


0.38 


0.2 


4.31 


0.5 


3.42 


1.0 


1.49 


2.0 


0.05 
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80. At a station P the angles between a straight line 
passing through P parallel to the axis of JT and the direc- 
tions from P of four points Pi, P^ Pj, P4, are measured. 
Having given the coordinates, (a, d), of the four points, 
find the coordinates of P. 

) 
Point. Oodrdinates. Angle. 





a 







Pl 


4.21 


3.24 


89» 18' 


p» 


1.21 


2.10 


147° 54' 


p» 


— 0.51 


0.22 


205» 24' 


p. 


2.50 


- 1.10 


277'' 15' 



If the coordinates of the point P are (a?, y), and the 
angle is denoted by A^ we have 

y — ft 

tan.^1 = 

X — a 

81. li asinhx = My and values of Jtf are observed 
for known values of a and ft, determine the most probable 
value of X. If x' is an approximate value of x found by 
trial, and m =z a ain bx' — M^ we shall have 

^ 2 a ft m cos bx' 
* ~ * "" S (a ft cos ft*')* 

82. If in one series of observations the value of h is 
twice what it is in another, what is the relative probability of 
the occurrence of an error of given magnitude a in the two 
series ? Show what the curves of error will be in the two 
cases. What error has the same probability for its occurrence 
in each series ? What is the relative probability of the occur- 
rence of an error not greater than a in the first case and not 
greater than 2a in the second case ? 

83. From 64 observations the latitude of a station was 
found to be 49** 10' 9'MIO ± 0".051. What was the prob- 
able error of a single observation ? 0".41 
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84. If twenty measurements of an angle give a result with 
an A,D, of 0".38, and it is required to find the angle so 
that the A.D. shall be only 0".25, how many more observa- 
tions must be made ? 27 

85. From the following determinations of the area of a 
field find the most probable area and its probable error. 

5674 ± 12, 5680 ± 4, 5685 ± 3, 5682 ± 1, 5678 ± 2 

A = 5681.41 ± 0.84 

86. From the following measurements by Fizeau and 
others, find the most probable value for the velocity of light 
together with its probable error. Measurements are in kilo- 
meters. 

298000 ± 1000 299990 ± 200 299930 ± 100 
298500 ± 1000 300100 ± 1000 

V = 299917 ± 88 

87. Two different instruments give for the value of an 
angle, 

34*^ 55' 33".0 ± 4M, 34** 55' 36".0 ± 6".3 

What is the best value to take for the angle ? 

34<* 55' 33".9 ± 3".4 

88. Determinations of the difference in longitude between 
Washington and Key West made on seven different days 
gave the results 

19"* 1'.42 ± 0'.044 19*" 1'.60 ± 0'.046 

1 .37 ± .037 1 .55 ± .045 

1 .38 ± .036 1 .57 ± .047 
1 .45 ± .036 

What is the best value and its probable error ? 

1'.460 ± 0'.016 
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89. In the triangulation of Lake Ontario two different 
instruments gave for an angle, 74** 25' 5".429 ± 0".29 from 
sixteen readings, and 74*' 25' 4".611 ± 0".22 from twenty- 
four readings. Find the most probable value of the angle 
and its probable error. 

90. In each of Examples 39-45 find the mean and prob- 
able errors and average deviation of each observation and of 
the most probable value, using formulas from (50) to (63) 
according as they apply. 

91. In Example 42 divide the observations in their order 
into six groups of four observations each and compute the 
mean of each group. Then determine the probable error of 
the first of these means : ( 1 ) considered as a single measure 
of four times the weight of those in Example 42 ; (2) directly 
as one of six observations of equal weight; (3) as a deter- 
mination from its four constituents. 0".67 ; 0".72 ; 1".00 

92. The following twenty-nine measurements on the den- 
sity of the earth, made by Cavendish, give as a mean result 
5.48. What is the probable error of an observation ? Solve 
by the usual method and also by taking the residual that 
occupies the middle position. 0.14 

5.50 5.55 5.57 5.34 5.42 5.30 



.61 


.36 


.53 


.79 


.47 


.75 


,88 


.29 


.62 


.10 


.63 


.68 


.07 


.68 


.29 


.27 


.34 


.85 


.26 


.65 


.44 


.39 


.46 





93. What is the probable error of the mean of two obser- 
vations which differ by the amount a ? 

94. A base-line is measured five times with a steel tape 
reading to hundredths of a foot, and five times with a chain 
reading to tenths of a foot, with results 

By tape, 741.17 741.09 741.22 741.12 741.10 
By chain, 741.2 741.4 741.0 741.3 741.1 
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Find the probable errors and weights for a single observa- 
tion in each case, and also the adjusted length of the line 
and its probable error. 741.146 ± 0.015 

95. Twenty-one determinations of a chronometer correc- 
tion gave results 

- 8.78 - 8.78 - 8.68 - 8.80 - 8.96 - 8.83 - 8.79 
.76 .51 .63 .75 .64 .70 .90 

.85 .64 .58 .78 .65 .64 .93 

Find the probable error of the mean by using both formulas 
(53) and (57), and also determine the probable error of a 
single observation by taking the middle residual. 

0.017; 0.018; 0.09 

96. In the following observations show that Mq = 49.64, 
;* = 1.95, r = 1.31, fjLo = 0.40, Tq = 0.27, /i, = 0.87, 
r, = 0.59. 

M = 48.81 48.76 49.53 51.56 50.88 49.84 
p = b 4 5 3 2 5 

97. Observations on the time of ending of a transit of 
Mercury are made by different observers with a variety of 
instruments and under more or less favorable circumstances. 
If the weights assigned by the computer are as indicated, find 
the best value for the time and its probable error. 

6* 38"* 23* wt. 1 38*" 26' wt 3 38*" 19* wt 8 

37 55 « 38 21 « 2 38 21 " 2 

38 10 « 1 38 18 " 2 38 15 « 2 

to = 5* 38*" 19'.9 

98. An angle is measured five times with a theodolite, and 
seven times with a transit, giving results 

Theodolite, 31".7, 39".8, 40".7, 28".6, 32".3 

Transit, 32 .8, 36 .7, 38 .2, 29 .3, 41 .6 35".3, 36".2 
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If the relative valaes of readings by the two instroments 
are as 3 to 2, what is the most probable valae of the angle ? 
What is the mean error of the result ? 

99. Given M^ = 65.58 ± .59, M^ = 35.15 ± .93, 
M^ = 49.64 ± .27, find the probable errors of 4 Jfi — 

3 JM; + 2 JKi and of ^ + ^ ^ ^. 3.69; 0.43 

100. The three angles of a triangle are measured, and the 
probable error of each observation is r . What is the prob- 
able error of the triangle error ? r ^T" 

101. The zenith distance of a star on the meridian is 
observed to be 2 = 21^ 17' 20".3 ± 2".3. The declina- 
tion of the star is given as rf = 19** 30' 14".8 ± 0".8. 
What is the latitude of the place and its probable error ? 

Z = z -\- d =z 40** 47' 35".l ± 2"4. 

102. The zenith distance a; of a star at upper culmina- 
tion is observed n times, and its zenith distance z' at 
lower culmination n' times. If the latitude is given by 
X = 90° — i (z -\- z')y and the probable error of an 
observation is r, what is the probable error of the latitude ? 

103. The horizontal force necessary to start a 100-pound 
weight sliding along a table is observed to be 15.5 ± 0.2 
pounds. Find the probable error of the coefficient of friction. 

104. If a line is measured by the continued application of 
a unit of measure, and r is the probable error of the placing 
and reading of this measure, what is the probable error of 
the length / ? ^ Jf 

105. If the average deviations of 2^, z^ 2,, are a, 6, c, 
respectively, what is the average deviation of Zi^ +25,^+28*? 

106. If the radius of a circle is measured with result 
1000.0 ± 2.0, how should the circumference and area be 
expressed ? 

107. Two sides, a and by and the included angle C of 
a triangle are measured with results a = 252.52 ± .06 
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feet, h = 300.01 ± .06 feet, C = 42^ 13' 00" ± 30". 
What is the area and its probable error ? 25452 ± 9 

108. Measarements of adjacent sides of a rectangle gave 
a ± rj, and b ± r,. What is the probable error of the 
area, and for what kind of a rectangle will this probable 
error be the least ? 

109. If the measured sides of a rectangle have the same 
a.d.^ what is the a,d. of the diagonal determined from 
them ? Same 

110. If the sides of a rectangle are measured in the 
manner indicated in Example 104 and found to be a and d, 
what is the probable error of the area ? 

111. The correction to be applied to a chronometer is 
found to be -|- 12"* 13'.2 ± 0*.3. Ten days later the cor- 
rection is again determined and found to be 12"* 21^.4 ± 0'.3. 
What is the mean daily rate and its probable error ? 

0'.820 ± 0'.042 

112. Measurements of the compression of the earth's 
meridian have resulted in 

^ ± .000046 



294 



What is the probable error of the denominator 294 ? 3.98 

113. The current flowing in a circuit is due to two 
sources whose electromotive forces are determined to be 
«! = 200 ± 2, e, = 400 ± 3. The resistance of the 
circuit is 30 ± 1. Find the current and its probable error. 

20 ± 0.68 

114. The side b and angles B and C of a triangle are 
measured with results b = 106 ± .06 metres, jB = 
29° 39' ±1', C — 120^ T ± 2'. What is the most 
probable value of the angle A and of the side c ? 

A = 30° 14' ± 2'.2; c = 185.35 ± 0.15 

115. The distance between two divisions on a graduated 
scale is measured by a micrometer. Show that the average 
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deviation of the mean of two results is the same as the aver- 
age deviation of a single reading. 

116. If the weights of the determinations of three angles 
Ay JB, Cy are 3, 3, 1, respectively, what is the weight of 
the sum of the three angles ? 0.6 

117. If the weight of x is />, what is the weight of 
loga X ? 

118. If 35 = — and the weight of y is />, what is the 
weight of 85 ? c^p 

119. In Example 107, how closely mast the parts be 
measured in order to obtain the area within 0.5 per cent ? 

120. From observations on / and t the value of ^ is 
to be computed by the pendulum formula 



t 



V g 



What changes in g will be produced by changes in / and 
^ of Si and 82 units, respectively, and what are the allow- 
able errors in / and ^ if ^ is to be determined within 
1 per cent ? 

121. The moment of inertia of a cylindrical bar is to be 
obtained from measurements on its mass m, its length A, 
and its diameter d. The error in the determination of m 
is negligible, the precision of the determination of d is four 
times that of h. If the measurements give m = 48, 
h = 8.000, d = 1.200 ± 0.10, and 



\l2 ^ 16/ 



what is the probable error of J^ and what should be the 
ratio of d to A to determine jT most accurately ? 

d: h = 256 : 9 



\ 
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122. If observations give for a certain quantity x the 
value 303, with a mean error of 2, what is the mean error 
of the expression 3 as -f- ^^gio 2 a;? 

123. The probable error of the determination of the 
angle A is 20". What is the maximum probable error of 
sin A -f- cos A ? 

124. If the probable error of an observation on an angle 
is 10", is there any difference between the probable error 
of the function sin A -\- cos -4 -f- sin C and of the func- 
tion sin A + cos B -j- ^^^ ^^ supposing A and B are 
of the same magnitude ? 

125. Given the observations, 

21-222+ 23-3 = 32i+ 22 + 228-17 = 

22i + 322 - 42s - 2 = - 2l + 42a + 328-10 = 

Find the most probable values of 2^, 23, 29, and also their 
weights and precision measures. 

2i = 3.541; p^ = 29; r,, = .024 

126. Find the weights and precision measures of the 
unknowns in Examples 48 to 57. 

127. Determine the probable errors of the constants in 
Examples 72 to 79, inclusive. 

128. The length of a pendulum which beats seconds is 
given by 




where I' is the length at the equator, q the ratio of 

^o9 

the centrifugal force at the equator to the weight, and 8 the 

compression of the meridian regarded as unknown. Putting 

r = 991'""*- + a, l—q—8]rz=y. 




observations in different latitudes gave in millimetres the 
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foliowii^ equatiomg, from which we are to determine / and 
8 together with their probable errors : — 

X + 0.969y = 5.13 X 4- 0.152y = ^0.77 

X 4- 0.749y = 3.97 x -f 0.327y = 1.70 

X 4- 0.426y = 2.24 x -f 0.685y = 3,62 

X 4- 0.095y = 0.56 x \~ 0.793y = 4.23 

X = 0.19 

V = 991.669 ± .02^; s = — ± 0.00046 

294 

129. Find the weights and precision measures of the 
unknowns in Examples 64 to 70, and also in Examples 59 to 
63, and in 71. 

130. In Example 95 the probable error of a single obser- 
vation is 0.08 seconds. Find the number of errors which 
should fall between 0.^0 seconds and 0.10 seconds, between 
0.10 seconds and 0.20 seconds, and also the number tbat 
should be over 0.20 seconds. CJompare the results with the 
number flctoally found. 

131. In 470 determinatione of Xkte rv^X asceiMionB of 
Sirius and Altair made by Bradley, the probable error of a 
single observation was 0".2637. Tlie number of errors falling 
between specified lindtB was as -siiown below. Compare this 
result with the distribution of errors called for by the theory. 



Limits. 


Errors, 


Limits. 


Errors. 


0'\0 to 0".l 


94 


0".6 to 0".7 


26 


.1 to .2 


88 


.7 to .8 


14 


.2 to .3 


78 


.8 to .9 


10 


.3 to .4 


58 


.9 to 1 .0 


7 


.4 to .5 


51 


Over 1 .0 


8 


.5 to .6 


36 







132. What is the probability that the error of a single 
observation will be as large as twice the probable error? 
As large as five times the pro^bable error ? 
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133. On the average how many observations must be made 
before an error as large as three times the mean error will 
occur ? 

134. In Example 46, assuming that all errors between any 
two limits fall half way between those limits, compute the 
average deviation and mean error of an observation and com- 
pare their ratio with the theoretical value given in the table 
in paragraph 55. 

135. A line is measured 500 times and the probable error 
of each observation is 0.6 cm. How many errors should 
occur between 0.4 cm. and 0.8 cm. ? 

136. Show how the value of v could be determined 
experimentally from observations such as those in Example 
131. 

137. In a system of observations all equally good, r being 
the probable error of a single observation, if two observations 
are taken at random, what quantity is their difference as 
likely as not to exceed, and what is the probability that the 
difference will be less than r ? |. ^^; 0.367 

138. In the following measurements of an angle, ought 
any t)f the observations to be rejected ? 

12' 51".75 47".85 47".40 48".90 44".46 
48 .45 51 .05 48 .85 50 .95 
50 .60 47 .75 49 .20 50 .55 

139. Determine whether any of the observations in 
Example 44 should be rejected. 

140. A quantity M is measured with the results given 
below. Ought all the observations to be retained ? 

M = 236, 251, 249, 252, 248, 254, 246, 257, 243, 274 

141. A certain angle has been laid out with such accu- 
racy that its true value may be taken as exactiy 90°. Twenty- 
five observations are made upon it with a transit that it is 
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desired to test, and the result obtained is Sd"* 59' 57" ± 0".8. 
What are the odds in favor of a constant error in the instru- 
ment between — 1" and - 5"? Between 0" and — 6"? 

908:92; 86:1 

142. Repeated measurements of a standard metre bar 

with a decimetre scale gave a result 10.032 ± 0.010. What 

are the odds in favor of a constant error in the scale between 

0.010 , 0.054 « .„ - 
and ? 43 : 7 

Vio vn[o 

143. Two determinations of the length of a line gave 
683.4 ± 0.3 and 684.9 ± 0.3, respectively. Show that the 
best value for the length is 684.15 ± 0.51, and that the 
probable systematic error of each determination is 0.65. 

144. Two men A and ^ observe an angle repeatedly 
with the same instrument with results 

A. B. 

47° 23' 40" 23' 35" 47° 23' 30" 24' 00" 

23 45 23 40 23 40 23 20 

23 30 23 50 

Is there any relative personal error, and what is the best 
final value? 47° 23' 38".2 ± 1".6 

145. Three independent determinations of the capacity 
of a condenser made with three different instruments gave 
results 42.22 ± 0.21, 43.40 ± .15 and 44.20 ± 0.18. 
What is the most probable value of the capacity ? 

For extended treatment of the subject illustrated in Exam- 
ples 143 to 145 see Johnson^ "The Theory of Errors and 
Method of Least Squares," chap. vii. 

146. In an estimation of tenths what is the probable error 
of an observation ? What is the average deviation ? 0.025 

147. In obtaining the angle of deflection of the needle of 
a tangent galvanometer by the usual method what is the 
probable error of the result ? 
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148. If all the errors of a series of observations must fall 
between and a, and the frequency of any error is pro- 
portional to its magnitude, what is the Curve of Error? 

What are the values of r, a.<?., and /a ? r = 

149. In Example 148 what is the probability that the error 
of a single observation will be as large as 0.5 a. 

150. If all values of x between and a are possible, 
and their probabilities are proportional to their squares, find 
the mean value of x and the probability that x will be as 
large as 0.5 a. Also draw the Curve of Error. 

151. What is the greatest probable error of a logarithm 
found by interpolation in a seven place table ? .000000016 

1 52. Given the following set of Normal Equations, together 
with [mm] = 1.3409, find the most probable values of 
the unknowns and their weights and probable errors. There 
were sixteen observations. 

3.1217 2i + .5756 z^ - .1565 Zj, - .0067 z^ 

- 1.5710 = 
.5756 2i 4- 2.9375 z^ + .1103 2, — .0015 z^ 

+ .9275 = 

- .1565 2i + .1103 2a + 4.1273 a;, + .2051 z^ 

+ .0652 = 

- .0067 2i - .0015 2, -f .2051 «, + 4.1328 z^ 

4- .0178 = 

2i = 0.583 ± 0.018 

24 = -. 0.004 ± 0.016 

Pz, = 4.12 

153. From ten observation equations, for which was found 
[mm] = 2.6322, there resulted the normal equations 

5.2486 2i - 1.7472 2, - 2.1954 2, -f 0.6399 = 

- 1.7472 2i -f 1.8869 2, -f 0.8041 2, - 1.4493 = 

- 2.1964 2x -j- 0.8041 2, -f 4.0440 2, - 1.8681 = 
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Find the most probable values of Zi, z^ and z^ together with 
their probable errors. z^ = 0.42 ± 0.11 

154. Find the most probable values of the unknowns in 
the normal equations 

459 2i - 308 Zj -- 389 2, + 244 z^ - 507 = 

- 308 «i 4- 464 «, -f- 408 s, ~ 269 z^ + 695 = 

- 389 J5i -f 408 2, -f 676 s, - 331 z^ + ^53 = 
244 afi — 269 z^ — 331 2, + 469 z^ — 283 = 

[mm] = 1129 

24 = — 0.488 ; p^^ = 281 

155. If thirteen observation equations give rise to the 
result [mm] = 100.34 and to the normal equations 

17.50 «i — 6.50 2, — 6.50 2, — 2.14 = 

- 6.50 2i + 17.50 2a — 6.50 2, — 13.96 = 

- 6.50 2i - 6.50 2, + 20.50 2, + 5.40 = 

show that the most probable values of the unknowns are 
2i = 0.67 ± 0.60, 2, = 1.17 ± 0.60, 2, = 0.32 ± 0.55 
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BLBMBHTS OF TBB THEORY OF PBOBABZUTT. 

200. Definition. If an event can happen in a ways, and 
fail in b ways, and all these ways are equally likely to oconr, 

the probability of the happening of the event is — ir~X> 

and the probability of its failure is — TTX* 

Since the event must either happen or fail, the sum of the 
above probabilities must represent a certainty. But 



+ -:rnr-T = i 



a ^ b ^ a -i- b 

That is, the probability of a certainty is expressed by unity. 
Also, if the probability, P, of the happening of an event is 
known, the probability of its failure is given at once by 

201. Example A, A single throw is made with a pair of 
dice. What is the probability that the sum of the spots 
turned up will be 5 ? 

Number of ways of throwing the dice is 6 X ^ = 36 
Number of ways of throwing five is 4 

4 1 

.'. Probability of throwing five is — = — 

Example B. A coin is tossed up six times. Find the 
chance that three heads and three tails will be the result. 



ii METHOD OF LEAST SQUARES. 

Number of ways of throwing the coin is 2* = 64 

Number of ways of throwing three heads is = 20 

1X2X3 

20 5 

Probability of throwing three heads is — = — 

64 16 

202. Compound Events. A certain event can happen in 
a ways and fail in b ways : a second independent event can 
happen in a' ways, and fail in b' wskja^ all of these ways 
being equally likely to occur. 

To find the probability of the simultaneous occurrence of 
the two events. 

The total number of ways in which the events can take 
place together is (a + ^) («' + ^') 

(1) Both events can happen in a a' ways. 

(2) Both events can fail in b b' ways. 

(3) First event can happen and second fail in a b' ways. 

(4) First event can fail and second happen in a' b ways. 



The probability of (1) is 



a a' 




(a + b) (a' 


+ *') 


bb' 




(a + b) (a' 


+ *') 


ab' 




(a + *) (a' 


+ «') 


a' b 





The probability of (2) is 



The probability of (3) is 



The probability of (4) is , , ,. , . , ,,^ 

(« + ^) (« + * ) 

But the probability of the happening of the first . event is 

ft Oi 

-. — 5^, and of the second event is -; — 5 — =-/> etc. Hence it 

a -|- " a + b' 

will at once be seen that the probability of the simultaneous 
occurrence of two independent events is equal to the product 
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of the probabilities of the occurrence of the component events. 

Or, in general, if Pi, Pjj • • • -^n are the probabilities of 

the occurrence of any number, n, of independent events, the 

probability of the simultaneous occurrence of all the events is 

Pi X A X ... Pn (A) 

By independent events is meant those such that the manner 
of occurrence of one has no influence upon the manner of 
occurrence of the others. 

203. Example C. The chance that A can solve a cer- 

2 
tain problem is — , and the chance that J3 can solve it is 

—. Find, 
12 

(a) The probability that both will solve it. 

(b) The probability that the problem will be solved. 

For (a). This is a question as to the probability of the 
concurrence of two independent events. Therefore by an 
application of (A), the probability that both will solve the 
problem is 

2 A _ A 

3 ^ 12 ~ 18 

For (b). The problem will be solved unless both fail. 

17 7 

The probability that both will fail is — v — = — 

^ ^ 3 ^ 12 36 

7 29 

The probability of getting a solution is 1 = — 

36 36 

Example D. A pack of cards is cut, and those taken off 
then replaced. In how many trials will it be an even wager 
that an ace will be cut ? 
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Let n be the number of trials. Then n is to be found 
from 



\62/ 2 



where the first member of the equation represents the proba- 
bility that we shall not fail n times in succession. 

Solving for n, 

log 2 
n = 5 

log 52 ~ log 48 

= 8.7 

In nine trials then there is a little more than an even chance 
of cutting an ace. 

204. Dependent Events. If we have a number of events 
whose modes of occurrence are dependent one upon another, 
the probability of their concurrence will be found by the same 
method as in paragraph 202 ; a' now denoting the number 
of ways in which after the first event has happened the 
second will follow, and h' the number of ways in which after 
the first has happened the second will not follow, etc. Accord- 
ingly, the general formula (A) of paragraph 202 applies to 
dependent events as well as to independent ones. 

Also, if an event can take place in a variety of ways, the 
total probability of its occurrence will be the sum of the 
probabilities of its occurrence in each of the different ways. 

205. JSxample E. Suppose two purses contain respect- 
ively five dimes and a copper, and six dimes. A coin is taken 
at random from the first purse and placed in the second, and 
then a coin is transferred from the second to the first. What 
is the probability that the copper will remain in the first purse ? 

The probability that the copper will be taken from the first 
purse and placed in the second, and then returned to the first 
purse is 

1 1 _ i_ 
6^7 ""42 
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and the probability that the copper will not be taken from 
the first parse at all is 

5 
"6 

Therefore the probability that the copper will finally remain 
in the first purse is 

J_ 5^ _ 36 _ 6^ 
42 "^ 6 ""42 ""7 



FUNCTIONS OF SEVERAL VARIABLES. 

206. For the application of Taylor's Theorem to the 
expansion of a function of several independent variables, see 
Osborne^ 8 << Differential and Integral Calcalus," page 145. 
And for the conditions that lead to maxima and minima 
values of such functions, see page 165 of the same work. 
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For a complete list of works on the Method of Least 
Squares published up to 1876, see 

Merrimarij ** A List of Writings relating to the Method of Least 
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1877. 
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TABLES. 



vu 



Values of the Integral 



TABLE I. 



Vtt 



jy 



^dt for Argnment 



t a 

or — 

0.4769 r 



r 



0.0 
0.1 
0.2 
0.3 
0.4 

0.5 
0.6 
0.7 
0.8 
0.9 

1.0 
1.1 
1.2 
1.3 
1.4 

1.6 
1.6 
1.7 
1.8 
1.9 

2.0 
2.1 
2.2 
2.3 
2.4 

2.5 
2.6 
2.7 
2.8 
2.9 

8.0 
3.1 
3.2 
3.3 
3.4 

8. 

4. 

6. 



r 



0.0000 0.0054 0.0108 0.0161 0.0216 

0538 0691 0645 0699 0762 

1073 1126 1180 1233 1286 

1603 1656 1709 1761 1814 

2127 2179 2230 2282 2334 

0.2641 0.2691 0.2742 0.2793 0.2843 

3143 3192 3242 3291 3340 

3632 3680 3728 3775 3823 

4105 4162 4198 4244 4290 

4662 4606 4661 4696 4739 

0.6000 0.6043 0.6085 0.5128 0.5170 

6419 6460 6600 5640 6681 

6817 6866 6894 5932 6970 

6194 6231 6267 6303 6339 

6650 6684 6618 6662 6686 

0.6883 0.6915 0.6947 0.6979 0.7011 

7196 7226 7265 7284 7313 

7486 7612 7640 7667 7694 

7763 7778 7804 7829 7854 

8000 8023 8047 8070 8093 

0.8227 0.8248 0.8270 0.8291 0.8312 
8433 8463 8473 8492 8611 
8622 8639 8657 8674 8692 
8792 8808 8824 8840 8866 
8946 8960 8974 8988 9002 

0.9082 0.9095 0.9106 0.9121 0.9133 
9205 9217 9228 9239 9250 
9314 9324 9334 9344 9364 
9410 9419 9428 9437 9446 
9496 9603 9511 9619 9526 



6 



8 



9 



0.9570 
9635 
9691 
9740 
9782 

0.9570 
9930 
9993 

1.0000 



0.9577 
9641 
9696 
9744 
9786 

0.9635 

894S 

wra 



0.9583 

9647 
9701 
9749 
9789 

0.9691 

9954 

QQQA 



0.9690 
9662 
9706 
9763 
9793 

0.9740 

9963 

9997 



0.9697 
9668 
9711 
9767 
9797 

0.9782 

9970 

9997 



0.0269 0.0323 0.0377 0.0430 0.0484 

0806 0660 0913 0966 1020 

1339 1392 1445 1498 1561 

1866 1918 19n 2023 2076 

2386 2436 2488 2689 2690 

0.2893 0.2944 0.2994 0.3043 0.3093 

3389 3438 ^487 3636 3683 

3870 3918 8966 4012 4069 

4336 4381 4427 4472 4617 

4783 4827 4860 4914 4967 

0.6212 0.5264 0.5295 0.5337 0.5378 

6620 6660 5700 6739 6778 

6008 6046 6083 6120 6167 

6376 6410 6446 6480 6615 

6719 6753 6786 6818 6851 

0.7042 0.7073 0.7104 0.7134 0.n65 

7342 7371 7400 7428 7457 

7621 7648 7676 7701 7727 

7879 7904 7928 7962 7976 

8116 8138 8161 8183 8205 

0.8332 0.8353 0.8373 0.8394 0.8414 

8530 8649 8667 8585 8604 

8709 8726 8742 8769 8775 

8870 8886 8901 8916 8930 

9016 9029 9043 9056 9069 

0.9146 0.9158 0.9170 0.9182 0.9193 

9261 9272 9283 9293 9304 

9364 9373 9383 9392 9401 

9464 9463 9471 9479 9487 

9534 9541 9548 9656 9563 

0.9603 0.9610 0.9616 0.9622 0.9629 

9664 9669 9675 9680 9686 

9716 9721 9726 9731 9735 

9761 9766 9770 9774 9778 

9800 9804 9807 9611 9614 

0.9818 0.9648 0.9874 0.9896 0.9915 

9976 9981 9986 9968 9990 

9998 9998 9999 9999 9999 



Diif. 



54 
54 
58 
52 
51 

50 
49 
46 
45 
43 

41 
89 
37 
85 
32 

30 
28 
26 
24 
22 

19 
18 
17 
15 
13 

12 

10 

9 

8 

7 

6 
5 

5 

4 
4 



6 



8 



9 



Difl. 
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TABLE II. — Common Logarithms. 



n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


Diil. 


10 


0000 


0048 


0066 


0128 


0170 


0212 


0263 


0294 


0384 


0374 


42 


11 


0414 


0463 


0492 


0531 


0669 


0607 


0646 


0682 


ono 


0765 


38 


12 


0792 


0628 


0864 


0699 


0934 


0969 


1004 


1038 


1072 


1106 


86 


13 


1139 


1178 


1206 


1239 


1271 


1308 


1335 


1367 


1399 


1430 


32 


14 


1461 


1482 


1628 


1663 


1684 


1614 


1644 


1673 


1703 


1732 


80 


15 


1761 


1790 


1818 


1847 


1875 


1903 


1931 


1960 


1987 


2014 


28 


16 


2041 


2068 


2096 


2122 


2148 


2176 


2201 


2227 


2253 


2279 


27 


17 


2304 


2330 


2365 


2880 


2405 


2430 


2455 


2480 


2604 


2629 


25 


18 


2663 


2577 


2601 


2626 


2648 


2672 


2696 


2718 


2742 


2765 


24 


19 


2788 


2810 


2833 


2856 


2878 


2900 


2928 


2945 


2967 


2989 


23 


20 


3010 


3032 


3064 


3076 


3096 


3118 


8139 


3160 


3181 


8201 


21 


21 


3222 


8243 


3263 


3284 


3304 


3324 


3345 


3366 


3385 


3404 


20 


22 


3424 


3444 


3464 


3483 


3502 


3522 


3641 


3660 


8570 


8698 


18 


23 


3617 


8636 


3665 


3674 


3692 


3711 


3729 


3747 


8766 


3784 


18 


24 


3802 


3820 


3838 


3856 


3874 


3892 


3909 


3927 


3946 


3962 


18 


26 


3979 


8997 


4014 


4031 


4048 


4065 


4062 


4099 


4116 


4133 


17 


26 


4160 


4166 


4183 


4200 


4216 


4232 


4249 


4265 


4281 


4296 


17 


27 


4314 


4330 


4346 


4362 


4378 


4393 


4400 


4425 


4440 


4466 


16 


28 


4472 


4487 


4602 


4518 


4533 


4648 


4664 


4679 


4694 


4609 


16 


20 


4624 


4639 


4654 


4669 


4683 


4696 


4713 


4728 


4742 


4767 


16 


30 


4771 


4786 


4800 


4814 


4829 


4848 


4867 


4871 


4886 


4900 


14 


31 


4914 


4928 


4942 


4965 


4969 


4983 


4997 


5011 


5024 


5038 


14 


32 


6061 


5065 


5079 


6092 


5105 


5119 


5132 


5146 


5160 


5172 


13 


88 


6186 


5198 


5211 


5224 


6237 


6260 


5263 


5276 


6288 


5302 


13 


34 


6316 


6328 


5340 


5363 


5366 


5378 


5391 


5403 


5416 


5428 


13 


35 


6441 


6453 


5466 


5478 


5490 


5602 


5514 


5527 


5638 


5651 


12 


36 


6563 


6676 


5687 


6699 


6611 


6623 


5636 


6647 


6658 


5670 


12 


37 


6682 


6694 


6705 


6717 


6729 


5740 


5762 


5763 


5776 


5786 


12 


38 


6798 


6809 


5821 


5832 


6843 


6866 


5866 


6877 


5888 


6899 


11 


39 


6911 


5922 


6933 


5944 


5965 


5966 


5977 


6088 


5999 


6010 


11 


40 


6021 


6031 


6042 


6063 


6064 


6075 


6085 


6096 


6107 


6117 


11 


41 


6128 


6138 


6149 


6160 


6170 


6180 


6191 


6201 


6212 


6222 


11 


42 


6232 


6243 


6263 


6263 


6274 


6284 


6294 


6304 


6314 


6325 


10 


48 


6336 


6346 


6365 


6365 


6376 


6885 


6395 


6405 


6416 


6426 


10 


44 


6435 


6444 


6454 


6464 


6474 


6484 


6493 


6603 


6613 


6522 


10 


45 


6532 


6642 


6561 


6661 


6671 


6680 


6690 


6699 


6609 


6618 


10 


46 


6628 


6637 


6646 


6666 


6665 


6675 


6684 


6693 


6702 


6712 


9 


47 


6721 


6730 


6739 


6749 


6768 


6767 


6776 


6785 


6794 


6803 


9 


48 


6812 


6821 


6830 


6839 


6848 


6857 


6866 


6875 


6884 


6893 


9 


49 


6902 


6911 


6920 


6928 


6937 


6946 


6965 


6964 


6972 


6981 


8 


60 


6990 


6998 


7007 


7016 


7024 


7033 


7048 


7050 


7060 


7067 


8 


61 


7076 


7064 


7093 


7101 


7110 


7118 


7126 


7136 


7143 


7152 


8 


62 


7160 


7168 


7177 


7186 


7193 


7202 


7210 


7218 


7226 


7235 


8 


63 


7243 


7261 


7259 


7267 


7276 


7284 


7292 


7300 


7308 


7316 


8 


64 


7324 


7332 


7340 


7348 


7366 


7364 


7372 


7380 


7388 


7396 


8 


n 





1 


2 


3 


4 


6 


6 


7 


8 


9 
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IX 



TABLE II. — Ck>mmon Logarithms. 



n 





1 


2 


3 


4 


6 


6 


7 


8 





Difl. 


65 


7404 


7412 


7410 


7427 


7485 


7443 


7461 


7460 


7466 


7474 


8 


66 


7482 


7490 


7497 


7506 


7513 


7620 


7528 


7536 


7648 


7661 




67 


7650 


7566 


7674 


7682 


7580 


7607 


7604 


7612 


7610 


7627 




68 


7634 


7642 


7649 


7657 


7664 


7672 


7679 


7686 


7604 


7701 




60 


7700 


7716 


7723 


7731 


7788 


7745 


7762 


7760 


7767 


7774 




60 


7782 


7780 


7796 


7808 


7810 


7818 


7826 


7882 


7830 


7846 


7 


61 


7868 


7860 


7868 


7876 


7882 


7889 


7896 


7003 


7010 


7817 




62 


7024 


7931 


7938 


7046 


7962 


7960 


7966 


7073 


7080 


7087 




63 


7003 


8000 


8007 


8014 


8021 


8028 


8036 


8041 


8048 


8065 




64 


8062 


8060 


8076 


8082 


8080 


8096 


8102 


8100 


8116 


8122 




65 


8120 


8136 


8142 


8140 


8156 


8162 


8168 


8176 


8182 


8180 


7 


66 


8196 


8202 


8200 


8216 


8222 


8228 


8235 


8241 


8248 


8264 




67 


8261 


8267 


8274 


8280 


8287 


8288 


8299 


8306 


8312 


8810 




68 


8825 


8831 


8338 


8344 


8361 


8367 


8363 


8370 


8376 


8382 




60 


8388 


8806 


8401 


8407 


8414 


8420 


8426 


8432 


8438 


8445 




70 


8461 


8457 


8468 


8470 


8476 


8482 




8404 


8600 


8006 


6 


n 


8613 


8510 


8625 


8531 


8637 


8548 


8649 


8566 


8561 


8667 




72 


8673 


8579 


8585 


8691 


8607 


8603 


8600 


8616 


8621 


8627 




T8 


8R33 


8630 


8646 


8651 


8657 


8663 


8600 


8675 


8681 


8686 




74 


8602 


8686 


8704 


8710 


8716 


8722 


8727 


8783 


8730 


8746 




76 


8761 


8756 


8762 


8768 


8774 


8770 


8785 


8701 


8707 


8808 


6 


76 


8808 


8814 


8820 


8825 


8831 


8837 


8842 


8848 


8864 


8860 




77 


8866 


8871 


8876 


8882 


8887 


8803 


8800 


8904 


8010 


8016 




78 


8921 


8\)27 


8032 


8938 


8943 


8940 


8864 


8960 


8866 


8071 




70 


8976 


8982 


8087 


8903 


8998 


0004 


0000 


0016 


0020 


0025 




80 


9081 


9036 


8042 


0047 


9053 


0068 


0068 


0060 


0074 


0070 


6 


81 


9065 


0000 


0096 


0101 


9106 


0112 


0117 


0122 


0128 


0138 




82 


9138 


0148 


9140 


0164 


9160 


0166 


0170 


0176 


0180 


0186 




83 


9191 


0196 


0201 


0206 


0212 


0217 


0222 


0227 


0232 


0238 




84 


9243 


9248 


9258 


0258 


0263 


0260 


0274 


0270 


0284 


0280 




85 


0204 


0299 


0304 


0600 


0316 


0320 


0325 


0380 


0336 


0840 


6 


86 


0346 


9360 


0855 


0860 


0366 


0370 


0376 


0380 


0885 


0890 




87 


0305 


0400 


0406 


0410 


0416 


0420 


0425 


0430 


0436 


9440 




88 


0446 


0450 


0465 


0460 


0466 


0460 


0474 


0470 


04fli 


0480 




88 


0494 


0400 


0504 


0600 


0618 


0518 


0623 


0628 


0533 


0688 




00 


0642 


0547 


0668 


0657 


0662 


0666 


0571 


8576 


0581 


0686 


6 


01 


OBOO 


0606 


0600 


0606 


0600 


0614 


0610 


0624 


0628 


0633 




92 


0638 


0643 


0647 


0668 


0657 


0661 


0666 


0671 


0675 


0680 




08 


0686 


0680 


0604 


0600 


0703 


0708 


ono 


0717 


0722 


0727 




04 


0731 


0736 


0741 


0745 


0750 


0764 


0760 


0763 


0768 


ons 




06 


0777 


0782 


0786 


0701 


0705 


0800 


0805 


0600 


0814 


0818 


4 


06 


0823 


0827 


0832 


0886 


0841 


9845 


0850 


0654 


0850 


0868 




07 


0868 


0872 


0877 


0881 


0886 


9880 


0804 


0800 


0003 


0008 




08 


0012 


0017 


0021 


0026 


9030 


0934 


0080 


0043 


0048 


0052 




00 


0956 


0061 


0066 


0060 


0074 


9978 


0068 


0087 


0001 


0006 




n 





1 


2 


8 


4 


6 


6 


7 


8 
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TABLE III. — Squares of Numbers. 



n 





1 


2 


3 


4 


6 


6 


7 


8 


9 


Diff. 


1.0 


1.000 


1.020 


1.040 


1.061 


1.082 


1.103 


1.124 


1.146 


1.166 


1.188 


22 


1.1 


1.210 


1.232 


1.264 


1.277 


1.300 


1.323 


1.346 


1.369 


1.392 


1.416 


24 


1.2 


1.440 


1.464 


1.488 


1.613 


1.638 


1.563 


1.588 


1.613 


1.638 


1.664 


26 


1.3 


1.690 


1.716 


1.742 


1.769 


1.796 


1.823 


1.850 


1.877 


1.904 


1.932 


28 


1.4 


1.960 


1.988 


2.016 


2.046 


2.074 


2.103 


2.132 


2.161 


2.190 


2.220 


80 


1.6 


2.260 


2.280 


2.310 


2.341 


2.372 


2.403 


2.434 


2.466 


2.496 


2.628 


32 


1.6 


2.660 


2.692 


2.624 


2.667 


2.690 


2.723 


2.756 


2.789 


2.822 


2.856 


34 


1.7 


2.890 


2.924 


2.968 


2.993 


3.028 


3.063 


3.098 


3.133 


3.168 


3.204 


36 


1.8 


3.240 


3.276 


3.312 


3.349 


3.386 


3.423 


3.460 


3.497 


8.534 


3.572 


38 


1.9 


3.610 


3.648 


3.686 


8.726 


8.764 


3.803 


3.842 


3.881 


3.920 


3.960 


40 


2.0 


4.000 


4.040 


4.060 


4.121 


4.162 


4.203 


4.244 


4.286 


4.326 


4.368 


42 


2.1 


4.410 


4.462 


4.494 


4.537 


4.680 


4.623 


4.666 


4.709 


4.752 


4.796 


44 


2.2 


4.840 


4.884 


4.928 


4.973 


6.018 


5.063 


5.108 


5.153 


6.198 


6.244 


46 


2.3 


6.290 


6.336 


6.382 


6.429 


6.476 


5.523 


6.570 


6.617 


6.664 


5.712 


48 


2.4 


6.760 


6.808 


6.866 


6.906 


6.964 


6.003 


6.052 


6.101 


6.160 


6.200 


60 


2.6 


6.260 


6.300 


6.360 


6.401 


6.462 


fiJSOS 


6.664 


6.606 


6.666 


6.708 


62 


2.6 


6.760 


6.812 


6.864 


6.917 


6.970 


7.023 


7.076 


7.129 


7.182 


7.236 


64 


2.7 


7.290 


7.344 


7.398 


7.463 


7.608 


7.668 


7.618 


7.673 


7.728 


7.784 


56 


2.8 


7.840 


7.896 


7.962 


8.009 


8.066 


8.123 


8.180 


8.237 


8.294 


8.352 


68 


2.9 


8.410 


8.468 


8.626 


8.686 


8.644 


8.703 


8.762 


8.821 


8.880 


8.940 


60 


3.0 


9.000 


9.060 


9.120 


9.181 


9.242 


9.303 


9.364 


9.426 


9.486 


9.648 


62 


3.1 


9.610 


9.672 


9.734 


9.797 


9.860 


9.923 


9.986 


10.06 


10.11 


10.18 


6 


3.2 


10.24 


10.30 


10.37 


10.43 


10.60 


10.56 


10.63 


10.69 


10.76 


10.82 


7 


3.3 


10.89 


10.96 


11.02 


11.09 


11.16 


11.22 


11.29 


11.36 


11.42 


11.49 


7 


8.4 


11.66 


11.63 


11.70 


11.76 


11.83 


11.90 


11.97 


12.04 


12.11 


12.18 


7 


3.6 


12.26 


12.32 


12.39 


12.46 


12.53 


12.60 


12.67 


12.74 


12.82 


12.89 


7 


3.6 


12.96 


13.03 


13.10 


13.18 


13.26 


13.32 


13.40 


13.47 


13.54 


13.62 


7 


3.7 


13.69 


13.76 


13.84 


13.91 


13.99 


14.06 


14.14 


14.21 


14.29 


14.36 


8 


3.8 


14.44 


14.62 


14.69 


14.67 


14.76 


14.82 


14.90 


14.98 


15.06 


15.13 


8 


3.9 


16.21 


16.29 


16.37 


16.44 


15.52 


16.60 


16.68 


16.76 


16.84 


15.92 


8 


4.0 


16.00 


16.06 


16.16 


16.24 


16.32 


16.40 


16.48 


16.66 


16.66 


16.73 


8 


4.1 


16.81 


16.89 


16.97 


17.06 


17.14 


17.22 


17.31 


17.39 


17.47 


17.56 


8 


4.2 


17.64 


17.72 


17.81 


17.89 


17.98 


18.06 


18.16 


18.23 


18.32 


18.40 


9 


4.3 


18.49 


18.68 


18.66 


18.76 


18.84 


18.92 


19.01 


19.10 


19.18 


19.27 


9 


4.4 


19.38 


19.46 


19.64 


19.62 


19.71 


19.80 


19.89 


19.98 


20.07 


20.16 


9 


4.6 


20.26 


20.34 


20.43 


20.52 


20.61 


20.70 


20.79 


20.88 


20.96 


21.07 


9 


4.6 


21.16 


21.26 


21.34 


21.44 


21.53 


21.62 


21.72 


21.81 


21.90 


22.00 


9 


4.7 


22.09 


22.18 


22.28 


22.37 


22.47 


22.66 


22.66 


22.76 


22.86 


22.94 


10 


4.8 


23.04 


23.14 


23.23 


23.33 


23.43 


23.52 


23.62 


23.72 


23.81 


23.91 


10 


4.9 


24.01 


24.11 


24.21 


24.30 


24.40 


24.60 


24.60 


24.70 


24.80 


24.90 


10 


6.0 


26.00 


26.10 


26.20 


26.30 


26.40 


25JS0 


26.60 


26.70 


26.81 


26.91 


10 


6.1 


26.01 


26.11 


26.21 


26.32 


26.42 


26J}2 


26.68 


26.73 


26.83 


26.94 


10 


6.2 


27.04 


27.14 


27.26 


27.36 


27.46 


27.66 


27.67 


27.77 


27.88 


27.98 


11 


6.3 


28.09 


28.20 


28.30 


28.41 


28.52 


28.62 


28.78 


28.84 


28.94 


29.06 


11 


6.4 


29.16 


29J27 


29.38 


29.48 


29.69 


29.70 


29.81 


29.92 


30.03 


30.14 


11 


w 





1 


2 


3 


4 


6 


6 


7 


8 


9 


Diff. 
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TABLE III.— Squares of Numbers. 



n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


Dili. 


6.6 


30.26 


30.36 


30.47 


30.68 


30.69 


30.80 


30.91 


81.02 


31.14 


31.26 


11 


6.6 


31.36 


31.47 


31.68 


31.70 


31.81 


31.92 


32.04 


32.16 


32.26 


32.38 


11 


6.7 


32.49 


32.60 


32.72 


82.83 


32.96 


33.06 


33.18 


33.29 


33.41 


33.52 


12 


6.8 


33.64 


33.76 


33.87 


33.99 


34.11 


34.22 


o4.«>4 


34.46 


34.67 


34.69 


12 


6.9 


34.81 


34.93 


36.06 


35.16 


36.28 


36.40 


36.62 


36.64 


36.76 


36.88 


12 


6.0 


36.00 


36.12 


36.24 


36.36 


36.48 


36.60 


36.72 


36.84 


36.97 


37.09 


12 


6.1 


37.21 


37.33 


37.46 


37.68 


37.70 


37.82 


37.96 


38.07 


38.19 


38.32 


12 


6.2 


38.44 


38.66 


38.69 


38.81 


38.94 


39.06 


39.19 


39.31 


89.44 


39.56 


13 


6.3 


39.69 


39.82 


39.94 


40.07 


40.20 


40.32 


40.46 


40.68 


40.70 


40.83 


13 


6.4 


40.96 


41.09 


41.22 


41.34 


41.47 


41.60 


41.73 


41.86 


41.99 


42.12 


13 


6.6 


42.26 


42.38 


42.61 


42.64 


42.77 


42.90 


43.03 


48.16 


43.30 


43.43 


13 


6.6 


43.66 


43.69 


43.82 


43.96 


44.09 


44.22 


44.36 


44.49 


44.62 


44.76 


13 


6.7 


44.89 


46.02 


46.16 


46.29 


46.43 


46JS6 


46.70 


46.83 


46.97 


46.10 


14 


6.8 


46.24 


46.38 


46.61 


46.66 


46.79 


46.92 


47.06 


47.20 


47.38 


47.47 


14 


6.9 


47.61 


47.76 


47%89 


48.02 


48.16 


48.30 


48.44 


48.68 


48.72 


48.86 


14 


7.0 


49.00 


49.14 


49.28 


49.42 


49.66 


49.70 


49.84 


49.96 


60.13 


60.27 


14 


7.1 


60.41 


60.66 


60.69 


60.84 


60.98 


61.12 


61.27 


61.41 


61.66 


61.70 


14 


7.2 


61.84 


61.98 


62.13 


62.27 


62.42 


52.66 


62.71 


62.86 


63.00 


63.14 


16 


7.8 


63.29 


63.44 


63ii8 


63.73 


63.88 


64.02 


64.17 


64.32 


64.46 


64.61 


16 


7.4 


64.76 


64.91 


66.06 


66.20 


66.36 


66.60 


66.66 


66.80 


66.96 


66.10 


16 


7.6 


66.26 


66.40 


66.66 


66.70 


66.86 


67.00 


67.16 


67.30 


67.46 


67.61 


16 


7.6 


67.76 


67.91 


68.06 


68.22 


68.37 


68.62 


68.68 


68.83 


68.98 


69.14 


16 


7.7 


69.29 


69.44 


69.60 


60.76 


69.91 


60.06 


60.22 


60.37 


60.63 


60.68 


16 


7.8 


60.84 


61.00 


61.16 


61.31 


61.47 


61.62 


61.78 


61.94 


62.09 


62.26 


16 


7.9 


62.41 


62.67 


62.73 


62.88 


63.04 


63.20 


63.36 


63.62 


63.68 


63.84 


16 


8.0 


64.00 


64.16 


64.32 


64.48 


64.64 


64.80 


64.96 


66.12 


66.29 


66.46 


16 


8.1 


66.61 


66.77 


66.93 


66.10 


66.26 


66.42 


66.60 


66 76 


66.91 


67.08 


16 


8.2 


67.24 


67.40 


67.67 


67.73 


67.90 


68.06 


68.23 


68.39 


68.66 


68.72 


17 


8.3 


68.89 


69.06 


69.22 


69.39 


69JS6 


69.72 


69.89 


70.06 


70.22 


70.39 


17 


8.4 


70.66 


70.73 


70.90 


71.06 


71.23 


71.40 


71.67 


71.74 


71.91 


72.08 


17 


8.5 


72.26 


72.42 


72J» 


72.76 


72.93 


73.10 


73i!7 


73.44 


73.62 


73.79 


17 


8.6 


73.96 


74.13 


74.30 


74.48 


74.66 


74.82 


76.00 


76.17 


76.34 


76.62 


17 


8.7 


76.69 


76.86 


76.04 


76.21 


76.39 


76.66 


76.74 


76.91 


77.09 


77.26 


18 


8.8 


77.44 


77.62 


77.79 


77.97 


78.16 


78.32 


78.60 


78.68 


78.86 


79.03 


18 


8.9 


79.21 


79.89 


79.67 


79.74 


79.92 


80.10 


80.28 


80.46 


80.64 


80.82 


18 


9.0 


81.00 


81.18 


81.36 


81.64 


81.72 


81.90 


82.08 


82.26 


82.46 


82.63 


18 


9.1 


82.81 


82.99 


83.17 


83.36 


83.64 


83.72 


83.91 


84.09 


84.27 


84.46 


18 


9.2 


84.64 


84.82 


86.01 


86.19 


86.38 


85.56 


86.75 


86.93 


86.12 


86.30 


19 


9.8 


86.49 


86.68 


86.86 


87.06 


87.24 


87.42 


87.61 


87.80 


87.98 


88.17 


19 


9.4 


88.36 


88.66 


88.74 


884»2 


89.11 


89.30 


89.49 


89.68 


89.87 


90.06 


19 


9.6 


90.26 


90.44 


90.63 


90.82 


91.01 


91.20 


91.39 


91.68 


91.78 


91.97 


19 


9.6 


92.16 


92.36 


92.64 


92.74 


92.93 


93.12 


93.32 


93.61 


93.70 


93.90 


19 


9.7 


94.09 


94.28 


94.48 


94.67 


94.87 


96.06 


96.26 


96.46 


96.66 


96.84 


20 


9.8 


96.04 


96.24 


96.43 


96.63 


96.83 


97.02 


97.22 


97.42 


97.61 


97.81 


20 


9.9 


96.01 


96.21 


98.41 


96.60 


96.80 


99.00 


99.20 


99.40 


99.60 


99.80 


20 


n 





1 


2 


8 


4 


6 


6 


7 


8 


9 


Diff. 



i . A 



i'- 



i» 






